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Notation 



- Starting from Chap. 2 natural units with h = 1 and c = 1 are used. 

- The Minkowski metric signature is (1,-1,-1,-1). 

- The Greek indices are spacetime indices running from to 3. The ones 
at the beginning of the Greek alphabet like a, (3, 7 or 5 label quantities 
referring to Minkowski spacetime while the others in the middle of the 
Greek alphabet like A, fj,, v or p label quantities referring to a generic 
curved spacetime. 

- The Latin indices at the beginning of the Latin alphabet like a, b or 
c are spinor indices running from 1 to 4. The Latin index i is a space 
index running from 1 to 3. Finally, the Latin indices at the end of the 
Latin alphabet like u or v are coordinate indices standing for x, y or 
z. 

- The index j embodies all the quantum numbers of a relativistic elec- 
tron (positron) in the presence of a constant and uniform magnetic 
field in Minkowski spacetime. The index J embodies the correspond- 
ing quantum numbers but in curved spacetime. Instead, the index j 
indicates a set of quantum numbers to be determined. 

- The sans serif letters like x or y embody the four general coordinates 
of a spacetime fourpoint. 

- The subscript "T" refers to quantities concerning the physical situation 
in which the magnetic field is always directed along the z axis. 

- The subscript "/" refers to quantities concerning the physical situation 
in which the magnetic field lies in the y—z plane with a nonzero y 
component. 

- The subscripts "_L" and "||" refer to quantities perpendicular and par- 
allel to the magnetic field respectively. 

- The superscripts "(ann)", "(syn)" and "(dir)" refer to photon pro- 
duction mechanisms and stand for "annihilation" , "synchrotron" and 
"direct" respectively. 
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Notation 



- The superscripts "lin" and "exp" refer to the magnetic field time vari- 
ation and stand for "linear" and "exponential" respectively. 

- The superscript "(R)" refers to quantities calculated in Rindler space- 
time. 

- The primed classical and quantum fields, the primed propagators and 
the primed 5-matrix elements refer to the physical situation in which 
the magnetic field is always directed along the z axis. 

- The capital calligraphic letters like V or H refer to one-particle electron 
and positron quantum operators. 

- The variant capital calligraphic letters like «Sf or Jt? refer to field 
density quantities. 

- The "Euler" font letters like A, F or e refer to quantities concerning 
the photon radiation field. 
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During the years Quantum Electrodynamics (QED) has received a large 
amount of experimental confirmations and today it is considered the most 
reliable microscopic physical theory we have. It is enough to think about the 
excellent agreement between the predicted [1] and the measured [2] anoma- 
lous magnetic moment of the electron. The main subject of this thesis 
concerns a very interesting and fascinating branch of QED: Quantum Elec- 
trodynamics in the presence of External Fields (QEDEF) OHj. As the name 
itself suggests, QEDEF concerns those quantum electromagnetic processes 
happening when a classical (nonquantized) field is present. The external 
fields are typically electromagnetic fields and they have two main general 
characteristics: 

1. they are produced by sources not belonging to the system under study 
and their spatiotemporal evolution is assigned; 

2. they are so intense that: 

2a. quantizing them would be useless; 

2b. trying to compute their effects by means of standard perturbative 
techniques would be in most cases conceptually wrong. 

Now, many processes studied in pure QED can also be analyzed in the 
presence of external fields such as particle scattering and so on. But there 
is a process that is typical of QEDEF: the production of particles from vac- 
uum. The reason is understandable because pure QED describes a closed 
system whose energy is conserved and then particles can not be created from 
vacuum. Instead, in QEDEF just the external fields can supply the energy 
necessary to create electron-positron pairs or photons directly from vacuum. 
The seminal work about this subject is that by Schwinger [5] in which he 
shows that in the presence of a constant and uniform electric field, real pair^] 
can be produced directly from vacuum. The production probability results 
significantly different from zero only for electric field strengths much larger 
than E cr = m 2 c 3 /(he) = 1.3 x 10 16 V/cm where m and — e < are the mass 

1 From now on, when it is not specified, it is understood that "pair" ("pairs") stands 
for "electron-positron pair" ("electron-positron pairs"). 
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and the charge of the electron respectively. This value of E cr represents a 
sort of benchmark over which the effects of an external electric field become 
important. Nevertheless, we will see in the next Chapter that E cr is much 
larger than the electric fields that today can be produced in terrestrial lab- 
oratories and the experimental confirmations concerning QEDEF are less 
numerous than those concerning pure QED. 

The expression of E cr can be seen qualitatively as that of the electric 
field strength whose energy in a volume with typical length of the order of 
the Compton length X = h/ (mc) is large enough to produce an electron and 
a positron at rest. By using the same energetic argument it can be seen that 
a magnetic field strength stronger than B cr = rr?(?l(he) = 4.4 x 10 13 gauss 
is capable to "break" the vacuum and create a pairo The value of B cr is also 
much larger than that of the strongest, steady magnetic field ever produced 
in a terrestrial laboratory which is of the order of 10 5 gauss [6]. Instead, we 
will see in the next Chapter that there are various indirect evidences that 
around astrophysical compact objects (neutron stars, black holes) magnetic 
fields much stronger than B cr are present. Actually, this is still not a good 
enough reason to start studying QED in the presence of such strong mag- 
netic fields if an (at least) indirect experimental effect of their presence is 
not at hand. Now, one of the most intriguing and mysterious astrophysical 
phenomena seems to be originated just around neutron stars or black holes 
surrounded by an accretion disk. I am referring to the so-called Gamma-Ray 
Bursts (GRBs) that are huge pulses of soft gamma-rays that our satellites 
register on average once a day. The exact GRBs production mechanism is 
still not completely well understood. Actually, following the widely accepted 
fireball model of GRBs, the photons making a GRB are produced, through 
standard electromagnetic mechanisms such as synchrotron emission, by a 
"fireball" in turn made mostly of electrons, positrons and photons them- 
selves. Nevertheless, it is not still clear how the fireball itself is generated. 
On the other hand, as we will see in the next Chapter, it seems almost sure 
that the fireball is produced near a neutron star or a black hole surrounded 
by an accretion disk and that the huge magnetic fields that are there play 
a fundamental role in this process. In this respect, trying to find a mi- 
croscopic underlying mechanism responsible of the formation of the fireball 
represents for a theoretical physicist a very challenging reason to study those 
electromagnetic processes that can be primed in the presence of such strong 
magnetic fields. In fact, in this way, some general experimental features of 
GRBs could be explained and interpreted. 

In this thesis I will focus in particular on particle production processes 
in the presence of strong (i.e., much larger than B cr ), uniform and slowly- 

2 It must be pointed out that, although from an energetic point of view pairs can be 
created in the presence of a constant and uniform magnetic field with strength larger than 
Bc r , it is impossible, as Schwinger also showed in [5], from a dynamical point of view 
essentially because the Lorentz force does not do any work. 
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varying magnetic fields. The theoretical justifications of these assumptions 
about the structure of the magnetic fields considered in the present work 
will be given in the following Chapters that are organized as follows. Chap. 
1 is divided into two independent parts: in the first one I briefly review some 
developments and theoretical predictions about the production of pairs in 
the presence of external electromagnetic fields; instead, in the second one 
I give a brief introduction about the phenomenological characteristics of 
GRBs and about the fireball model of GRBs. In Chap. 2 I just outline the 
following three "textbook" subjects: the motion of a relativistic charged par- 
ticle in the presence of a constant and uniform magnetic field, the adiabatic 
perturbation theory in quantum mechanics and the quantum field theory in 
curved spacetime. Even if these subjects (at least the first two) are well- 
known I treat them for the sake of completeness because they represent the 
theoretical tools I have used to derive my own results. These last are pre- 
sented in Chaps. 3-5. In particular, in Chap. 3 I study the production of 
an electron-positron pair from vacuum in the presence of a strong, uniform 
and slowly-varying magnetic field. In particular two different magnetic field 
time evolutions are considered. Since, as I have said, a possible application 
of these calculations concerns the mechanisms accounting for the production 
of a GRB, in Chap. 4 I study some processes through that the electrons 
and the positrons already created in a strong, uniform and slowly-rotating 
magnetic field can produce photons. In particular, I calculate the energy 
spectra of photons produced through pair annihilation and as synchrotron 
radiation. Also, by using the so-called effective Lagrangean method I obtain 
the energy spectrum of the photons produced directly from vacuum in the 
presence of a strong, uniform and slowly-rotating magnetic field. Finally, in 
the astrophysical scenario I have just sketched, it is interesting to examine 
what the role of the compact object's gravitational field is. I devote the 
last Chapter to this subject by analyzing, in the context of quantum field 
theory in curved spacetime, how the presence of a weak and of a strong grav- 
itational field modifies some results on the production of electron-positron 
pairs already obtained in Chap. 3. 
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Chapter 1 

Physical background and 
motivations 

As I have sketched in the Introduction, despite the possible application to 
the study of GRBs, this work is principally theoretical. As a consequence, 
in this Chapter I want to firstly describe the scientific "landscape" where 
this work is included with particular reference to pair production in the 
presence of strong electromagnetic fields (Sect. Il.ip . Also, in order to make 
the cited astrophysical application clear I give a brief description of the 
general experimental features of GRBs, of the physical scenario where they 
are supposed to be originated and of the widely accepted fireball model of 
GRBs (Sect. OP . 

1.1 Pair production in the presence of strong elec- 
tromagnetic fields: a short review 

This work will concern mainly the theoretical study of the production of 
pairs in the presence of a strong, uniform and time- varying magnetic field in 
the framework of QEDEF. Then, it is worth giving a brief review about, in 
general, the production of pairs in the presence of external electromagnetic 
fields. Actually, books [3j 0] and conferences [7] have been devoted to this 
subject and then my review will be unavoidably incomplete. 

As I have said in the Introduction, the possibility that real pairs can 
be produced from vacuum was first investigated by Schwinger in [5J. In 
particular, he discussed the pair production process in the presence of an 
external constant and uniform electric field. The calculations are rather 
difficult but the production mechanism, sometimes called Schwinger mecha- 
nism, is qualitatively easy to be understood. In fact, it is known that virtual 
pairs are spontaneously produced and then annihilated as "vacuum fluctu- 
ations". The electron and the positron making a virtual pair "live" at a 
very short distance between each other of the order of the Compton length 
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A = 3.9 x 10 cm. Nevertheless, if the external electric field strength 
is larger than E cr = 1.3 x 10 16 V/cm, then the electron and the positron 
forming the virtual pair can be separated becoming real particles. In terms 
of the Dirac picture of the quantum vacuum, in the presence of an electric 
field E with strength of order of E cr the electron levels are so distorted that 
an electron with negative energy can "tunnel" to the positive energy levels 
leaving a positively charged "hole" (see Fig. II. 1|) . The interpretation of the 
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Figure 1.1: Tunneling of an electron from a negative to a positive energy 
level in the presence of a strong electric field E. 

pair production mechanism in terms of the tunnel effect is confirmed by the 
typically nonperturbative expression of the production probability. In fact, if 
the external field strength E is much lower than E cr , it results proportional 
to exp(— -kE ct /E) [5]. The same previous qualitative interpretation but also 
rigorous calculations show that the production of pairs from vacuum is im- 
possible in the presence of a pure, constant and uniform magnetic field. In 
fact, unlike the electric force, the Lorentz force can not do the work neces- 
sary to "separate" the electron and the positron forming a virtual pair and 
to make them real particles. 

The extraordinarily high value of E cr makes impossible to test in a ter- 
restrial laboratory the pair production Schwinger mechanism. A further 
difficulty is represented by the fact that in the Schwinger formalism the 
electric field is assumed to be constant and uniform while from an experi- 
mental point of view it is easier to produce strong alternating electric fields 
(lasers). A wide literature has been devoted to the theoretical study of 
pair production in the presence of uniform, alternating fields. The first 
work was that by Brezin and Itzykson [8]. In that paper it is found that 
if r?o is the electric field rotational frequency then the pair production is 
significantly different from zero if the electric field strength E is such that 
E > mQ,Qc/[esinh[hQ,Q/(4:mc 2 )]]. Even using the X-ray lasers actually at 
our disposal, the effect is still too tiny to be observed. In fact, the peak 
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value of the electric field strength at the focus of an X-ray laser necessary 
to measure some effects of pair creation is of the order of 10 15 V/cm, that is 
five orders of magnitude larger than those of the now available X-ray lasers 

Due to this unavoidable "inadequacy" of the experimental resources, var- 
ious electromagnetic field configurations have been analyzed by theoretical 
physicists in order to find the most efficient one from the point of view of 
particle production. For example, in [HI [12] the production of pairs in the 
presence of a uniform and time-dependent (but nonoscillating) electric field 
is considered. The authors use the effective Lagrangean method [similar to 
that used by Schwinger in [5]] to calculate approximatively the number of 
pairs per unit volume and unit time created from vacuum as twice the imagi- 
nary part of the effective Lagrangean density of the system |131I14|. Instead, 
in [T5] the same technique is used to estimate numerically the production of 
electrons and positrons in the presence of electric and magnetic fields cou- 
pled in various configurations: an alternating electric field superimposed to 
a uniform magnetic field or an alternating electric field superimposed to an 
also alternating magnetic field. Finally, in another class of papers the pos- 
sible multiple pair production in the electromagnetic field of two colliding 
heavy ions is studied [16\ [T7] but till now the experimental evidence of this 
mechanism of particle production is inconclusive and controversial. 

Finally, a different theoretical method has been used to predict, in gen- 
eral, particle production in the presence of external time-dependent fields: 
the Bogoliubov transformation. In order to describe the physical content of 
the Bogoliubov transformation, I assume, as usual in quantum field theory, 
to work in the Heisenberg picture. Now, when a time-dependent external 
field is present the second quantized number operators associated with a 
generic quantum field change nontrivially with time. It may happen that 
even if the time-independent state of the quantum field system is the vacuum 
state of the number operators calculated at t — > — oo it is not the vacuum 
state of the number operators at t — ► ooQ In this case, one just concludes 
that while in the far past no particles were present, they are there in the 
far future and then that the external time-dependent field induced their 
production directly from vacuum. As it is evident from this qualitative de- 
scription, the use of the Bogoliubov transformation technique is not limited 
to the study of electron-positron pairs production in the presence of time- 
dependent electromagnetic fields. In fact, it has also been used, for example, 
in [181 QjJJ [201 ED (221 E3] to calculate the production rate of different kinds of 
particles in the presence of time-dependent cosmological gravitational fields. 



It can be shown that these two classes of operators are connected by means of a linear 
transformation called Bogoliubov transformation. 
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1.2 GRBs: a brief introduction 

In this Section I want to give a brief introduction concerning GRBs, their 
general features and so on. The literature about GRBs is really endless but 
the present introduction results mostly from my reading of some nonspecial- 
istic papers (21 [251 [261 EZ] and some review papers [281 ESI Efl M\ ■ In the 
following I will not quote the previous references anymore but only other 
ones where I found specific information. 

I. 2.1 General experimental characteristics of GRBs 

The first GRB was detected by chance on July 2, 1967 by a military satel- 
lite watching for nuclear tests in space. The experimental data concerning 
that burst were published much later and, actually, the scientific community 
started studying systematically GRBs only after the launch of the Comp- 
ton Gamma-Ray Observatory (CGRO) and its detector BATSE (Burst and 
Transient Experiment). BATSE made a crucial contribution in establishing 
the distance scale of GRBs. In fact, it detected GRBs from all directions 
in the sky with a completely isotropic distribution (see Fig. 11.2)) strongly 
supporting the idea that GRBs are originated at cosmological (~ 10 28 cm) 
distances from usJl The cosmological origin of GRBs has tremendous conse- 
quences about the total energy carried by them. In fact, the detected GRBs 
fluences range typically from 10~ 7 erg/cm 2 to 10 -4 erg/cm 2 [30J, then, as- 
suming an isotropic emission, the total energy carried by a GRB is of the 
order of 10 50 erg-10 53 erg. For this reason, GRBs have been characterized 
as the brightest explosions in the Universe after Big Bang. 

Now, thousands of GRBs have been detected by high-energy astrophysics 
satellites and, as we will see in the next Paragraph, a feature that makes 
them hard to be understood is just their wide variety. For example, Fig. 

II. 31 shows very different GRBs time profiles, i. e. the number of photons 
detected per unit time as a function of time. In this respect, GRBs typically 
last from 0.01 s to 100 s and are generally classified into two large categories: 
short bursts lasting less than 2 s and long bursts the others. We know very 
little about short bursts. In general, their energy spectra are harder than 
those of the long GRBs and they are supposed to be originated nearer with 
respect to long bursts [32] ■ Instead, our knowledge about long bursts has 
been made deeper thanks to the discovery of their X-ray afterglows by means 
of the Italian-Dutch satellite Beppo-SAX. Actually, the data from Beppo- 
SAX showed that long GRBs are followed by a multiwavelength "afterglow" 
made of X-ray photons, ultraviolet photons down to radio photons. On 
the one hand, the optical observations confirmed the cosmological origin of 

2 By contrast, the isotropic distribution of GRBs rules out the possibility that, as it 
was firstly believed, GRBs are produced in our galaxy because it is not spherical and we 
do not occupy a privileged position in it. 
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1707 BATSE Gamma-Ray Bursts 



+90 




Galactic Coordinates 



Figure 1.2: Spatial distribution in the sky of 1707 GRBs detected by BATSE. 

GRBs. On the other hand, the presence of iron lines in the afterglows was 
the first indication of a connection, that has become stronger and stronger, 
between GRBs and supernovae explosions. In fact, it is known that during 
these explosions large amount of iron atoms are synthesized and ejected. 

I want to conclude this Paragraph by quoting three other general features 
of GRBs that give quite striking conditions mostly on the mechanism(s) 
responsible of their production: 

• the energy spectrum of GRBs is not thermal; 

• the time variability of the GRBs signals is observed down to time scales 
of less than 10 ms (see Fig. II. 3D : 

• the prompt gamma-ray emission of many GRBs shows a high linear 
polarization degree [33] . 

Concerning the first point, Fig. 11.41 shows two typical GRBs energy spectra. 
In general, called n(u) the number of photons per unit area and unit time 
as a function of the photon energy uo (indicated as E in Fig. I1.4p . a typical 
experimental GRB energy spectrum is well fitted by the following piecewise 
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Figure 1.3: Various GRBs time profiles. 
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with a ~ 1 and (3 ~ 2-3 and with u)b a break energy. Eq. (jl.2.1[) represents 
essentially a double power-law smoothly joined at w = by a decreas- 
ing exponential. The break energies uj}, typically lie between 0.1 MeV and 
0.3 MeV, even if bursts with > 1 MeV have been detected [3^| 155]. 



1.2.2 The fireball model of GRBs 

In this Paragraph I want to show how the previous experimental evidences 
led astrophysicists to work out the so-called fireball model [36l E7] that is 
the most widely accepted model of long GRBsH 

In the fireball model the photons forming a GRB are thought to be 
emitted by a "fireball" , a plasma made essentially of electrons, positrons and 



3 Actually, due to the great variety of GRBs, many models of GRBs have been pro- 
posed during the years describing GRBs as jets from pulsars [38] or as bursts emitted by 
cannonballs in turn originated during supernovae explosions [35] etc.... 
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E, keV E, keV 

Figure 1.4: Two typical GRBs energy spectra. 

photons, that expands high relativistically [10] @ In general, a direct emission 
of the photons forming a GRB would be inconsistent with the experimental 
evidence about the nonthermal spectrum of GRBs. In fact, as we have seen, 
during a burst 10 50 erg-10 53 erg are released on average in less than 100 s. 
Also, the fact that GRBs show in the soft-gamma region variabilities at time 
scales of the order of 1 ms implies, by using the causality limit, that their 
source should have a linear length of the order of 1 ms x c = 3.0 x 10 7 cm. 
In these extreme conditions, even if initially only photons were produced 
then their energy (< 1 MeV) and their density would be high enough to 
create electron-positron pairs through the reaction 77 —* e~e + . In turn, the 
strong electromagnetic interactions among the electrons, the positrons and 
the photons would allow these last to escape only after so much time that 
they would thermalize. 

As I have only mentioned at the beginning of this Paragraph, the fire- 
ball model is intended to describe only long GRBs. In fact, the feature that 
mostly has contributed to its success is just the prediction of the existence of 
the afterglow. Actually, the fireball model gives quite naturally a mechanism 
for producing the true burst and the afterglow (see Fig. I1.5jl . In general, 
in this model a GRB and its afterglow are produced at the expense of the 

4 One of the strength points of the fireball model is its almost complete independence 
from the nature of the central engine that produces the fireball (see Sect. I1.2.3[l and from 
the mechanism(s) giving rise to the fireball itself. 
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Figure 1.5: Schematic representation of the mechanisms producing a long 
GRB and its afterglow in the framework of the fireball model. The left part 
of the Figure and the inset refer to two possible progenitors of the fireball 
that are discussed in Par. 11.2.31 

kinetic energy of the fireball. In particular, the true burst is produced as 
a consequence of the internal shocks among different layers of the fireball 
travelling at different speeds |41|. Of course, although these shocks may 
last hours, the resulting burst can last only few seconds because the fire- 
ball expands high relativistically with a Lorentz factor also more than 100. 
Also, the fact that the electrons and the positrons in the fireball are very 
energetic explains why during a burst photons can be emitted with energies 
of the order of 1 GeV or more and why the time profile of a burst can also 
show structures at the millisecond scale. Instead, the afterglow of a GRB is 
produced as a consequence of the shocks between the external layer of the 
fireball and the surrounding medium |42j. The two mechanisms producing 
the true burst and the afterglow are independent from each other but, in 
general, the internal shocks and the expansion slow down the fireball before 
it encounters the interstellar medium. In this way the photons forming the 
afterglow result, in general, to be less energetic than those forming the true 
burst and the afterglow itself lasts much longer than the true burst H 



Nevertheless, some GRBs have been detected having an afterglow "harder" than the 
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If we want to go a bit more deeply into the electromagnetic processes re- 
sponsible of the production of the photons forming a GRB and its afterglow 
we have to take into account the indications suggested by the experimental 
data (see the previous Paragraph). Firstly, as I have just mentioned before, 
the photons forming a GRB can not be those present in the fireball because 
their energetic spectrum would be thermal. Moreover, all the experimental 
indications support the idea that the main mechanisms responsible of the 
production of the photons forming a GRB are the inverse Compton effect 
and, mostly, the synchrotron emission. The inverse Compton effect, that 
is the scattering between an energetic electron (positron) and a soft pho- 
ton resulting in a less energetic electron (positron) and a harder photon, is 
important in explaining the high-energy (of the order of GeV) part of the 
photon spectrum of a GRB [33]. Instead, the photons forming the other less 
energetic parts of GRBs photon spectra are thought to be produced as syn- 
chrotron radiation by the shock-accelerated electrons and positrons forming 
the fireball. Firstly, this photon production mechanism explains the pres- 
ence of the break energy uj, in the energy spectra (see Fig. I1.4H and the 
power-law behaviour below uj^ [see Eq. (|1.2.1|) ]. Most important, the recent 
discovery of the highly linear polarization degree [II = (80 ± 20) %] of the 
gamma-ray emission of the GRB detected on December, 6 2002 [33J strongly 
supports the idea that those photons are emitted through synchrotron radi- 
ation just because synchrotron radiation is theoretically known to be high 
linearly polarized [44 \ H5"] . Obviously, in order that synchrotron radiation 
can be emitted, a magnetic field must be present in the emission region. 
The high polarization degree itself suggests that the magnetic field in the 
emission region should be nearly uniform and strong enough to "order" the 
motion of the emitting particles in the fireball |33|. 146] . This is in turn a 
clear indication that the true GRB is generated in regions where the mag- 
netic field is that produced by a macroscopic object that, as we will see in 
the next Paragraph, is identified with the central engine that powers the fire- 
ball |47j . By contrast, the polarization degree of the afterglow photons has 
been measured to be of the order of 10 %. This much lower value indicates 
that the afterglow is produced in regions where the magnetic field can also 
have a "turbulent" spatial distribution and can be produced by the shock- 
accelerated electrons and positrons themselves in the fireball [HJ SSI EQ] . 

1.2.3 Central engine 

As I have said before, the fireball model is quite independent from the nature 
of the source, the central engine that primes the formation of the fireball. 
Obviously, there is no direct evidence of what the central engine is and, 
for this reason, many candidates have been proposed [51]. In any case, 
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the fact that during a GRB a so large amount of energy is ejected in a so 
relatively short time leads to believe that their sources should be systems 
storing very much energy and that a sudden, explosive event should be 
at the origin of the burst. Two kinds of astrophysical systems have been 
proposed as the possible progenitors powering GRBs (see Fig. I1.5D : on the 
one hand, binary merging systems such as neutron star /neutron star systems 
or neutron star /black hole systems [52\ [3D] and, on the other hand, very 
massive (more than 20 solar masses) stars with a collapsing core, called failed 
supernovae or collapsars [53jH The result of the merging in the first case and 
of the collapse in the second case is the same: a massive rotating black hole 
surrounded by an accretion disk0 Clearly, the observation that many GRBs 
have been localized in star forming regions is an evidence against the first 
kind of progenitors. In fact, the merging of two compact objects comes even 
billions of years later than their formation, then the composite system has 
enough time to drift far from its birth location. Also, the fact that iron (but 
also silicon, sulfur and so on) lines have been detected in GRBs afterglows 
corroborates the hypothesis of the supernova explosion as the "dramatic" 
event producing the fireball because during supernovae explosions such kinds 
of atoms are synthesized and ejected into the interstellar medium. Finally, 
there are now also many direct evidences that GRBs are detected where also 
a supernova is "seen" (55] [56l E7J [58] . It is worth noting that all the previous 
indications against the binary-merging model are obtained from afterglow 
measurements. For this reason, since no afterglows have been so far detected 
for short GRBs, the hypothesis that the binary-merging model can describe 
the formation of short GRBs has been carried out |59[ [60] (since, from now 
on, I will not deal with short GRBs anymore but only with the long ones, I 
will omit the adjective "long" that, actually, will be understood). 

It is clear that a massive rotating black hole surrounded by an accretion 
disk is a huge container of rotational and gravitational energy. But, how this 
energy can be "extracted" from the black hole-accretion disk system to form 
a fireball and then a GRB? Also concerning this subject, there are no direct 
observations and many extraction mechanisms have been proposed. I quote, 
for example, the possibility that the electric fields around charged black holes 
"break" the vacuum producing large amounts of pairs [61]. Another mech- 
anism proposed in the framework of the collapsar model involves the anni- 
hilation into electron-positron pairs of the numerous neutrino-antineutrino 



The adjective "failed" refers to the fact that initially it was supposed that such su- 
pernovae explosions did not eject outside the envelope of the star. 

7 A variant of the collapsar model has been proposed in [54]. In this model, called 
supranova model, the supernova explosion gives rise to a supermassive, rapidly rotating 
neutron star that, in turn, collapses producing a black hole and the fireball. With respect 
to the collapsar model, in the supranova model the fireball is generated much later than 
the supernova explosion (approximatively a week later) in a "cleaner" environment poor 
of heavy baryons. In this way, the energy transfer to the GRB results much more efficient. 
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pairs that are produced during the supernova explosion but quantitative es- 
timates show this mechanism is not efficient enough. Another popular mech- 
anism that has been invoked to explain the formation of the fireball is the 
Blandford-Znajek mechanism [62, 163] [M]. This is a very complicated mech- 
anism |65] but it essentially works as a dynamo that transforms rotational 
energy into electromagnetic energy|f| Now, as in a dynamo, a magnetic field 
is necessary in order that this transformation can happen. In addition, in the 
Blandford-Znajek mechanism a huge, overcritic (3> B cr = 4.4 x 10 13 gauss) 
magnetic field is needed. In fact, the energy extraction rate is proportional 
to the square of the magnetic field strength and theoretical calculations es- 
timate that in order to extract 10 53 erg in less than 1000 s, magnetic field 
strengths of the order of 10 15 gauss are needed [63] . Moreover, the possible 
presence of such strong magnetic field around massive rotating black holes 
surrounded by an accretion disk is also confirmed by numerical simulations 
[see [60] and Ref. [9] in [63]]. In this respect, I want to mention two other 
models proposed in [66] and in [67] 168] where just the magnetic energy as- 
sociated with such strong magnetic fields is "directly" transformed to power 
the fireball. In these models the central engines of GRBs are ultramagne- 
tized rapidly rotating neutron stars called magnetars that, in fact, are able 
to produce dipole magnetic fields up to 10 15 gauss [Ml EH EI] • In particu- 
lar, in [66] the magnetic energy release to the fireball is supposed to happen 
during the formation of the magnetar as a consequence of the collapse of a 
white dwarf, while in [67, 68j it is supposed to happen during the collapse 
of the magnetar into a black hole. 

As a conclusion of this Paragraph and of the previous one, I want to stress 
the fundamental role that the presence of strong magnetic fields around the 
central engine has in the production of a GRB: 

• it is necessary to explain the high linear polarization degree of the 
GRBs gamma-ray spectrum; 

• it is invoked to account for the energy "extraction" from the central 
engine to power a GRB or, even, to account for the energy itself to 
power a GRB. 

I have stressed the importance and the role of such strong magnetic fields 
because, as we will see in the following Chapters, their existence around 
astrophysical compact objects represents one of the fundamental hypotheses 
of this work and its most important experimental counterpart. In fact, I will 
perform all the calculations by assuming to deal with overcritic magnetic 
fields and having in mind the astrophysical scenario I have just described. 



In the black hole-accretion disk system also a large amount of gravitational bind- 
ing energy is converted into electromagnetic energy by means of the Blandford-Znajek 
mechanism. 
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Chapter 2 

Theoretical tools 



In this Chapter I want to resume the theoretical background I needed to 
obtain the main results of the thesis. The subjects can be found in many 
textbooks and they are also, unavoidably, very different from each other. 
Nevertheless, the goal of the Chapter is to give the theoretical tools to 
derive the final results and, also, to fix the notation that will used in the 
rest of the thesis. 

The Chapter is divided into three Sections. In the first one I quote 
the main results about the motion of a charged relativistic particle in the 
presence of a constant and uniform magnetic field: I treat both the case of 
a classical and of a quantum particle. In the second Section I discuss about 
the effects and the transitions induced on a quantum system by an adiabatic 
perturbation. Finally, the third Section is devoted to a brief introduction to 
quantum field theory in curved spacetime with particular attention to the 
general covariant formalism to deal with a spinor field in the presence of a 
classical background gravitational field. 



2.1 Motion of a charged relativistic particle in the 
presence of a constant and uniform magnetic 
field 

The problem of a charged relativistic particle in the presence of a constant 
and uniform magnetic field can be solved exactly both at a classical and at a 
quantum level. Since this physical system is discussed in many textbooks I 
limit myself to a quotation of the results I will use in the following Chapters. 
A rigorous derivation of these results can be found in [72\ [73"] for what it 
concerns the classical case and in (73J El] for what it concerns the quantum 
case. 

As usual, I choose the reference system, whose coordinates are indicated 
as x, y and z, in such a way the magnetic field lies in the positive z direction 
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that is 




(2.1.1) 



with B > (I use the primed notation for later convenience). For definite- 
ness, the charged particle is assumed to be an electron with rest mass m 
and electric charge — e < 0. Finally, natural units with h = 1 and c = 1 are 
used throughout. 



2.1.1 Classical mechanics 

From a classical point of view, it is well known that the electron has a 
constant z component v z of the velocity while it rotates uniformly and anti- 
clockwise in the x—y plane. The resulting trajectory is, in general, an helix 
with the axis parallel to the z axis (see Fig. 12. ip . The axis of the helix 
intersects the x—y plane at a point Q whose nonzero coordinates have to be 
given as initial conditions and are indicated in Fig. 12.11 as xq and yo. In 
this way the square distance R xy of the axis of the helix from the origin is 
given by R xy = Xq + y\ and it is a constant of motion. Since the Lorentz 
force does not do any work, the energy of the electron is also a constant of 
motion. Instead, the conservation of the z component 

h = xp y - yp x (2.1.2) 

of the electron angular momentum depends on the gauge one chooses for 
the vector potential A'(r) corresponding to the magnetic field B' (the scalar 
potential can be assumed to vanish) In fact, if v x and v y are the x and 
y components of the electron velocity, the momenta p x and p y appearing in 
Eq. (I2.1.2p are defined as 

Px = mv x - eA' x (r), (2.1.3a) 
p y = mv y - eA' y (r) (2.1.3b) 

and, in general, l z is not conserved. Moreover, it can easily be shown that 
if one chooses the so-called "symmetric" gauge in which 



A'(r) = -i(rxB') = f | , J . (2.1.-1) 



-y 

X 





then L is also a constant of motion. 



1 This is not surprising because l z is not in the present problem a true physical quantity. 
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Figure 2.1: Classical trajectory of an electron in the presence of a constant 
and uniform magnetic field directed along the positive z axis. 

2.1.2 Quantum mechanics 

The quantum description of the motion of a relativistic electron in the pres- 
ence of the magnetic field B' is very different from the classical one and it 
depends on the complete set of commuting observables one chooses in the 
Hilbert space of the system. Now, the Hamiltonian of the system is, obvi- 
ously, a conserved quantity and, in the symmetric gauge, it is given by the 
Dirac Hamiltonian 

H' = a - [V + eA'(r)] + (5m (2.1.5) 

where 'P = —id is the linear momentum vector operator and a and (3 are 
the 4x4 Dirac matrices defined by the relations 



u = x,y,z, 



(2.1.6a) 
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{a u , a v } = {a u , (3} = u,v = x,y, z and u / v (2.1.6b) 

where / is the 4x4 unit matrix and the braces indicate the anticommutator. 
In particular, I choose the matrices a and (3 in the Dirac representation 

where 1^ and cr^ with 

■?>-(!;)■ -f-G^)- ^-(i-'O 

are the 2x2 unit and Pauli matrices respectively. 

We have seen before that in classical mechanics the z component of the 
orbital angular momentum is a conserved quantity in the symmetric gauge. 
In quantum mechanics this in not true because the electron spin has to be 
taken into account. In fact, the z component of the electron total angular 
momentum defined as 

=£ Z + ?A = xVy - yVx + 1± (2.1.9) 

with 




(2.1.10) 



is, actually, a constant of motion in the symmetric gaugeo In order to 
check that [H',J^ /2) ] = one only needs the well-known commutator rules 
among the operators r and T 3 : 

[u, V v ] = iS u , v u,v = x,y,z, (2.1.11a) 

[u,v] = [r u ,r v }=0 u,v = x,y,z (2.1.11b) 

and the anticommutator rules (12.1.61) . 

By observing from Eqs. ([2X5]) . flgXjl and (|2"X9|) that both W and 
Jz 1 ^ 2 ^ do not depend on the z coordinate, one realizes that the commuta- 
tors [V Z ,H'] and [V z , ji 1/2) ] also vanish. Now, since a relativistic electron 
has also an internal degree of freedom connected with its spin, the three 
operators Tt' , J^^ 2 ^ and V z are not enough to build up a complete set of 
commuting observables and then to describe in a quantum complete way 
the motion of the electron. Instead, another operator has to be added and 
a good candidate is represented by the quantum operator corresponding to 
the quantity R^. y = Xq + yfi that is, as we have seen, a constant of motion 



2 In the future there will be no possibility of confusion between 2x2 and 4x4 matrices 
and I will omit the superscript "(2)" from the 2x2 unit and Pauli matrices. 
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in classical mechanics. In fact, it can be shown that the operators corre- 
sponding to xq and yo (that will be indicated with the same symbols) can 
be written in terms of the fundamental dynamical operators r and as 

xo = I - (2.1.12a) 

+ § (2.1.12b) 

and, by using these equations and the commutators (|2.1.1ip . it is not difficult 
to show that the operator 



_ Vy\ 2 (y Vx\ 
\ 2 ill) ' \l ' eBj 



Rlv=il- LJ k) + (2-1-13) 



commutes both with TL' and with J^^ 2 ^ . Finally, since R 2 xy does not de- 
pend on z then [V z ,Rl y ] = and I can conclude that the operator set 

S' = {TL 1 ,Jz l ^ 2 \v z ,R^.y\ is a complete set of commuting observables in 
the Hilbert space of the system "relativistic electron in the presence of the 
magnetic field B"'. 

By concluding, I want to give here an alternative expression of the op- 
erator yo that will be useful in Chap. 5. In fact, from Eqs. fj2 . 1 . 121) and 
(I2.1.11|) it can easily be shown that 

[*o,yo] = ^ (2.1.14) 

in such a way in the representation in which xq is simply a multiplicative 
operator the operator yo can be written as 

yo = -^d Xo . (2.1.15) 
Electron and positron modes 

Once the complete set of commuting observables S' has been determined, the 
next task is to find the orthonormal basis of the Hilbert space built up by the 
common eigenstates or "modes" of S' . It is well known that in doing that in 
the relativistic domain one has to face the problem of the appearance of the 
negative-energy modes. This problem is solved through the so-called second 
quantization procedure that gives the possibility to interpret the pathological 
negative-energy modes as positive-energy antiparticle (positron) modes. The 
subject is well known and, for this reason, I will deal here directly with 
electron and positron modes without introducing the negative-energy modes, 
the charge-conjugation operator and so on. 

If the symbol "j" embodies all the needed quantum numbers, the electron 
and the positron modes can be indicated as u'j(r) and v'Ar) respectively and 
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they are, from a mathematical point of view, fourdimensional spinors. Now, 
it can be shown that four quantum numbers n d , k, a and n g have to be 
introduced and that they can assume the following values: 

n d ,n g = 0,1, ... , (2.1.16a) 
k = any real number, (2.1.16b) 
a = ±l. (2.1.16c) 

The physical meaning of these quantum numbers can be understood by 
looking at the eigenvalue equations that the spinors u'j(r) and v'j(r) satisfy: 



7~C'u'j = WjUj, (2.1.17a) 

V z u'j = ku'j, (2.1.17b) 

J^u'j = (n d - n g + ~) u'j, (2.1.17c) 

R >'i = (2 - L17d) 



and 



Ti'v'j = -wjv'j, (2.1.18a) 

V z v'j = -kv'j, (2.1.18b) 

Jz 1/2)v 'j = -(n d -n g + |) v' j: (2.1.18c) 

*M = 2 -^ v 'r ( 2 - L18d ) 
In these equations I introduced the energies 



Wj = yjm 2 + k 2 + eB(2n d + l + a), (2.1.19a) 



wj = ^m 2 + k 2 + eB(2n g + l-a) (2.1.19b) 

that are called electron and positron Landau levels respectively!! I point 
out that the electron (positron) Landau levels do not depend on n g (n^)0 
Now, from Eqs. ()2.1.17bj) . (|2.1.17cl) . (|2.1.18b|> and (|2.1.18cl) one realizes 
that k (—k) and a (— a) are the linear momentum and the polarization of 



3 1 mention the fact that even if I will work in the strong magnetic field regime in which 
B/B cr 3> 1, it can be shown that the radiative corrections to the electron and positron 
Landau levels are logarithmic in the ratio B/B cr [ZS] and then that they can be safely 
neglected. 

4 The electron (positron) Landau levels have a further kind of degeneration, in fact, 
if 3+ = {n d ,k,+l,n g } and j- = {n d + l,k, -l,n s } (J+ = {n d ,k,+l,n g + 1} and j_ = 
{rid, k, —1, n g }) then Wj + = Wj_ (tDj = w~-_ ). I will have to take into account this further 
degeneration in Chap. 5. 
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the electron (positron) along the magnetic field or, in other words, the lon- 
gitudinal linear momentum and the longitudinal polarization of the electron 
(positron). Instead, Eqs. (|2.1.17cj) . <j2.1.17dj> . (I2.1.18c|) and (|2.1.18dl) show 
that the quantum numbers n d and n g are connected with the transverse 
motion of the electron and of the positron that is with the motion of the 
electron and of the positron in the plane perpendicular to B'. It is worth 
giving here some details on how the quantum numbers n d and n g appear 
in the electron eigenvalue equations (I2.1.17ap . (|2.1.17cp and (|2.1.17dj) [their 
appearance in the corresponding positron eigenvalue equations (|2.1.18ap . 
(I2.1.18cp and (|2.1.18dp can be understood in a completely analogous way]. 
The essential point is that in order to solve these equations one introduces 
the operators ad, a\ and a g , a\ defined as 



9 2 




x — iy) + 




(Vy + lP x ) 



— (x + iy) + \—(V y -iV x ) 
en 



x + iy) 




\Py iV x 



— (x - iy) - J — (V y + iV x 



(2.1.20a) 
(2.1.20b) 
(2.1.20c) 
(2.1.20d) 



By means of these definitions the transverse position and momentum oper- 
ators x, y, V x and V y can be written as 



2^1 ^B^ 9 + a ' 9 + ad + a d^ 
a\ + a d - 4), 




a d + a d ), 



a d ~ a\] 



(2.1.21a) 
(2.1.21b) 
(2.1.21c) 
(2.1.21d) 



(1/2) 



and Rl„, one 

•by 



By performing these substitutions on the operators TC 1 , J- A 

sees that Tt', J^^ 2 ^ and R%. y themselves depend on a^, a g , a\ and a\ only 

through the quadratic operators Ma = &\a>d and M g = a\a g - Now, starting 
from Eqs. (I2.1.20p and from the commutators (12.1. lip it can be seen that 



[a d ,aJ d ] = 1, 
[a g ,a d ] = [a g ,a d ] = 



(2.1.22a) 
(2.1.22b) 
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and then that a<j, a\ and a g , a\ can be interpreted as two independent 
pairs of ladder operators. The corresponding operators Nd and M g can be 
interpreted consequently as number operators and their eigenvalues are just 
the previously introduced nonnegative integers numbers rid an d n g . 

I want to discuss now the structure of the electron and positron modes 
u'j(r) and v'Ar) and their orthonormalization relations. By solving step by 
step Eqs. (I2.1.17P and (I2.1.18P one finds that u'j(r) and v'j(r) are given by 



(2.1.23a) 
(2.1.23b) 



In these expressions I introduced the operator V' and the two twodimensional 
spinors ^p'Av) and x'j( r )- On the one hand, the operator V' is defined as 




P-^(rxB') 



(2.1.24) 



V = a- [V + eA'(r)} = a 

and, by substituting the expressions (|2.1.21j) . it can written as 

V = iV2eB(a a a- - a d a + ) + a z V z (2.1.25) 

with a± = (a x ± ia y )/2. On the other hand, the twodimensional spinors 
ip'j (r) and x'j ( r ) are given by 



<PjK r ) = fa°n d ,n g [X,y)- 



X' j (r)=f_j' w (x,y) 



where 



f' + l 



exp(-ikz) 

7r z 



f'-x 



(2.1.26a) 
(2.1.26b) 

(2.1.27) 



and where the scalar functions 



eB 1 1 

2^ij y. 



(a d ) h (al) l >exp 



eB(x 2 + y 2 



(4) h 



Vh ] - V V 2 J 



1 feB\ l2+ \ 

I — I (x — ly) 2 exp 



eB(x 2 + y 



2M 



(2.1.28) 
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depend only on the coordinates x and y [in this equation the operators a) d 
and a g are supposed to be expressed as in Eqs. (|2.1.20bh and (|2.1.20dj) ]. 

In Eqs. (|2.1,26p I also introduced the symbol L z which is the length 
of the quantization volume in the z direction. In fact, in the course of 
the calculations it is easier to deal with normalizable electron and positron 
modes. To do this, I assume the whole space to be confined between the 
two planes z = —L z /2 and z = L z /2 and the modes to satisfy the periodic 
boundary conditions 

u' J {x,y,-L z /2)=u' J {x,y,L z /2), (2.1.29a) 
v' j (x,y,-L z /2)=v' j (x,y,L z /2). (2.1.29b) 

In this way, the longitudinal momentum k assumes only the discrete values^ 

2tt£ 

k = + — £ = 0,1... (2.1.30) 

L z 

and only in the final results the continuum limit L z — » oo is performed. 
By using the boundary conditions (j2.1.29j) it can be shown that the twodi- 
mensional spinors ^(r) and x'j{ r ) satisfy the following orthonormalization 
relations 

' dr^(r)^,(r) = J dr X '] (r) X ' jf (r) = 5 jd , (2.1.31) 

where Sjji = <5 n(i ,n' ^k,k'^a,a'^n g ,n' g an d where the integrals on z are intended 
to be performed from —L z /2 to L z /2. Analogously, if we calculate the square 
of the operator V' as given in Eq. (|2,1.25p we have 

V' 2 = 2eB [N d (a-a + + a+a_) + a+a.} + V 2 Z = eB (2N d + I + a z ) + V 2 Z , 

(2.1.32) 

and then, from the expressions (12.1.19j) of the electron and positron Landau 
levels, we obtain 

dr V J(r)— ^—(p' jr {r) = Wj - - '- Sj r , (2.1.33a) 

drx'jir) „ V ' „ V [ X'yir) = '±-^<\,y- (2.1.33b) 

J Wj+m Wji + m J Wj+m 

Finally, from these equations and from Eqs. (|2.1.31l) and ()2.1.23p the or- 
thonormalization relations 

j dr^(r)^,(r) = / dru^(r)^-,(r) = S jtj >, (2.1.34a) 
dr^(r)^-,(r) = (2.1.34b) 
can be immediately found. 



1 For notational simplicity, I do not indicate the dependence of k on 
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Transverse ground states 



By looking at the Landau levels (12.1,190 . one sees that there is a class of 
electron (positron) states characterized by the quantum numbers = 
and a = — 1 (n g = and a = +1) whose energies do not depend on the 
magnetic field strength B. Prom a physical point of view, when the electron 
or the positron is in one of these states the energy associated with the 
interaction of the particle spin with the magnetic field compensates for the 
energy associated to the particle transverse rotational motion. These states 
will play a fundamental role in this work and they will be called Transverse 
Ground States (TGSs) to distinguish them from the other ones that will be 
generically indicated as excited Landau levels. In fact, as I have said at the 
end of the previous Chapter, I will deal mostly with pair production in the 
presence of strong magnetic fields such that B/B cr 3> 1. In this regime a 
TGS characterized by a given value k of the longitudinal linear momentum 
has an energy much smaller than that of the other excited Landau levels with 
the same k. In this respect, it is intuitively understandable that if there is 
no any other dynamical constraint the electrons and the positrons are more 
likely to be produced in TGSs than in other states. For this reason, it is 
worth giving here the explicit expression of the electron and positron TGSs. 
In order to simplify the notation, I will label them only by two indices n 
and k (since there is no possibility of confusion I omit the indices "d" and 
u g v on n d and n g ) that is [see Eqs. (12.1.231) and (12.1.261) ] 



u n,k( r ) — u 0,k,-l,n( r ) 



v 'n,k( r ) = v n,k,+l,o( l 



£k + m 
2e k 



( o \ 

1 



k 



2e k 



\ e k + mj 

( \ 

k 



Ek + m 


V 1 J 



exp(ikz) 
VIZ. 



(2.1.35a) 
exp(— ikz) 



(2.1.35b) 



where [see the second equality in Eq. ([2.1.280 ] 



n+l 



9' n (x,y) = 9' nfi (x,y) = J ^ ) (.r - ///)" cxp 



--r(x +y ) 



and where 



£k = WQk-m — W nt k,+ 



io = \J m 2 + k 2 



(2.1.36) 
(2.1.37) 
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are the energies of the TGSs that have the same expression for the electrons 
and the positrons. Another feature of the TGSs that it is worth stressing 
is that they are all eigenstates of the operator o z with the same eigenvalue 
— 1. I will show in the next Chapter how this feature will give the possibility 
to state some very important selection rules concerning the pair production 
process from vacuum in the presence of a strong, uniform and slowly- varying 
magnetic field. 



Rotated magnetic field 

All the previous results and equations are valid when the magnetic field is 
given by Eq. (|2.1.ip that is when it lies on the positive z direction. Actually, 
I will also deal with a magnetic field B that, more generally, lies on the y-z 
plane that is 

/ \ / \ 

(2.1.38) 



with 









( o 


H 




HI 


sin "d 




KB.) 







B = ^Bl + B 2 z , (2.1.39a) 
tantf=^. (2.1.39b) 

The magnetic fields B and B' have been assumed to have the same strength 
in order to exploit the results obtained in the previous Paragraphs. In fact, 
the magnetic field B can be obtained from B' by rotating the sources of B' 
clockwise around the x axis by an angle i?. For this reason, by introducing 
the rotation unitary operator 

n%M($) = expHtfJ-j 1/2) ) (2.1.40) 

with J^t^ 2 ^ = yV z — zV y +a x /2 the x component of the electron total angular 
momentum operator, some assertions can be immediately stated: 

• a complete set of commuting observables S in the Hilbert space of the 
system "relativistic electron in the presence of the magnetic field B" 
is built up by 

1. the Hamiltonian 

H = nWW{$)H'K£ /2 \&) =a[P + eA(r)] + /3m (2.1.41) 

with 

A(r) = ~(rxB) (2.1.42) 

the vector potential in the symmetric gauge corresponding to the 
magnetic field B, 
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2. the longitudinal (with respect to B) linear momentum 

V\\ =^ 1/2)t (tf)P z ?4 1/2) (tf), (2.1.43) 

3. the longitudinal total angular momentum 

^(1/2) = ^(V2)t WJ a/2)^(l/2) (2L44) 

4. the square transverse (with respect to B) distance 

Rl=n^(mly^W; (2.1.45) 

• the spinors 

Uj (r) = n£W(M(*)> ( 2 - L46a ) 

Uj(r) =^ 1 / 2 )t( 1? )^.( r ) (2.1.46b) 

build up the common eigenstates of the operators in S and are an 
orthonormal basis of the Hilbert space of the system; 

• the spinors uj (r) and Vj (r) satisfy the eigenvalue equations 

TLuj = WjUj, (2.1.47a) 

V\\Uj = kuj, (2.1.47b) 

J^ /2) u 3 = (n d - n g + |) Uj , (2.1.47c) 

,2„.._ 2n fl + l.,. (2.1.47d) 



and 



= -tOj-Uj, (2.1.48a) 

= —kvj, (2.1.48b) 

J^' 2) Vj = -{n i -n g + ^Vj, (2.1.48c) 

= ^j^vj, (2.1.48d) 



and the orthonormalization relations 



/^IVW./^W,)^, (2.1.49a, 
y drut(r)uj-/(r) = 0. (2.1.49b) 
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2.2 Adiabatic perturbation theory 

The adiabatic perturbation theory concerns the effects induced on a physi- 
cal system by an external perturbation whose typical time evolution is much 
larger than the typical free time evolution of the system itself. I am inter- 
ested here in adiabatic perturbations acting on a quantum system whose 
time evolution is described by the Schroedinger equation [THl [77]. I also 
suppose that 

1. the system under study is characterized by a Hamiltonian H(£(t)) that 
depends on time through a parameter £(t) representing the external 
perturbation; 

2. for any fixed value £ the eigenvalue equation 

H(0\n,0 = en(0\n,0 (2-2.1) 

can be solved exactly and all the resulting quantum numbers embodied 
in the symbol n are discrete. 

Now, the evolution of the system when varies with time is determined 
provided the Schroedinger equation 

4| = (2-2-2) 

is solved. Since, at any time t the states \n, build up an orthonormal 
basis of the Hilbert space of the system, I can write the state \t) in the form 



\t) = ^a n (t)exp -i dt' e n (£{t')) 
n L Jt 



\n,m) (2-2.3) 



where it has been assumed that the perturbation starts changing with time 
at to an d where the complex coefficients a n (t) have to be determined. By 
substituting Eq. (|2.2.3f) in Eq. (|2.2.2f) and by projecting the resulting equa- 
tion on the state |m, £(£)), I obtain the following differential equation for 
the coefficient a m (t): 



da n . 
dt 



+ ^2(m,C(t)\d ! :\n,C(t)}i(t)a n (t)exp if dt' Ae mn (£(t')) 

n L Jt 



= 0, 
(2.2.4) 

where 

Ae mn (£(t)) = e m (£(i)) - e n (£(t)). (2.2.5) 

Obviously, Eq. (|2.2.4|) is equivalent to the Schroedinger equation (|2.2.2|) . 
Nevertheless, since the parameter £(t) changes slowly with time, its time 
derivative £(t) is a small quantity (in a sense that will be specified at the 
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dt"Ae mi (£(*")) 



end of this Section) and then a perturbative solution of Eq. (12.2,40 in powers 
of can be built up. Up to zero order in £(t) one simply obtains am if) = 
a m (to). Also, by assuming that the system is at time to in the state labeled 
by the index i and that aj(io) = 1, then a™ (f) = S m> i and the first-order 
coefficients a$ (t) are given by 

rt' 

I to |_ Jto 

(2.2.6) 

This expression can be put in a more useful form. In fact, by deriving Eq. 
(12.2. ip with respect to the parameter £ and by projecting on \m, £), one 
obtains 

<m,e|ftff(0]|n,0 = [% n (0]<5n,m " ^mn(0 <m, £\d 6 \n, 0- (2-2.7) 

By assuming that the eigenstates \n, £} are nondegenerate the previous equa- 
tion is equivalent to the following ones: 

(n,{|S s |n,{) = is { (n,{|n,{) = 0, (2.2.8a) 
K^K9=- h y* fl (2.2.8b) 

I point out that in writing the first equation one tacitly exploits the following 
facts: 

1. being the states |n, £) nondegenerate, they can be chosen to be real; 

2. being the indices n discrete, the states \n, £) are normalizable. 

In conclusion, by substituting the previous equations in Eq. (|2.2.6|) one 
obtains 

a\ 1 \t) = 0, (2.2.9a) 



a m \t) = I dt — t f/j.i\\ eX P 



to 



Ae mi (£(t')) 



f 

dt"Ae mi (at")) 



to 



iim^i (2.2.9b) 



with 



iWO = (m,£(f)|#(e(0)|i>£(0) = <m,£(Olfttf(£(O)]lU(O>£(*0- 

(2.2.10) 

In particular, Eq. (|2.2.9ap implies that the coefficient dj(i) is up to first 
order in £(t) equal to one. In other words, in the previous approximations 
the depletion of the initial state is at least a second-order effect in £(i). 

Finally, it must be pointed out that the first-order perturbative approxi- 
mation can be safely used only if the variation of the coefficients oS' (t) with 
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m 7^ i is much smaller than one in the typical time given by the inverse of 
the Bohr frequency Ae m j (£(£)), that is only if 



Hfni (i) 



[Ae mi (£(t))] 2 
at any t > in and for any m ^ i. 



<€. 1 



(2.2.11) 



2.3 Quantum field theory in curved spacetime 

It is well known that, despite many efforts, a completely satisfying quantum 
theory of the gravitational field is still missing. Consequently, the interaction 
of a matter field and of the gravitational field has not been described in a 
completely quantum way. A less ambitious attempt to take into account the 
effects of the gravitational field on a matter field consists in |78} 1791 [80] 

1. treating classically the gravitational field itself as a modification of the 
metric properties of the spacetime; 

2. quantizing the matter field in the resulting curved spacetime. 

I do not want to go into the conceptual difficulties that quantum field theory 
in curved spacetime involves and the exact physical conditions that have to 
be satisfied in order that it can be safely applied [see in particular [75]]. In 
general, it can be said that if the typical lengths involved in the problem 
under study are much larger than the Planck scale Lp = \[G = 1.6 x 
10~ 33 cm with G the gravitational constant, then the quantum effects of 
the gravitational field can be neglected and it can be treated classically. In 
this dynamical regime the gravitational field is described mathematically by 
means of the metric tensor g l _ LV {y) of the spacetime where x = . . . , x 3 ) 
characterizes a generic point of the spacetime itself. When no gravitational 
field is present the spacetime is the Minkowski spacetime and, by definition, 
it is possible to choose the coordinates in such a way the metric tensor 
has the simple structure r} a p = diag(l, —1, —1, —1)11 The form itself of the 
metric tensor r} a p suggests a natural way to distinguish a time coordinate 
from the other three space coordinates. For this reason I will indicate the 
four coordinates of a generic event in Minkowski spacetime by using the 
noncovariant notation (t, r). 

Now, in this work I will deal essentially with electrons and positrons 
then I will be interested only in spinor matter fields. The treatment of a 
spinor field in general relativity is more complicated than that of a scalar 

6 For future notational convenience, the tensor indices of quantities referring to 
Minkowski spacetime are indicated with the first Greek letters a, /3, 7 or 8, while those of 
quantities referring to a generic curved spacetime are indicated with the Greek letters in 
the middle of the Greek alphabet such as A, /i, v or p. 
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field or, in general, of a tensor field and it requires a detailed analysis which 
is independent of the quantization procedure. For this reason in the next 
Paragraph I will discuss the "pedagogical" case of the quantization of a free, 
real scalar field in curved spacetime and then in Par. 12.3.21 1 will outline the 
formalism to describe a spinor field in a general covariant way. 

2.3.1 Quantization of a free, real scalar field in curved space- 
time 

The study of quantum field theory in curved spacetime gives the possibility 
to understand that the quantization of a matter field in Minkowski space- 
time and the consequent interpretation in terms of identical particles can 
be carried out in a coherent way only because of the particular metric fea- 
tures of this spacetime. Actually, in a completely generic spacetime it is 
impossible to characterize what a quantum particle is and how to detect it 
|78j . In the present Paragraph I firstly review the main steps one follows in 
quantizing a free, real scalar field in Minkowski spacetime and then I show 
which features a curved spacetime has to share in order that the same real 
scalar field can also be quantized coherently in it. In this respect, I observe 
that a scalar field in curved spacetime is intended as a general scalar that is 
a scalar quantity under general coordinate transformations, while a scalar 
field in Minkowski spacetime is intended as a Lorentz scalar that is a scalar 
quantity under Lorentz transformations. Nevertheless, with an abuse of no- 
tation I will indicate both these fields and their related quantities with the 
same symbol. 

If I call <f)(t, r) a free, real and Lorentz scalar field, then its Lagrangean 
density ££ is given by 

if = rf\d^){d^) - ™l<f (2-3.1) 

where, after the quantization of the field itself, ttlq will be interpreted as 
the mass of the resulting identical particles. The Lagrangean density if is 
assumed to be a Lorentz scalar in such a way the action 

S = j dtdv^ 

is also a Lorentz scalar and the equation of motion of the field 

(<9f - d 2 + m 2 ) <f> = 0, (2.3.3) 

obtained from the stationary condition 5S = 0, is covariant. The next step 
in the quantization procedure consists in finding the general solution of Eq. 
(|2.3.3p as a sum of normal modes that build up, in fact, an orthonormal 
basis of the space of the solutions of Eq. (|2.3.3p itself. As it follows from 
the general theory of partial differential equations, the normal modes can be 



(2.3.2) 
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unambiguously determined only after imposing some boundary conditions 
at a fixed time t = to- If one chooses to impose periodic boundary conditions 
in a finite cubic volume V = L 3 , the modes 4>k(t, r) are characterized by the 
three discrete momenta k = ±2£ir/L with £ a vector of natural numbers. In 
particular they can be written as 

^' r) = eXP[ "^ k " r)] ^ 



with iOk = yrr^+W and they are orthonormal in the sense that they satisfy 
the relations 

(<Ak,<Ak') = #k,k' (2.3.5) 

with the scalar product between two scalar functions (/>i(t, r) and (f)z(t,r) 
defined as 



I dr[<h(t,r)dt<&(t,r) - <fc(t,r)dt<h(t,r)l (2.3.6) 
Jv 



By using the modes (f>u.(t, r), the general solution of Eq. (12.3.30 can be 
written as 

<f>(t, r)=J2 [aicMt, r ) + 40k(*, r)] (2.3.7) 

k 

and the quantization of the field consists in transforming the complex num- 
bers ak into operators and in imposing the quantization rules 

[ak,oL] = 4,k', (2.3.8a) 
[a k ,a k /]=0. (2.3.8b) 

Finally, the particle interpretation of the quantized field is achieved by in- 
troducing the total linear momentum P of the field: 

P = - / drirdcf) (2.3.9) 
Jv 

and the total Hamiltonian H of the field: 

H= I dr{-ad t (j) - JSf ) (2.3.10) 
Jv 

where 

is the momentum field conjugated to <p(t,r). In fact, by substituting the 
expansion (|2.3.7j) in Eqs. (|2.3.9p an d (|2.3.10p and by exploiting the or- 
thonormalization relations (I2.3.5P and the quantization rules (12.3. 8p . one 
easily finds that 

P = k («kal + 4 fl k) = Yl kiV k, (2.3.12a) 
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H = - "hiatal + ala k ) = £ uj k N k + W (2.3.12b) 

k k 

with Nk = a^ak and with Wq = 1/2 ^ k Wfc an infinite constant that will 
be discussed below. Now, by using again Eqs. (|2.3.8p one sees that the 
operators commute among them and that they have nonnegative integer 
eigenvalues n^. Their generic common eigenstate {n^) is such that 

P K> = (^2 kn ^j K), (2.3.13a) 

H\n k ) = <^«k + W ^j |n k ) (2.3.13b) 

and it can be interpreted as the state in which n k relativistic particles with 
mass mo, i. e. particles with fourmomentum (wfc,k) = (y 7 m\ + k 2 , k) are 
present. In this picture the infinite energy Wq is interpreted as the energy 
of the vacuum state |0) characterized by = for all k. In Minkowski 
spacetime and, in general, when the time evolution of the spacetime metric 
is assigned, the presence of Wq in Eq. (j2.3. 13bf) does not play any role 
because it is just a constant zero-point energy and it is neglected 

In what follows, I want to do the analogous steps I have already done 
in Minkowski spacetime but in a spacetime with general coordinates x = 
(x°, . . . , x s ) and with metric tensor g^ u (x). Firstly, I generalize the La- 
grangean density (|2.3.1|) to transform it into a general scalar Lagrangean 
density. The easiest way to do this is to replace the Minkowski metric ten- 
sor r] al3 with g^ u {y) and to multiply the resulting Lagrangean density by 
V~9(x) witn 0( x ) = det(g^(x)): 



J? = v^xj [<T (x)(0„0)(&0) - mfo 2 ] ■ (2.3.14) 

Actually, other terms proportional, for example, to the scalar curvature 
of the spacetime could be added to this Lagrangean density but I am not 
interested in them here. Analogously, the action is defined here as 



S = J d 4 xi? (2.3.15) 

in order that it is also a general scalar and the equation of motion 

J_ a M (V-g(x)^) +m 2 <f> = (2.3.16) 
V~9W 



7 The presence of Wo must be taken into account if one wants to compute the backre- 
action effects that the matter field can have on the time evolution of the spacetime metric 
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is generally covariant. Now, analogously to what I have done in Minkowski 
spacetime, I have to determine the normal modes building up an orthonor- 
mal basis of the solutions of the partial differential equation (12.3. 16)) . Prelim- 
inarily, I also have to define a threedimensional hypersurface, corresponding 
to the ordinary threedimensional space in Minkowski spacetime, where 

1. the scalar product between two general scalar functions can be defined; 

2. some boundary conditions can be imposed in such a way a set of normal 
modes can be unambiguously determined. 

Now, a rigorous treatment of these subjects can be found in [81j and here 
only some stated results will be quoted. In particular, it can be shown 
that a Cauchy hypersurface, that is a spacelike hypersurface such that any 
timelike or lightlike curve intersects it only once, shares both the previous 
features (SH [73]. A generic spacetime does not admit the existence of a 
Cauchy hypersurface. Instead, it can be shown that the so-called global 
hyperbolic spacetimes not only admit a Cauchy hypersurface but they can 
also be "foliated" by means of Cauchy hypersurfaces. By ignoring here the 
exact mathematical definition of a global hyperbolic spacetime [it can be 
found, for example, in [81]], one can conclude that the topological structure 
of a globally hyperbolic fourdimensional spacetime is M x M 3 where E will 
act as the time axis and M 3 as the threedimensional manifold representing 
the ordinary space. If one also assumes that each "leaf" of the foliation can 
be represented by an equation like x° = const., one can consider x° as the 
time coordinate and ( X if X ■ X ^ clS the three space coordinates. In this way, 
to work in a globally hyperbolic spacetime allows to 

1. define, analogously to Eq. (|2.3.6|) . the scalar product between two 
general scalar functions <^>i(x) and <^2( x ) as 



where E is, in fact, a Cauchy hypersurface (it can be shown that the 
value of the scalar product does not depend on the Cauchy hypersur- 
face one chooses) and n^(x) is the normal versor to S at x; 

2. assume that the normal modes of Eq. (|2.3.16p can be determined. 

Now, I remind that also in Minkowski spacetime when a matter field 
has to be quantized in the presence of an external (for example, electric or 
magnetic) time-dependent field, then the definition of the "energy" of the 
particles and then of the particles themselves is obscure (3j 2] . An analogous 
problem arises in quantum field theory in curved spacetimes when the metric 
tensor g^(x) depends on the time coordinate x°. It can be shown that, in 
order to avoid these further difficulties, one has to assume the spacetimes to 




(2.3.17) 
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admit a Killing vector which is everywhere timelike |82} I78j. In particular, 
this feature implies that the spacetime is static meaning that it is possible 
to choose the coordinates x in such a way that goi( x ) = and dog fMU (x) = 0. 
At this point the quantization procedure is completely analogous to that in 
Minkowski spacetime and I will not repeat the remaining steps. I only point 
out that the definition of the Hamiltonian density of a quantum field is a 
controversial operation in curved spacetime [H31 [HH ES] ■ For definiteness, I 
shall adopt the same definition given in [86] that is, actually, the same one 
gives in Minkowski spacetime. In particular, in the present case one firstly 
introduces the momentum canonically conjugated to the field 0(x) as 



2.3.2 Spinor fields in general relativity: the tetrad formalism 

We know that by using the equivalence principle and the principle of general 
covariance, a recipe can be given in order to transform a Lorentz covariant 
equation into a general covariant equation |82j : 

1. replace the metric tensor r] a p with g^ u (x); 

2. replace the derivatives with the covariant derivatives. 

Actually, in doing so one tacitly assumes to interpret all the Lorentz tensor 
fields as general tensor fields and this can be done only because a tensor 
under general coordinate transformations is also a tensor under Lorentz 
transformations. Now, this procedure can not be carried out when equations 
contain spinor quantities because there is not a representation of the group of 
the general coordinate transformations that behaves like a finite dimensional 
spinor representation under the subgroup of the Lorentz transformations [82\. 
For this reason the previous recipe to transform a Lorentz covariant equation 
into a general covariant equation becomes inapplicable if the equation itself 
contains spinor quantities and another procedure has to be followed. This 
procedure is called vierbein or tetrad formalism and it consists, firstly, in 
introducing normal coordinates at any spacetime point X. By means of 
these coordinates the metric tensor g^ u (x) can be written in terms of another 
set of coordinates x = (x°, . . . , x 3 ) as 




(2.3.18) 



and then defines the Hamiltonian density as 



(2.3.19) 




(2.3.20) 



where the quantities 




x=X 



(2.3.21) 
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build up the so-called tetrad field. The tetrad field has the peculiar property 
to transform as a fourvector both under a general coordinate transformation: 



^-^(x) =► e «(x) ^ ^e«(x) (2.3.22) 
and also under a local Lorentz transformation A^(x): 

yS^A$(x)y? =*► e«(x) ^ A^(x)e£(x). (2.3.23) 

For this reason, by contracting e^(x) with a contravariant general fourvector 
F^(x) as 

V a (x) = e£(x)V^(x) (2.3.24) 

one obtains a set of four general scalars. Analogously, by using the metric 
tensor r) a p (r/ a/3 ) to lower (rise) the indices a, (3 and so on, and the metric 
tensor ^(x) (g^ u (x)) to lower (rise) the indices /U, z/ and so on, that is by 
putting 

e a/1 (x) = r/ Q/3 e^(x), (2.3.25a) 
e°"»(x) = ^(x)e«(x), (2.3.25b) 
eS(x) = »/^(x)^(x), (2.3.25c) 

one can also transform, for example, general covariant tensors T^ u (x) into a 
set of general scalars T a p(x) by simply defining 

r^(x)^e^(x)e^(x)T^(x). (2.3.26) 

Now, let consider only the particular problem of transforming a man- 
ifestly Lorentz scalar Lagrangean density into a manifestly general scalar 
Lagrangean density. At first sight, one can simply assume that all the fields 
appearing in the initial Lagrangean density are sets of general scalars fields 
like V a (x) or T Q/ g (x) , but this is not enough because a problem arises when 
derivative terms are present. In fact, let suppose that 

1. V&(x) is a generic multicomponent field that 

la. is a set of general scalar fields [like the fields V a (x) or T Qj g(x)]; 

lb. transforms as ^(x) — ► D(A(x))\I/(x) under the local Lorentz trans- 
formation A2(x) with D(A(x)) a certain matrix representation of 
the Lorentz group; 



dx v 



2. Vl/(x) appears in the Lagrangean density to be transformed through a 
term like d a ^(x) = e%(x)d^(x). 
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Since ^^(x) transform as a fourvector under general coordinate transforma- 
tions then d a ^f(x) is a set of general scalars. But, the equivalence principle 
requires that special relativity should apply in locally inertial reference sys- 
tems regardless the local inertial reference system one chooses. In other 
words, the final Lagrangean density must also be invariant under a local 
Lorentz transformation. Now, the quantity d a ^>(x) transforms under the 
local Lorentz transformation A?(x) as 

0«¥(x) - A£(x)£>(A(x))fy¥(x) + Ag(x)[fy£>(A(x))]tt(x) (2.3.27) 

and the presence of the second term would destroy the invariance of the 
Lagrangean density under the same transformation. For this reason, one 
introduces the "covariant" derivative 

@ a = eg(x) [fy + r M (x)] = e»(x)% (2.3.28) 

where the quantities r^(x) are matrices of the same dimension of the repre- 
sentation JD(A(x)) and are called connections. By assumption, the connec- 
tions r^(x) transform under the local Lorentz transformation A^(x) as 

r M (x) -> D(A(x))r M (x) J D- 1 (A(x)) - [^Z)(A(x))] J D- 1 (A(x)) (2.3.29) 

in such a way that [see Eq. (|2.3.27|) ] 

9^{x) -> Ag(x)D(A(x))%tt(x). (2.3.30) 

As a conclusion, fry assuming the fields appearing in the initial Lorentz scalar 
Lagrangean density as sets of general scalars and by substituting the ordinary 
derivatives d a with the covariant derivatives @ a , the resulting Lagrangean 
density is a manifestly general scalar quantity. Actually, one requires that a 
Lagrangean density is a general scalar density quantity that integrated on dx 
gives a general scalar action, then, after all, the resulting Lagrangean density 
must be multiplied by y 7 — g(x) with g(x) = det(g fll/ (x)) [57]. Obviously, by 
excluding this last point, the previous recipe can be followed to transform 
any Lorentz covariant equation into a general covariant one. In particular, 
it can be shown that when one deals only with tensor fields it coincides with 
the more usual recipe given at the beginning of the Paragraph. 

Now, if the generic field ^(x) is, actually, a fourcomponent spinor field 
then it can be shown that the connections T^(x) are given by [82] 

IV(x) = -ia^<(x)[^ + r^(x)]e M (x) (2.3.31) 

where 

rjUx) = [d^g vp (x) + d v9lv (x) - d p9llu (x)} (2.3.32) 



y'x a 



A£(x)y: 



Chapter 2 



35 



are the Christoffel symbols [57] and where a a ^ = i[^ a , 7*/2 with 7 = 
0, 7 1 = f3a x , 7 2 = (3a y and 7 3 = /3a z are the covariant Dirac matrices 
satisfying the anticommutation relations {7°, 7^} = 2rf l P . In this way, by 
reminding that the Dirac Lagrangean density in Minkowski spacetime can 
be written as 

if = 1 [% Q (d a ^) - (9 a ^)7 Q *] - mf ^ (2.3.33) 

and by following the rules just stated, the general scalar Dirac Lagrangean 
density is given by 

^= y Z 5M{^[^7 Q ^(x)(^)-(*i M )e^(x)7 a ^] -m**| 



(2.3.34) 



It is worth doing some observations about Eq. (12.3.340 : 

1. as in quantum field theory in Minkowski spacetime, the adjoint field 
'f(x) and the Hermitian conjugated field ^(x) are connected here by 
the relation *&(x) = ^(x)j°; 

2. the "left" covariant derivative is defined as 

i^^-iyx) (2.3.35) 

where the symbol "3^" means that the partial derivative acts on the 
spinor on its left and the minus sign comes from the fact that the 
Lagrangean density must be real and that rj^(x) = — 7°r^(x)7°; 

3. the matrices 7**(x) = e«(x)7 a are called general covariant Dirac ma- 
trices and they satisfy the anticommutation rules {7^(x), 7 iy (x)} = 
2<T(x). 

By applying the usual variational method to the action S = J dxif , one 
obtains the general covariant Dirac equation in the form 

iy*(x)^* - m* = 0. (2.3.36) 

Finally, the scalar product between two spinors ^i(x) and ^(x) is defined 
in a general covariant way as 

(*!, * 2 ) = J dSn^(x)^ 1 (x) 7M (x)^ 2 (x) (2.3.37) 

with S a Cauchy hypersurface, n^(x) the normal versor to £ at x and 7^(x) = 
5'^(x)7 l/ (x)- 
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Pair production in a strong, 
uniform and slowly- varying 
magnetic field 

In this Chapter I present my original results about the electron-positron pair 
production in the presence of a strong, uniform and slowly-varying magnetic 
field. The Chapter is divided in two parts. In the first one the theoretical 
assumptions and the general features of my approach are given (Sect. 13.11) 
while in the second one (Sect. I3.2p the presence probability of a single pair 
is calculated by considering two different magnetic field time variations. 

3.1 Theoretical model 

As I have said in the Introduction and in Chap. 1, my work mainly concerns 
the production of pairs in the presence of strong (3> B cr = 4.4 x 10 13 gauss), 
uniform and slowly- varying magnetic fields. I will show in the next Para- 
graph that these assumptions about the structure of the magnetic field are 
justified in the astrophysical scenario I imagine to apply to my calculations 
and that I have described in Sect. 11.21 Since the magnetic field is assumed 
to be slowly- varying, I will perform in Par. 13.1.21 some preliminary calcula- 
tions in the framework of the first-order adiabatic perturbation theory that 
I will use in Sect. 13.21 

3.1.1 General assumptions 

As I have said in Par. 11.2.31 I imagine the magnetic field I deal with to 
be produced by a (forming or collapsing) magnetar or by a forming black 
hole. In particular, in this Chapter I want restrict my considerations to the 
production of pairs around magnetars in order to simplify the treatment 
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by neglecting all the effects due to the compact object gravitational fieldQ 
Now, we have seen that indirect estimates suggest that the magnetic field 
of a magnetar can be safely assumed to be one or two orders of magnitude 
larger than B cr . In general, the magnetic field produced by a magnetar is a 
dipole field and then it is not uniform in space. Also, since the production of 
the fireball happens during the formation or the collapse of the magnetar it 
is reasonable to assume the magnetic field produced by the magnetar itself 
to be time-varying. Now, the physical problem I want to solve clearly has 
two different scales of length and time: a macroscopic one related to the 
typical linear dimension of a magnetar (~ 10 6 cm) and to its typical time 
formation or collapse (~ 1 s) and a microscopic one related to the elementary 
particles to be produced which is given by the Compton length of the electron 
A ~ 10~ n cm. Clearly, the order of magnitude of the macroscopic scale is 
much larger than that of the microscopic one and this allows me to consider 
during the pair production process the magnetic field of the magnetar as 
uniform in space and slowly-varying in time. Below, I will consider two 
different kinds of magnetic field time evolutions and they are both particular 
cases of the following oneH 

B(t) = f B y (t) ) = B(t) f sintf(i) ) (3.1.1) 
\B z (t)J \costf(i)/ 

where 



B(t) = JBm)+Bl(t), (3.1.2a) 



y 

tantf(t) = g|j. (3.1.2b) 

From a macroscopic point of view the time variation of B(t) implies the 
presence of a nonuniform electric field E(i, r). In the symmetric gauge in 
which 

A(t,r) = ~[rxB(t)], (3.1.3) 
the electric field is given by 

E(i,r) = -d t A(t,r) = ~[r x B(t)]. (3.1.4) 



1 The gravitational force exerted on an electron (positron) by a typical 2 solar masses 
magnetar is several orders of magnitudes weaker than that exerted by its magnetic field. 
Nevertheless, the gravitational effects will be taken into account in Chap. 5 where the 
pair production around massive black holes is considered. 

2 Concerning the magnetic field and its components, 1 use, apart from the time depen- 
dence, the same notation 1 have used in the previous Chapter just to have the possibility 
to exploit many results obtained there. 
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Apart from the particular case in which the magnetic field depends linearly 
on time, the electric field also depends on time. Now, since I am interested in 
the production of electrons and positrons from vacuum I assume the calcu- 
lations to be applied to regions where conduction currents are not present]! 
In this case, the Ampere-Maxwell equation d x B — <9(E = would require 
the existence of a correction to B(i) proportional to B(i) and nonuniform 
in space. But, reminding that B(t) is assumed to be slowly- varying, I am 
allowed to carry out all the calculations up to first order in B(t), so, consis- 
tently, to neglect any contribution proportional to B(i). On the other hand, 
as we will see, the presence of the electric field E(t, r) plays an important 
role in the pair production process. Actually, it can be said that just the 
presence of the electric field E(i, r) accounts for the pair production process 
because it supplies the energy to create the pair. Nevertheless, in this re- 
spect, a conceptual difference between the present model and the Schwinger 
pair production mechanism described in Sect. Il.ll must be pointed out: here 
the electric field (I3.1.4P is rotational, then it does not admit a scalar poten- 
tial and the interpretation of the pair production process as a tunnel effect 
is not straightforward. 



3.1.2 Application of the first-order adiabatic perturbation 
theory to the pair production process: preliminary cal- 
culations 

I have just said that the interpretation of the pair production process in 
terms of the Schwinger mechanism is not straightforward in the problem 
under consideration. Also, since the astrophysical models generally refer 
to the strong magnetic fields present around magnetars and not to their 
induced electric fields, I want to put into evidence in the final results the 
role of the magnetic field B(i). For this reason and being the magnetic field 
B(t) a slowly- varying quantity, I will interpret the pair production process 
as a transition induced by the external adiabatic perturbation B(t) on the 
quantum system represented by a second-quantized Dirac field ip{t,r). The 
Hamiltonian of this system is given by 

H(t) = J dr^(t,r)H(t)i>(t,r) (3.1.5) 

with 

H{t) = ac[-id + eA(t,r)]+/3m (3.1.6) 

and it depends on time through the slowly- varying magnetic field B(i) hid- 
den in A(t, r). In order to apply the first-order adiabatic perturbation theory 
to the problem at hand, I have to determine the instantaneous eigenvalues 



3 The electric field (|3.1.4p is divergenceless then the electric charge density vanishes 
everywhere. 



40 



Section 3.1 



and eigenstates of the Hamiltonian (13,1.50 (see Sect. 12. 2p . As a first step, 
I address the problem of quantizing the Dirac field tp(t,r) in the presence 
of the magnetic B introduced in the previous Chapter [see Eq. (|2. 1.381) ] 
which is identical to B(t) just apart from the time dependence. Since I have 
already determined the one-particle orthonormal basis {uj(r), Vj(r)} with 
j = {n-d, k, a, rig] [see Eqs. (12.1.460 and fl2.1.23[> ] . I can expand the Dirac 
field tp(t,r) as 

tp(t, r) = y^[cj exp(— iwjt)uj(r) + d* exp(iwjt)vj(r)] (3.1.7) 

j 

and impose the usual anticommutation rules 

{cj,cj,} = {dj,d^,} = Sjj>, (3.1.8a) 
{cj,cj<} = {cj.d^,} = {cj,dj>} = {dj,dji} = (3.1.8b) 

among the coefficients Cj and dj that are now operators. By exploiting the 
orthonormalization rules (j2.1.49p among the spinors Uj(r) and Vj(r) and 
the previous anticommutators it is easy to see that the second quantized 
Hamiltonian H = J dr ^(^^(f.r) with H given by Eq. (I2.1.4ip [that 
is by Eq. (|3.1.6p without the time dependence] becomes 

H = ^(wjNj + WjNj) (3.1.9) 

3 

where Nj = c^Cj and Nj = d}-dj and where the zero-point energy has been 
set to zero. Obviously, if |0) is the vacuum state defined by the relations 
Cj|0) = dj\0) = 0, the eigenstates of the Hamiltonian (I3.1.9P are the Fock 
states 

IK-}; K» ^ [ci) nh (ip • • • {d]fh (dlfh . . . |o), (3.1.10) 

while the corresponding eigenvalues are given by 

W = ^2(w Jl n Jl +w ] n- n ) (3.1.11) 
I 

with I a, generic integer index. At this point, the instantaneous eigenvalues 
and eigenstates of the time-dependent Hamiltonian (|3,1.5p are simply given 
by the previous expressions (|3.1.10p and (13.1. lip by adding the dependence 
on time. Also, the creation from vacuum of a pair with the electron and the 
positron in the states j and j' respectively at time t is represented by the 
adiabatic transition 

\0(t)}^\jjXt)) = cUt)dUt)\0(t)}. (3.1.12) 
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The probability associated with this transition is given, in the framework 
of the adiabatic perturbation theory up to first order in B(i), by Pjj'(t) = 
\jjj'(t)\ 2 where [see Eq. (I2.2.9bl) ] 

w«> = jf dt '^m%fr v [ l [ dt " K(t " ) + ti *' (t " )1 } (3L13) 

with 

^■'(0 = (jf(t')\H(t')\0(t')). (3.1.14) 

It is worth doing a couple of observations about this formula. Firstly, the 
initial time t = when the magnetic field is supposed to start changing 
represents, in the astrophysical scenario I imagine to apply my calculations, 
the beginning of the magnetar collapse into a black hole or the beginning of 
the magnetar (black hole) formation. Secondly, the usual physical interpre- 
tation of the square modulus of the coefficients a n (t) in Eq. (|2.2.3|) forces 
to interpret Pjj/(t) not as a creation probability at time t but as a presence 
probability at time t that is as the probability that, being the system in the 
vacuum state at the initial time t = 0, a pair is present at time t > with 
the electron in the state j and the positron in the state j' . 

Now, the time derivative of the Hamiltonian H{t) is given by [see Eqs. 
QSXg]) and (|3X6|) ] 

H(t) = B(t) • d B (H(t)) = B(t) ■ I dr^(t,r)d B (H{t))iP(t,r) 

(3.1.15) 



eB(t) 



drij}*(t,T)(r x ct)ip(t,r) 



2 

and it can be written in terms of the electric field (|3.1.4p as 

H(t) = ~ej dr^(t, r)a ■ E(t, r)if>(t, r). (3.1.16) 

This expression points out the fundamental role of the induced electric field 
E(t, r) in the production process and, by reasoning as in first quantization, 
it has a clear physical meaning. In fact, by reminding that the vector a. can 
be interpreted as the one-particle relativistic operator corresponding to the 
velocity of the electron, Eq. (|3.1.16p expresses the fact that the variation 
per unit time of the energy of the electron is equal to the mean value of 
the power supplied to the electron itself by the external electric field E(£, r). 
On the other hand, the time derivative in Eq. (|3.1.15|) can be written in 
terms of the rotated Dirac field ip'(t,r) = lZ^ 2 \ , d(t))')p(t,r), which is the 
field operator when the magnetic field is directed along the z axis and its 
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strength is B(t) [see Eqs. (|2.1.4fi|) ]: 



H(t) 



eB(t) 
2~ 

eB{t) 
2~ 

eB'(t) 



drifj<(t, r)(r x a)tp(t, r) 

dr^(t,r)1l£ /2) (#(t))(r x c*)74 1/2)t (tf (t))V>'(^ 
dr^' t (t,r)(r x a)V>'(t,r) 



(3.1.17) 

where I used the fact that both a and r are vector operators and where 







B'{t) 



cos$(t)B y (t) - sm$(t)B z (t) 
K sm$(t)B y (t) +cos^{t)B z (t) / 




(3.1.18) 



The form of the vector B'(t) shows that the effects of the change of the 
direction and of the change of the strength of the magnetic field have been 
completely disentangled. In this way, one is able to connect the time varia- 
tion of the magnetic field with some selection rules concerning the states in 
which the pair can be created. To do this I define the operator 

T'(t) = -\l dnf^(t,r)(r x a)ip'(t,r). (3.1.19) 

Since T'(t) is a vector operator and the states |0(i)) and \jj'(t)) are instan- 
taneous eigenstates of the total angular momentum operator j9~^ 2 \t) 

J dr^(t,r)J^ 1/2) (t)ip(t,r), then the Wi gner-Eckart theorem allows to state 
the selection rules resumed in Tab. 13.11 where j\\ = — n g + a/2 and 



n d -n ' + a'/2 [see Eq. (|3.1.18|) and the last line of Eq. (13.1.171) ]. 





Transition operator 


Selection rule 


= 
B(t) = 




in +if| =o 



Table 3.1: Selection rules 

The first selection rule in the previous table has an interesting classical 
counterpart. In fact, following the interpretation of the pair creation in the 
Dirac hole theory, this selection rule means that an electron in the presence 
of a slowly-varying magnetic field with constant direction conserves its total 
angular momentum going from a negative energy level to a positive one. 
Analogously, in classical mechanics the angular momentum of the electron 
is, in the same physical situation, an adiabatic invariant [72J. Apart from 
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those in Tab. 13.11 another important selection rule can be obtained. In fact, 
it can easily be shown that 

H jj ,(t')=B'(t)-(jf(t)\T'(t)\0(t)) 

= • J dru'j(t,r)(r x a)v' f (t,r) (3.1.20) 

= B / (t).T} J ,(0 

where I have defined the one-particle transition matrix elements 

T ii'(*) = -\ f dru 'H^ r )( r x a )4(*> r )- ( 3 - L21 ) 

If the magnetic field has a constant direction then = and only the 
matrix elements Tj-, x (t) play a role (see Tab. 13. ljl . Also, it can be seen that 
(r x a)g anticommutes with the one-particle spin operator cr z /2: 

{(r x a) z ,a z } = 0. (3.1.22) 

Then, if two states u'j(t, r) and Vj,(t, r) are eigenstates of a z with eigenvalues 
a and a 1 respectively then the transition matrix element Tjj, z (t) is different 
from zero only if 

a + cr' = 0. (3.1.23) 

Now, even if the magnetic field changes with time, the instantaneous TGSs 
u' n k (t, r) and v' n , k ,(t, r) given in Eqs. (|2. 1.35j) are eigenstates of a z with the 
same eigenvalue —1 and for every t [see the discussion below Eq. (|2.1.37j) ]. 
This implies that if B(t) changes only in strength the creation of a pair in 
which the electron and the positron are both in a TGS is forbidden because 
Eq. (|3.1.23p is not satisfied. In other words, up to first-order adiabatic 
perturbation theory only the rotation of the magnetic field may allow the 
creation of a pair with both the electron and the positron in a TGS. 



3.2 Pair production in the presence of a strong, 
uniform and slowly-rotating magnetic field 

In the first analysis on the production of pairs in the presence of strong, 
uniform and slowly-varying magnetic fields it was considered the case in 
which only the magnetic field strength changes with time [88]. In order to 
avoid the presence of nonuniform corrections to the magnetic field B(t) its 
strength was supposed to change linearly with time [see discussion below 
Eq. (13X41) ] that is 

B(t) = B^(t) =1 J - I J (3.2.1) 
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with Bq, b > and Bq 3> B cr . Since this case was analyzed before my own 
work, I only quote in my notation the final result that is the total probability 
per unit volume that a pair is present at a time t: 



dPf n (t) 1 



dV 16tt\/2 



C(l-5) 1 
256 15?r 



B l r(t) + 1 



? 3/2 



e 3 / 2 (bR ±M f 

(3.2.2) 

where C( x ) 1S the Riemann function with ^(1.5) = 2.61 . . . |89] and where 
the symbol ~ indicates that the result holds asymptotically at large times. 
In Eq. (13. 2, 2D the square of the quantity R±m also appears. Its origin 
and its exact meaning will be explained below when the analogous presence 
probability will be calculated in the presence of a purely rotating magnetic 
field. At the moment, I only point out that 

1. R±m can be interpreted as the typical radius of a circle in the plane 
perpendicular to B 1 t m (t) within which the magnetar magnetic field can 
be assumed as uniform; 

2. the appearance of R±m m Eq. (|3.2.2p can be understood in terms 
of the nonuniform electric field induced by the time- varying magnetic 
field B 1 T in (t). 

Finally, another observation concerns the dependence of the presence prob- 
ability (|3.2.2|) on the square of the time derivative b of the magnetic field 
B l l n (t) that is, in fact, typical when the first-order adiabatic treatment is 
used to calculate a transition probability (the transition amplitude clearly 
is proportional to the time derivative of the perturbation and it has to be 
squared to obtain the corresponding transition probability). 

Now, already in [HH] it was observed that if the magnetic field direction 
does not change with time the production of pairs in which both the electron 
and the positron are in a TGS is forbidden. As I have said in Par. 12.1.21 since 
in the strong magnetic field regime the energies of the TGSs are much smaller 
than those of the other excited Landau levels, it is reasonable to imagine that 
the production of pairs with both the electron and the positron in a TGS 
is strongly favoured (when it is not forbidden by any selection rule). For 
this reason, by remembering the last statement in the previous Paragraph 
[see the discussion below Eq. (|3. 1 .23j) ] . the most natural continuation of 
the analysis started in [88] has been the treatment of the pair production 
process in the presence of a changing-direction magnetic field. In particular, 
the purely rotating magnetic field configuration [90J : 



B r (t) = B; sinftt (3.2.3) 
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with Bf > and 5/ S> B cr allows an easier mathematical treatment of the 
problem. In fact, from Eq. (|3. 1.18j) we have that 

I ° \ 

B' ; (t) = I OB; I (3.2.4) 

and that only the matrix elements 

T 'n'y^ = -\j dru}(t,r)(r x cx) y v' f (t,r) 



-| J dru'j(t,r)(za x -xa z ) v v' j ,{t,r) 



(3.2.5) 



have to be calculated. Now, I have checked that if P gg '(t) = P(\0(t)) — * 
\gg'(t))) is the probability that a pair is present at time t with the electron 
and the positron both in a TGS and Pjj>(t) is the probability that a pair 
is present at time t with the electron and/or the positron in another state, 
then 

< -F « I- (3-2-6) 



For this reason, I am allowed in first approximation to neglect all the tran- 
sitions to pair states in which at least the electron or the positron is not in 
a TGS0 There are only two transition amplitudes different from zero that 
contribute to the creation of a pair with the electron and the positron both 
in a TGS. In fact, since the TGSs are eigenstates of a z the term in Eq. 
()3.2.5p containing a x vanishes. Also, reminding that the electron (positron) 
TGSs are those with = and a = — 1 (n g = and a = +1), the only 
allowed final pair states are the states [see also the expression (|2.1.21ap of 
the coordinate x in terms of a^, a d , a g , a\\ 

\n,k;n±l,-k) = |0, k, -1, n; n ± 1, -k, +1, 0) (3.2.7) 

where, for simplicity, I have omitted the time dependence. The transition 
amplitudes corresponding to these states are given by 



eBf{n + 1) mfl 
32 d 



7n,k;n+i,-k(t) = \/ — exp(ie k t) sme k t, (3.2.8a) 



, eBrn m£l 

7n,fe;n-l,-Jfc W = \/ — iL r exp(ze fc t) sin e k t (3.2.8b) 

61 £ k 



where [see Eq. (|2.1.37|) ] e k = Vm 2 + k 2 . Actually, for later convenience, 
I want to consider the probability that the electron and the positron are 



4 In any case, in Appendix A I study the general features of the amplitudes of these 
remaining transitions. 
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present in wave packets that are linear superpositions of TGSs with the 
same energy. In particular, I will consider the following pair state 

\a e a p ) = \a e )\a p ) (3.2.9) 

where 

oo oo 

\a e ) = ^J& n (a e )|n,fc) = y 6„(q e )|0, k, —1, n), (3.2.10a) 

n=0 n=0 

oo oo 

I«p) = b n>{a P )\n',-k) = ^ b n >(a p )\n', -k, +1, 0) (3.2.10b) 

n'=0 n'=0 

and where 

z n ( 1 \ 
& n (z)^— exp(^--|z| 2 J (3.2.11) 

with 2 a generic complex number. As it is evident from these equations, I 
have chosen as final states a double plane wave state concerning the longitu- 
dinal motion and a double coherent state concerning the transverse motion 
[73] . In the coherent states with \a e \ 1 (\a p \ 3> 1) the electron (positron) 
has a good spatial localization in the plane orthogonal to the magnetic field 
and this will give me the possibility to understand the role of the induced 
nonuniform electric field 

E / (t,r) = i[rxB / (t)] (3.2.12) 

in the process of pair creation. 

Now, it can easily be shown from Eqs. (|3.2.8p that the transition ampli- 
tude from vacuum to the state \a e a p ) can be written in the form 

lk{a e ,a p ;t) = T||fc(*)T-L («e, a p ) (3.2.13) 
with the longitudinal and the transverse amplitudes given by 



7|[fe(*) = V-^ r exp{i£ k t)sme k t, (3.2.14a) 

oo 

7_l (a e ,a p ) = [K( a e)K+i( a p)Vn + 1 + b;(a e )6*_ 1 (ap)v / rT| . (3.2.14b) 

ra=0 

In this way the longitudinal and the transverse part of the presence ampli- 
tude have been disentangled and only the longitudinal part depends on t. 
Obviously, the corresponding presence probability |7fc(a e , a p ; t)\ 2 can also 
be divided into a longitudinal part and a transverse part. The longitudinal 
part is given by 

d P]l (k,t) = | 71| ( M )| 2 ^dk = |^) 2 ^^^ (3-2.15) 
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where I multiplied by the number of states L z dk/(2ir) [see Eqs. (I2.1.30P ] and 
where — * e(k) for continuous k. The corresponding integrated probability 
per unit time is given by 

dP\\{t) _ eB r L z 2 f°° sm2e{k)t 



dt 64vr 



(m nf f dk^i (3.2.16) 



and I can give an asymptotic estimate of this integral by assuming to be 
interested only in times t such that mt = t/X ^> 1 (in the astrophysical sys- 
tem under consideration this assumption is very realistic). The asymptotic 
estimate can be performed by using the method proposed in [91] and the 
result is 



dPu(t) eB r L z fn\ 2 fir 



TT 



J — sin [2mt + - . (3.2.17) 



dt 64-7T \m I V m t ^ 4 



Now, I will calculate the transverse part of Eq. (|3,2.13p or, equivalently, 
its complex conjugate 

oo 

7*(a e ,a p ) = ^2 [b n (a e )b n+1 (a p )\/n + 1 + b n (a e )b n -i(a p )y/n\ . (3.2.18) 

n=0 

From Eqs. (|3.2.1ip I obtain 



oo 

>y ± {a e ,ap) = exp ^ - I ^ 

' n=0 



n! n! 



1 2 i U, 1 2 



(a e + a p ) exp ( - — — ^ apl + a e a p 



(3.2.19) 



then the transverse presence probability of a pair with the electron in a 
coherent state between \a e ) and \a e + da e ) and the positron in a coherent 
state between \a p ) and |a p + da p ) is 

<iP^(a e ,a p ) = |a e +aJ 2 exp[— (|a e | 2 + |aJ 2 )+2Re(a e a p )] — - — -. (3.2.20) 

IT IT 

It can be shown [73] that the phase of a e (a p ) is connected with the azimuth 
of the mean position of the electron (positron) in the plane orthogonal to 
B^(£) while the modulus \a e \ (\a p \) is connected with its mean distance from 
the origin in the same plane. Then, since I am not interested in the exact 
position of the pair, I put 

a e = \a e \ exp(i(f) e ) = rj e exp(i(f) e ), (3.2.21a) 
a p = \a p \ ex.p(icft p ) = rj p exp(i(p p ) (3.2.21b) 
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and I integrate on the angles (ft e and <f> p . The result is 

dP±(Ve,rj P ) = 4r]eVp(Ve + vl) exp[-(??e + vl)]M 27 leV P )dr]edr]p (3.2.22) 

where lo(x) is the zero-order modified Bessel function [89]. In order to obtain 
a more transparent formula I define the variables r\ and 4> by means of the 
equations 

r] e = r]cos(p, (3.2.23a) 
7] p = r]sm(j) (3.2.23b) 

and I integrate on the phase <ft. Since r/ e and rj p vary from to do, r\ and (j) 
vary from to oo and from to ir/2 respectively, then, after performing the 
integral on (j), I obtain the differential probability 

dP ± (r]) = rf [1 - exp(-2?? 2 )] dr). (3.2.24) 

Now, there is an interesting relation between rj and the mean value in the 
states \a e ) and \a p ) of the operator R\ defined in Eq. (|2.1.45p [see also Eq. 
(|2.1.13p ]. In fact, by using the eigenvalue equations (|2.1.17d[) . (|2.1.18d[) and 
the definitions ()3.2.10p . it is easy to show that 

Rl(a e ) = (a e \R 2 ± \a e ) = 2 ^ D + 1 , (3.2.25a) 

-2 



H + l 

eBt 



Rl(a p ) = (a p \R 2 ± \a p } = ~' ,p _ - . (3.2.25b) 



But, since the one-particle energy of the electron (positron) do not depend on 
the quantum number n g (n^) and then on the quantity r/ e (rip), there are no 
dynamical constraints on the maximum value that rj e (rj p ) can assume. Also, 
since the pair production process is imagined to take place in a macroscopic 
astrophysical environment, I can assume that R±(a e ) 3> A and R±(a p ) ^> A 
in such a way 

Ve - ^fR\{ae) » 1, (3.2.26a) 
r$ ~ ^-R\{a p ) » 1 (3.2.26b) 

where I exploited the fact that in the strong magnetic field limit y/eBj S> 
1/A. Even if the previous inequalities hold, the internal consistency of the 
model requires that the quantities R\(a e ) and R\(a p ) can not be too large. 
In fact, on the one hand the magnetar magnetic field must be approxima- 
tively uniform inside a region with typical linear length R±(a e ) or R±(a p ) 
than R±(a e ) <C 10 6 cm and R±(a p ) <C 10 6 cm. On the other hand, it can 



Chapter 3 



49 



be shown, by applying Eq. (|2.2.11f) to the present case, that the first-order 
adiabatic treatment can be safely used only if the conditions 

^l(«e)|^ < 1, (3.2.27a) 

JD cr 

nR±(a p )—^<^l (3.2.27b) 

are satisfied. I observe that these conditions do not constraint very much the 
allowed values of R±(a e ) and R±(a p ). In fact, by assuming Bf = 10 15 gauss 
and n = 1 s" 1 [see the discussion before Eq. ([3.1. ip ], then the conditions 
([3.2.271) are satisfied if R±(a e ) <C 10 9 cm and, analogously, if R±(a e ) <C 
10 9 cm while we already know that the strong inequalities R±(a e ) <C 10 6 cm 
and R±(a p ) -C 10 6 cm must hold. Also, it can easily be checked that these 
previous strong inequalities do not contradict Eqs. (|3.2.26j) . 
Now, if I define the "mean" quantity 



then [see Eqs. (13.2.261) ] 



eB r 



»7 = V ^ + ~ V 2 [R 2 ± («e) + R\ (op)] = \/eBjR± m (3.2.29) 
and I can write the transverse probability per unit area as [see Eq. ([3.2.24p ] 

dP±(R± m ) (eB r )\ 2 



dA± 2ir 



-R 2 ±m (3.2.30) 



where dA± = i^dR 2 ^_ m is the differential transverse area. Since Eq. ([3.2.30P 
does not depend on t, by putting together Eqs. ([3.2.171) and ([3.2.300 and 
by dividing by L z , I obtain the presence probability per unit time and unit 
volume dV = L z dA±, as 

iPrit) _ m >( B, y ^2n* + «/4.) 



dVdt 2vr \AB cr J K ^ m ' 

where the symbol ~ reminds that this is an asymptotic formula valid for 
mt 3> 1. As I have said below Eq. ()3.2.2[) . the fact that the probability 
(|3.2.3ip grows quadratically with R± m can be understood in terms of the 
induced electric field E/(t, r) [see Eq. ()3.2. 12H ] . In fact, for a purely rotating 
magnetic field 

&Bj „ 

Ej y (t, r) + E%(t, r) = ~^x 2 . (3.2.32) 
In order to connect this quantity with R± m I observe that 



x = x± = R±(t) cos cj>(t) 



(3.2.33) 
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where R±(t) and cp(t) are the polar coordinates in the plane orthogonal to 
B^(t). If I define the function Ej yz (R±) as the average of Ej (t, r)+Ej z (t, r) 
on the angle <f>(t), I obtain from Eqs. (13.2. 32p and I13.2.33|) 

Ej yz (R ± ) = (E% + Ej z )^ = ^-^Ri. (3.2.34) 

With this formula the probability (|3.2.31l) can be written in the different 
form 

dP r (t) e 2 B r Ej yz (R ±m ) sin (2mt + Tr/4) 
dVdt ~ 16vr B cr v / 7rmi 

that points out which is the role of the induced electric field in the pair 
production process and why the final presence probability per unit volume 
and unit time depends on R\ m - The fact that the probability itself depends 
only on Ef y (t, r) and Ef z (t, r) can be explained by rewriting the transition 
matrix elements (13. 1.20j) as [see Eq. (I3.1.16P ] 

H jf (t) = -e J drE' ; (t,r) ■ u'j(t, r)av' f (t, r) (3.2.36) 

where 

Ht ; (t,r) = ±[TxB' r (t)] (3.2.37) 

is the induced electric field seen from the frame which rotates around the 
x axis and whose z axis is instantaneously parallel to B^(i). Now, the 
particular structure of the TGSs [see Eqs. (|2.1.35j) ] makes so that only the 
a z term contributes to the transition and this term contains only E+ Z (t,r) 
which is a linear superposition only of Ef y (t,r) and Ef z (t,r): 

E' ;z (t,r) = costlt E rz (t,r) +smQt E ;y (t,r). (3.2.38) 

Obviously, it would be interesting to compare the final presence prob- 
ability (|3.2.3ip obtained here and that obtained in [88] [see Eq. (|3.2.2|) ]. 
Nevertheless, this comparison is made hard because 

1. the time evolutions of the magnetic fields (|3.2.1|) and (I3.2.3P are very 
different from each other; 

2. Eq. (I3.2.2P is a presence probability per unit volume while Eq. (I3.2.3ip 
represents a presence probability per unit volume and unit time; 

3. the physical meaning of the quantities R±m an d R± m appearing in 
Eq. (|3.2.2p and in Eq. ((3.2.31I) respectively is slightly different. 



For these reasons, in order to show explicitly that the pair production process 
is much more efficient in the presence of a changing-direction magnetic field, 
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I considered an analogous problem to that just treated but in the presence 
of the magnetic field [92j : 



Bf(t) 




(3.2.39) 



with b and Bq the same parameters as in Eq. (|3.2.1|) . This magnetic field 
depends linearly on time as B| m (t) but it also changes direction with time. 



Also, as in [88], I have calculated only the presence probability of a pair 
with the electron and the positron in pure Landau levels and not in coher- 
ent states as before. Obviously, unlike the case treated in [88], the dominant 
contribution to the presence probability is given here by the presence prob- 
ability that a pair is created with both the electron and the positron in a 
pure TGS. The calculations are similar to those I have just presented but it 
is instructive to quote some steps. As in the presence of the purely rotating 
magnetic field, if I take into account only the production of a pair with the 
electron and the positron in a TGS, the pair itself can be created here only 
in the states (|3.2.7|) . In this case the presence amplitudes corresponding to 
Eqs. (|3.2.8f) are given, in the asymptotic limit mt S> 1, by 



lin / e(n + 1) mb 



Jin 



en mb 



(3.2.40a) 
(3.2.40b) 



By using the previous amplitudes, I obtain the total presence probability in 
the form 



2 h 2 



em b 
128^0 



n=0 



dk 



2vr ./-oo im- 



+ k 2 )-' 



(3.2.41) 



where I have summed on the transverse quantum number n and I have 
integrated on the longitudinal momentum k. Now, the sum on the quan- 
tum number n can not be coherently extended to infinity. In fact, it must 
be stopped to a given value N l } n (t) corresponding, through the eigenvalue 
equations (|2.1.17d|) and (j2.1.18dj) . to the fixed value R 2 ±M [see also 



P 2 



2Nj n {t) + 1 
eBj n (t) 



(3.2.42) 



with B l f n (t) = \J1Jq~+ (bt) 2 . This value of R\ M must satisfy analogous 
upper limit conditions to those satisfied by R 2 l (a e ) and R 2 l (a e ), that is the 
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stellar magnetic field must be approximatively uniform inside a region with 
typical linear length R±m and 

VRxm^- « I- (3.2.43) 

JD cr 

Also, since the dominant contribution to Pl in (i) comes from the terms with 
large n, the sum on n can be calculated approximatively by means of the 
"semiclassical" substitutions 

n — ► — ~^ R ^ (3.2.44a) 



£ 

n=0 



eB l } n {t) [R 



f ii/ i / -J.M _ 
' I JDi 



(3.2.44b) 



Finally, by performing the remaining integral on fc in Eq. (I3.2.4ip and by 
dividing by the total volume V = L z 7rR^_ M , I obtain the final presence 
probability per unit volume as 



dPWt) 3 



dV vr 



Bfit) 

B rr 



bR±M S 
6AB 



meB . (3.2.45) 



Clearly, the procedure to obtain this quantity is less rigorous than that 
used to obtain Eq. (|3.2.3ip because the pure TGSs are not well-localized 
states while the definition of a probability per unit volume is, in general, 
local. Nevertheless, it can be shown that by using here the more rigorous 
procedure used in the case of a purely rotating magnetic field, the final result 
would be again Eq. (|3.2.45|) but with R± m instead of R±m- 

At this point, to show that the production of pairs is much more efficient 
in the presence of a changing-direction magnetic field, I divide Eq. (|3.2.2j) 
by Eq. (I3.2.45P and, apart from numerical factors of the order of one, I 
obtain 

dPf n {t) B cr B l ; n (t) [B^ ^fB cr \ 3/2 

which shows, in fact, that in the strong magnetic field regime the pair pres- 
ence probability is much larger in the presence of a changing-direction mag- 
netic field. 



3.3 Summary and conclusions 

In this Chapter I have continued the study started in [HS] about the pro- 
duction of electron-positron pairs in the presence of a strong, uniform and 
slowly- varying magnetic field. 
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We had already seen in the previous Chapter that, in the presence of a 
strong magnetic field, the so-called electron and positron TGSs are "priv- 
ileged" states because their energies are independent of the magnetic field 
and then are much less than the other excited Landau levels. Now, by using 
the first-order adiabatic perturbation theory, I have shown in Par. I3.1.2l that 
the time- variation of the magnetic field can be connected with some selection 
rules concerning the state in which the pair can be created [see Tab. 13. II and 
Eq. (|3.1.23p ]. In particular, I have obtained that only if the direction of the 
magnetic field changes with time it is possible to create a pair with the elec- 
tron and the positron both in a TGS. This fact gives the changing-direction 
magnetic field configurations a particular relevance. In other words, the 
possibility of producing pairs with the electron and the positron both in a 
TGS makes a changing-direction magnetic field much more efficient than a 
fixed- direction magnetic field from the point of view of pair production. This 
fact has been verified quantitatively by comparing the two pair presence 
probabilities per unit volume (|3.2.2I) and (|3.2.45l) . Both these probabilities 
have been calculated by means of the first-order adiabatic perturbation the- 
ory and in the presence of a magnetic field with a linear dependence on 
time. But, the probability (I3,2.2p refers to the fixed direction magnetic field 
(|3.2.1|) while the probability (|3.2.45|) refers to changing-direction magnetic 
field (|3.2.39p and, actually, the ratio (j3.2.46H shows explicitly how the pair 
presence probability in the second case is much larger than the correspond- 
ing quantity calculated in the first case. 

Finally, the pair presence probability has also been calculated in the 
presence of a purely rotating magnetic field. Actually, I have calculated 
the presence probability per unit volume and unit time of a pair with the 
electron and the positron both in a coherent TGS [see Eqs. (I3.2.9P and 
(I3.2,10p ]. These states are well spatially localized in the plane orthogonal to 
the magnetic field and this fact gave me the possibility to better understand 
the role of the nonuniform electric field induced by the rotating magnetic 
field. In particular, I have shown that the pair presence probability per unit 
volume and unit time depends on the square of the electric field [see Eq. 
(|3.2.35p ]. Actually, from a physical point of view it was clear a priori that 
the production of a pair would have vanished in the absence of an electric 
field. Nevertheless, in this respect, Eq. (|3.2.35p also shows an additional de- 
pendence on Bj- /B cr that can be interpreted as a direct effect of the presence 
of the magnetic field. 
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Photon production in a 
strong, uniform and 
slowly-rotating magnetic 
field 

Even if the results of the previous Chapter about the production of electron- 
positron pairs are self consistent, I want to use them now to study some 
processes through that the electrons and the positrons created may produce 
photons. In fact, as I have said in Sect. 11.21 a possible application of my 
theoretical calculations is the study of GRBs and of their properties. In 
particular, in Sect. 14.11 1 will calculate the energy spectrum per unit time of 
the photons produced through the annihilation of already created electrons 
and positrons [93J. Analogously, in Sect. 14.21 1 will evaluate the energy spec- 
trum of the photons emitted as synchrotron radiation by already created 
electrons and positrons [94] . Since one of the conclusions of the previous 
Chapter was that a rotating magnetic field primes very efficient mechanisms 
of electron-positron pair production, I will consider in this Chapter only the 
production of photons in the presence of this kind of magnetic field. Fi- 
nally, I have also calculated in [93] the energy spectrum of photons emitted 
directly from vacuum in the presence of the same magnetic field configura- 
tion. Nevertheless, this photon production mechanism is very inefficient and 
I will not report here all the calculations but only some intermediate steps 
and the final photon spectrum (Sect. 14. 3p . 

4.1 Photon production through pairs annihilation 

Before starting the calculations I want to make a general observation con- 
cerning the approximations made about the spatiotemporal structure of the 
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magnetic fields I deal with. As it will be clear, this observation is inde- 
pendent of the photon production process at hand. In fact, as I have just 
mentioned, I will consider here the production of photons in the presence of 
the already introduced strong, uniform and slowly-rotating magnetic field 
given in Eq. (|3.2.3f) and that, for the sake of clarity, I write again: 

( ° ^ 

B/(t) = Bf shim . (4.1.1) 
\ cos VLt J 



Now, in the previous Chapter I was allowed to assume the magnetar mag- 
netic field as uniform (in the production region) and slowly-varying (during 
the production process) because the nonvanishing mass m of the electron 
(positron) provided a microscopic length and time scale given by the Comp- 
ton length A = 1/m. Instead, here I will deal with photons that are massless 
then, in order to keep the previous approximations, / have to assume ex- 
plicitly to restrict my attention to photons with energies u such that £1 <C oj. 
Actually, this strong inequality does not constraint at all the following results 
from a physical point of view. In fact, we have seen that in the astrophysical 
context under consideration $7 ~ Is -1 while the energies of the photons in a 
GRB pulse are typically larger than 10 -2 MeV corresponding to an angular 
frequency of 1.6 x 10 19 s _1 . 

Now, in Sect. 13.21 1 have shown that pairs can be created in the presence 
of a strong, uniform and slowly-rotating magnetic field with the electron and 
the positron both in a coherent TGS. By using the same technique sketched 
at the end of the previous Chapter, one can also obtain the probability per 
unit volume that a pair is present at time t with the electron (positron) in 
any pure TGS with longitudinal momentum between k and k + dk {—k and 
—k — dk). By indicating this probability as f(k, t)dk it can easily be shown 
that 

3 

m 2 (ttR± M ) 2 sin 2 e{k)t (4.1.2) 

with s(k) = V m? + k 2 . It is worth pointing out the dependence of f(k, t) on 
s~ e (k). In fact, this means that in the physical situation under consideration 
the production of high-energy electrons and positrons is strongly suppressed 
and this fact will also affect the production of high energetic photons. 

Now, in Appendix B I show that the quantity f(k, t)dk can be interpreted 
as the mean number of electrons or, symmetrically, positrons per unit volume 
present at time t with a longitudinal momentum between k and k + dk. This 
allows me to use the distribution (I4.1.2P to calculate the annihilation photon 
spectrum per unit time, i. e. the number of photons per unit energy, unit 
volume and unit time that are produced as a consequence of the annihilations 
of the electrons and of the positrons previously created, by using the formula 



eB; 
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mm 



' dfcdfc' — v ' ' J v(k,k')f{k,t)f{k',t). (4.1.3) 



diudVdt J du> 

In this formula, da(k, k' ,uj)/duj is the cross section per unit of photon energy 
of the pair annihilation into two photons process and v(k, k') is the flux factor 
of the colliding particles defined as: 

v(k,k') = J{k,k')V (4.1.4) 

where J(k, k') is the flux density of the colliding particles and V is the quan- 
tization volume. I will give below a more precise definition of the quantity 
da(k,k' ,uj)/duj and its explicit form. At the moment, I observe that this 
picture is already approximated. In fact, I am considering the distribution 
f(k, t) as a given function while its time dependence is affected by the anni- 
hilation process itself. Nevertheless, I neglect this fact for the moment and 
I will deal with it at the end of this Section. 

In order to calculate the quantity da(k, k' , u)/du, I start by writing the 
cross section da of the annihilation of an electron with fourmomentum (e, k) 
and spin s and a positron with fourmomentum (e',k') and spin s 1 into two 
photons with fourmomenta (w,q) and (u/,q') and polarizations A and A' 
respectively. Actually, the process of pair annihilation will take place in the 
magnetic field, which therefore affects all the dynamical processes. However, 
in computing the photon production rate I consider that the main dynamical 
effect of the magnetic field is the pair production. The rest of the process 
will be therefore calculated neglecting the effects of the magnetic field both 
on the annihilation process and on the photon final states. Following this 
approximation, I need the cross section da of the pair annihilation into 
two photons process in the vacuum. This quantity may be found in many 
textbooks and I will quote (in my notation) Eq. (33.2) in 



d(J = p> ^| 2 7^( k + k' - q - q'We + e' - u - J) (4.1.5) 

where uQv is the annihilation matrix clement and ct em = e 2 /(4vr) the fine- 
structure constant. In particular, even if the electrons and the positrons in 
the present case have a preferential direction of the spin [see Eq. (]3. 2. 71) ]. 
I will also neglect this aspect and I will use the cross section summed over 
the photon polarizations and averaged over the electron and positron spins. 
Eq. (|4.1.5p remains unchanged but the quantity |"uQt>| 2 must be replaced 
by [95] 
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with 



x = -%(e - kcosO), (4.1.7a) 



2UJ {e' -k'cosO) (4.1.7b) 



m 2 



where I have assumed that the colliding particles have momenta along the z 
axis and that 9 is the angle between the photon momentum q and this axis 
(this means that k and k' can be positive or negative) 

By using the distribution ([4.1. 2p in Eq. ([4.1, 3D . I take into account, 
from a kinematical point of view, the anisotropy induced by the presence 
of the magnetic field. From the dynamical point of view, I should take into 
account this anisotropy by substituting the matrix element uQv with the 
corresponding one calculated in the presence of the strong magnetic field 

= (0,0, Bf). The calculation is very complicated and it involves the use 
of the Schwinger propagator pj, i. e. the electron propagator in the presence 
of a constant and uniform magnetic field. In Appendix C I will analyze the 
matrix element so calculated in order to give the correction induced by the 
use of the Schwinger propagator and to find the physical conditions in which 
the previous approximated treatment is correct. 

Now, the threedimensional 6 function in Eq. (|4.1.5f) can be exploited to 
perform the integral on q'. Also, since I need the cross section as a function 
of the energy of one of the photons created (the other being fixed by the 
energy 5 function), I integrate the cross section itself with respect to the 
angular variables. The result is 



dcr(k,k',u) wa^ f . . E \uQv\ 2 



duo v(k,k ; ) 



[ rf(cosfl) ^ U< ?^ 8(e + e' - lo - a/W (4.1.8) 

7-i 



where the integral on the azimuthal angle gives a factor tt to take into 
account that the two final photons are indistinguishable and where 



J = |q'| = |k + k' - q| = y/(k + k') 2 +lo 2 - 2uj(k + k') cos 9. (4.1.9) 

By substituting this expression of lo' in the argument of the 5 function in 
Eq. (14.1.8)) I can perform the integral on cos 9 and the final result is 



da(k,k',oo) ira 2 em E WQv\ 



cos 9=cos 9q 



oIlo v(k,k') \k + k'\ 

where is the step function and where 



i?(l-|cos0 o |) (4.1.10) 



( fc + A/)2+ W 2_( e + e /_ w )2 

cose ° = MkTV) • (4 ' L11) 



1 Since the time t in Eq. (I4.1.3P is fixed, I assume to work in a reference system which 
has the z axis aiong the instantaneous direction of the magnetic field B/(t). 
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At this point I have all the ingredients to calculate the photon energy 
spectrum per unit time by means of Eq. (|4.1.3p Fl The region of integration 
in Eq. (|4.1.3|) can be divided into the four sectors: 



k > 0, k' > 0, 

k > 0, k' < 0, 

k < 0, k' > 0, 

k < 0, k' < 



(4.1.12a) 
(4.1.12b) 
(4.1.12c) 
(4.1. 12d) 



and it can easily be seen that the integrals in the sectors I and II are equal 
to those in the sectors IV and III respectively. Also, the time dependence of 
the spectrum (|4,1.3p is carried by the oscillating functions in the distribution 
f(k,t) [see Eq. (|4.1.2j) ]. Since I am interested in macroscopic times such 
that mt S> 1, I can calculate directly the mean value of Eq. (|4.1.3p f° r 
large times by means of the substitutions (sin 2 e(£;)i) = (sin 2 e(k')t) ~ 1/2. 
Finally, by using the adimensional variables r\ = m/e = m/e(k) and rj = 
m/e' = m/e(k') that vary in the unit square, Eq. (|4.1.3p can be written in 
the form 



/ dN^(uj,t) 
\ dujdVdt 



a 2 em m 3 ( nR ±M y (B r \ 6 
V 2vr J \B cr ) 

l M(rj ) r/,w ) +l)+M(T),r/,w,-l) 



2vr 3 



o J o 



dr\dr\ 



V(l-^)(1-^' 2 ) 



(4.1.13) 



where w = tu/m and 
M(rj,r]',w,C) -- 



(W) 



A 6 



&(1 - | cos9o(rj,r]',w,C)\) 



x'(r],7]',w,() 
1 



+ 



1 

+ 



+ 4 



x(ri,7]',w,C) x'(ri,r)' ',w,() 

1 1 

+ 



with 



2w 



1 - \fl-rj 2 cos 6 (r),r)',w,C) 



(4.1.14) 



(4.1.15a) 



2 Note that the integral in Eq. (|4.1.3p seems to be divergent because of the factor 
l/|fc+fc'| in Eq. (|4.1.10p . but, actually, the constraint | cos 8q\ < 1 prevents this divergence. 
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2 W r 

x'(rj,ri',w,() = — 1 - C\A ~V' 2 cos 9 (t], 7]', w,() 

Tj I 



cos o (r),r)',w,C) 



w 



(V + rf) + CV!-^V!-V 2 



777/ 



it) 1 77 



(4.1.15b) 
(4.1.15c) 



The analytical integrations in Eq. (|4.1.13p are very difficult because the 
trivial condition | cos 6o(r/, rf, w, C)| < 1 is a complicated condition on the 
integration domain over r\ and rf , so I shall present the result of a numerical 
calculation (see Fig. UTTj) . Although, the quantity {dN^ nn \uj,t)/{duodVdt)) 




Figure 4.1: Photon spectrum per unit time (^dN^ nn \ui,t)/(dujdVdt)^ in 
arbitrary units. The dotted curve represents a function proportional to 

is a spectrum per unit time I want to compare it at least qualitatively with 
the two typical GRBs energy spectra shown in Fig. 11.41 Unfortunately, 
the shape of the theoretical spectrum ()4.1.13f) results also qualitatively very 
different from those of GRBs (see Fig. \1A\) and then I can conclude that the 
pair annihilation mechanism can not be the dominant one giving rise to a 
GRB. Nevertheless, similarly to experimental GRBs energy spectra, it shows 
a break energy w[ ann ^ of the order of 0.1 MeV (see the end of Par. I1.2.ip . 

Actually, the value of cj[ ann ^ is precisely around the value 0.51 MeV of the 
electron mass and this result can be interpreted in the following way. In 
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the center-of-momentum system the pair will yield two photons with equal 
energy and thus with oj > m. In order to get a soft photon one needs a 
large boost of the pair system and this would unavoidably result in a very 
energetic electron (or positron), but the distribution f(k, t) decreases rapidly 
for large k and so this process has small probability [see Eq. (I4.1,2p ], For the 
same reason the production of very energetic photons is also unlikely and 
the high-energy region of the spectrum in Fig. 14.11 decreases more rapidly 
than uj~ . I point out that the comparison with a function proportional 
to lu~ 3 is motivated because, as we have seen at the end of Par. 11.2.11 
the high-energy part of the experimental GRBs spectra are well fitted by a 
function proportional to u>~^ with j3 ~ 2-3. Finally, I want to observe that 
the sharpness of the peak at 0.51 MeV is also due to the fact that only the 
TGSs have been put into the electron and positron distributions f(k,t). In 
fact, the other possible states correspond to excited Landau levels and the 
energy of the electron (or of the positron) in such states has in the center- 
of-momentum system a minimum which is higher than m. Consequently, 
taking into account these states would make the maximum in Fig. 14.11 less 
sharp. 

As I have mentioned before, the production of photons through pair an- 
nihilation has as a consequence, of course, the depletion of the electron and 
positron populations. In fact, even if the two processes of pair production 
and pair annihilation are not disjoint in time, I performed the calculation 
keeping them separated because, in the presence of a purely rotating mag- 
netic field, the electron (positron) population does not grow continuously but 
it reaches a stationary density with superimposed rapid fluctuations [see Eq. 
(|4.1.2p ]. Now, the electron population is the same as the positron popula- 
tion and this equality is clearly kept by the annihilation process. However, 
the rate of change of these populations depends in general on k and the 
time evolution may be very complicated^ An overall indication of the time 
variation may be obtained by assuming that the momentum distributions 
of the electron and positron populations do not depend on time. In this 
approximation the electron (positron) distribution, that I also indicate as 
f(k,t), can be written in the factorized form 

f(k,t) = K(k)T(t). (4.1.16) 

Now, it can always assumed that J dkK(k) = 1 in such a way the function 
T(t) is given by 

T(t) = j dkf(k,t). (4.1.17) 

Also, I observe that the photons production rate is twice with respect to the 
rate of the electrons annihilation, then, introducing the quantity cr(k, k') = 

3 For example, in this case one should also take into account the scattering processes 
among the electrons and the positrons that do not change the number of particles but 
that do change their momenta. 
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J cLudcr(k,k' ,u))/cLu that is the total cross section of the annihilation of an 
electron and a positron with energies e(k) and s(k') respectively into two 
photons, I get from Eq. (14.1.30 

dN^ nn \cu,t) 



t {t) = /**/(*,*) = -\ J ^ = - 2T '^ < 4L18 > 



with 



dkdk'a(k, k')v(k, k')K{k)K{k'). (4.1.19) 
The solution of the differential equation (I4.1.18P is 



T{t) = 1 7° - (4.1.20) 
1 + s z T t 

with To = T(0), so I can conclude that the photon production rate —2T(t) 
decreases at large t as t~ 2 . It is not difficult to verify that, in the same 
hypotheses, if I started with unbalanced populations (more electrons than 
positrons, or the opposite), then the photons production rate would show 
an exponential decay with time. 



4.2 Photon production as synchrotron emission 

As I have said in Par. 11.2.21 there are strong evidences that the photons 
forming a GRB are emitted as synchrotron radiation by the electrons and 
the positrons in the fireball. I have also said that the fact that the gamma- 
ray radiation of the GRBs is high linearly polarized leads to think that the 
magnetic field in which the radiation itself is emitted is produced by the 
central engine of the GRB that is by a magnetar or a black hole. For this 
reason, I want to calculate here the energy spectrum of the photons emitted 
as synchrotron radiation by electrons and positrons in the presence of the 
rotating magnetic field ([4.1. ip [94j . As in the previous Section, the choice 
of the rotating field configuration ensures an easier mathematical treatment 
but it is also justified because the production of electron-positron pairs (and 
then of photons) is much more efficient with respect to the production in the 
presence of a magnetic field varying only in strength. By contrast, in the case 
of photon production through pair annihilation I have used the electron and 
positron distributions previously calculated to obtain, by means of the pair 
annihilation cross section, the photon spectrum per unit time. In this case 
it is easier to start from the beginning and to calculate the matrix elements 
corresponding to the whole process "pair production+photon emission" . 



4.2.1 Theoretical model 

By considering the process I want to study, a good theoretical starting point 
is the Lagrangean density Jz? of QED in the presence of an external electro- 
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magnetic field. If Ap a (t, r) = [0, — A/(i)] with 



A,(f,r) = ~[rxB ; (t)] (4.2.1) 
is the fourpotential describing the external rotating magnetic field then 

(4.2.2) 



1. 



2 = rp {^[ida + eA; a (t, r) + eA a ] -m}%b- ^ a p¥ a ^ 

where the radiation field A a (t,r) = [V(t, r), — A(t, r)] is assumed in the 
Coulomb gauge 



and where 



V(t,r)=0, 
d A(t,r) = 

F^i, r) = d a Ap(t, r) - dpAait, r). 



(4.2.3a) 
(4.2.3b) 

(4.2.4) 



In Eq. (j4.2.2H the two terms proportional to Ff a pFp and to F^ a pY a ^ with 
Ff a p(t,r) = d a Afp(t,r) — dpAf Q (t,r) have been omitted because they do 
not give any significant contribution to the equations of motion of the fields 
ip(t,r) and A(t,r) = [A°(t, r), . . . , A 3 (t, r)] and to the process I want to 
study. 

Now, the Lagrangean density (I4.2.2P is not in the most suitable form 
for the physical scenario I want to describe. In fact, as it stands it would 
be suitable for dealing with a weak external magnetic field because the 
usual perturbation theory could be used, while I am dealing with a strong 
magnetic field. Nevertheless, we also know that the perturbation induced 
by the macroscopic magnetic field (|4.1.ip can be assumed to be adiabatic. 
For this reason, I first perform the time-dependent rotation 

r'(t) = [x (t),y'(t),z (t)] = (x,ycosClt — zsmQt,ysinQt + zcosQt). (4.2.5) 

As a consequence, the spinor field ip(t, r) and the fourvector field A(t, r) 
transform as 



ip'(t,r'(t)) = exp (-»yfit) 1>{t,r) 
A'(t, r'(t)) = exp (-iS x M) A(t, r) 

where the matrices a x and S x are given by 



(4.2.6a) 
(4.2.6b) 



0", 



/0 1 0\ 

10 

1 

\0 1 0/ 



Sx 











°\ 























i 


Vo 





— i 


0/ 



(4.2.7) 
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I point out that, although a x and S x are two 4x4 matrices, they act on 
two different spaces: the first one acts on the spinor space and the second 
one acts on the fourvector space labeled by the Lorentz indices {0, . . . , 3}. 
Now, by using the equality 

exp (-i^-ntj = cos y^) ~ icr x sin (jpj (4.2.8) 

and by exploiting the usual commutation rules among the Dirac 7 matrices, 
it can easily be shown that 

exp y— i-^-QtJ 7 exp yi-^-Qtj = exp (iS x Qt) 7 (4.2.9) 

with 7 = (7 , . . . ,7 3 ). This equation is nothing but the mathematical ex- 
pression of the fact that the current density ip(t, r)7^(t, r) transforms under 
the rotation (|4.2.5|> as a fourvector [see Eq. (|4.2.6b|) ]: 

■ijjjip = if? exp f-i-^-nt) 7 exp (i-^Qi) i)J = exp (iS x Qt) tp'^yip' (4.2.10) 



where, from now on, the "primed" fields are always intended to be calculated 
at the "primed" coordinates r'(t). 

These previous equations can be exploited to rewrite the Lagrangean 
density (|4.2.2|) in terms of the primed variables and fields. In particular, it 
is evident from Eqs. (14.2. 6b j) and (|4.2.1()l) that 

^ a ^A a = i>'j a 4>'A' a . (4.2.11) 

The transformation of the terms involving the external magnetic field are 
more complicated. In fact, from Eqs. (|4.1.1|) . (|4.2.1|) and (|4.2.9|) I obtain 

ip 7 a M? a (t,r) = ?l [^y^y(t)-^yvv(i)] =^ 7 a ^4«(r'(t)) 

(4.2.12) 

where I have introduced the fourpotential A'^ Q (r'(t)) = [0, — Ay>(r'(i))] with 
A'^r'W) = - [r'(i) x B' ; ]/2 and B' ; = (0,0,B/). ' 

Now, I want to show in detail how the terms containing derivatives of 
the fields in the Lagrangean density (|4.2,2p transform. I first transform 
separately the time derivatives dtip(t,r) and dtA(t, r). From Eqs. (|4.2.6p 
and by reminding that the variables r'(i) depend on time [see Eq. (|4.2.5I) ]. 
I have 

d t iP(t, r) = exp (iy fit) [Wj 1/2) V (*, r'(t)) + d t ifj'(t, r'(i))] , (4.2.13a) 
d t A{t, r) = exp (iS x nt) \injW A' (t, r'(t)) + d t A'(t, r'(t))] . (4.2.13b) 
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where I have introduced the x components 

jj 1/2) '=4+y, (4.2.14a) 
JP' = 4 + S x (4.2.14b) 

of the one-particle electron (positron) and photon total angular momentum 
operators respectively. On the one hand, by means of Eq. (|4.2.13ap it can 
be seen that 

^ a 8 a ^ = 4>' exp (-iy Oi) 7% exp (i y fit) V' + $V Sj^' 

+ (7 2 cos Qt + 7 3 sin tit) (cos fit ^ + sin D,t d' 3 ) ip' 
+ ip' (-7 2 sin ft£ + 7 3 cos Qt) (- sin fit d' 2 + cos Sit d^) if;' 

(4.2.15) 

where, for notational simplicity, I do not indicate the dependence on time of 
the derivatives with respect to the primed coordinates. On the other hand, 
Eq. (14.2. 13bj ) can be used to see how the free radiation field Lagrangean 
density in Eq. (|4.2,2p transforms. Firstly, I write it in the more useful 
noncovariant form 

" = \ [ (dtA) 2 - (d x A) 2 ] . (4.2.16) 

Now, it is evident that 

(d x A) 2 = (d' x A') 2 . (4.2.17) 

Also, in the Coulomb gauge where the scalar potential vanishes identically 
one has that (d t A) 2 = -(d t A a )(d t A a ), then, by using Eq. (|4.2.13hj) . I 
obtain 

(d t A) 2 = (d t A'f + 2in(d t A') ■ jpA' + o(n 2 ) (4.2.18) 

where, being in my approximations the rotational frequency fi a small quan- 
tity, I have neglected the terms proportional to O 2 . By substituting Eqs. 
(|4.2.17|) and (14.2.180 in Eq. (14.2.161) I can write it as 

- = -\^' aP - in(d t A f a )jP'A"* + o(tf) (4.2.19) 

with 

P aP (t,r'(t)) = d' a A'^v'(t)) - & p M a {t,v'{t)). (4.2.20) 

By collecting Eqs. (I4.111|) . (j4.2M2j) . (|4.115l) and (14.2M9j) and by per- 
forming the remaining trivial transformation of the mass term in the La- 
grangean density (|4.2.2p , it can be written in terms of the primed variables 
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If' P a P 



(4.2.21) 



and fields as 

jSf = ^ {y + e4 Q (r'(t)) + eA' Q ] - m) $ 

+ in^jjV^y - ^(dtA'JJ-W'A' + 0(ft 2 ) 
or, by removing the now useless prime and the time dependence on r'(i)0 



as 



with 



and 



JSJg = $ {7 a [id a + eA' ;a (r)] - m} $ - \V a ^ 



(4.2.22) 
(4.2.23) 



££[ = e^j a ^'A' a + infi^jW 2 ^' - i^l{d t M a )J^M a + 0(ft 2 ). (4.2.24) 

In this way the original time-dependent Lagrangean density ()4.2.2p has been 
transformed into an effective Lagrangean density that does not depend ex- 
plicitly on time and that embodies the effects of the rotation of the ex- 
ternal magnetic field in the interaction terms proportional to its rotational 
frequency f2. I note that the Lagrangean density (|4.2.22|) is just the La- 
grangean density of QED in the presence of the external static magnetic 
field = (0,0, Bf ) plus other extra interaction terms that are propor- 
tional to O (or to 2 )H 

In order to build the Hamiltonian density I calculate now the momenta 
conjugated to the Dirac and to the radiation field. From Eqs. (|4.2.22|i - 
(I4.2.24P I obtain 



7r' A ,(t,r) 



t yt(t,r), 
" ir -d t A'(t,r)+mj^A'{t,r) 



d(d t A'] 

and the Hamiltonian density can be written in the form 
X s = n'^' + 7r' A , • d t A' - Xfi = + 

with 

jr ' = ^ {a ■ [-id + eA' r (r)] + [3m} $ + \ 



tv% + (dx A') 



(4.2.25a) 
(4.2.25b) 

(4.2.26) 
(4.2.27) 



4 The elimination of the time-dependence on the variables r'(t) can be safely done 
because it does not change the equation of motion of the radiation field. In fact, on 
the one hand, being the transformation r'(t) = r'(t, r) a rotation, then dr = dr'(t). On 
the other hand, the integral of the Lagrangean density Jzf' on r'(t) to obtain the total 
Lagrangean (and then the action) extends over the whole space and the coordinates r'(t) 
are dumb variables. 

As in the rest of the thesis, I continue to use the "prime" to indicate all the fields and 
the related quantities when the external magnetic field lies in the z direction. 
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and 

3% = -enj/^aif/ ■ A' + ifi^* J-j 1/2) - iflw'^ ■ J^A' + 0(n 2 ). (4.2.28) 

Now, since the perturbation induced by the magnetic field B^(£) is adi- 
abatic, I can consider all the terms in the interaction Hamiltonian density 
(j4.2.28jl as small perturbations of the free Hamiltonian density (|4.2.27jl . At 
this point, I can use the machinery of the ordinary perturbation theory to 
calculate the matrix elements corresponding to the process under study: 
the emission of photons by the electrons and the positrons created in the 
strong slowly-rotating magnetic field B^(i). By neglecting all the radiative 
corrections and taking into account only the tree-level contributions, the 
Feynman diagrams accounting for the mentioned process are those shown 
in Fig. 14.21 where the lower vertices actually represent the interaction with 
the time derivative of the external magnetic field and correspond to the 
term iQip 1 in Jtfj. Now, as in ordinary QED, in order to calculate 





Figure 4.2: Tree-level Feynman diagrams of the photon emission by an elec- 
tron [part (a)] or by a positron [part (b)] created in the presence of the strong, 
uniform and slowly-rotating magnetic field B^(i) given in Eq. ([4.1.1)) . The 
thick fermion lines indicate that the calculations of the corresponding S- 
matrix elements have been performed by using the electron and positron 
one-particle states in the presence of the magnetic field. 



the S'-matrix elements corresponding to the Feynman diagrams in Fig. 
I quantize the Dirac field and the photon field in the interaction picture. 
Since, as I have said, the Lagrangean density (|4.2.23|) is the free Lagrangean 
density of QED in the presence of a uniform and static magnetic field in the 
z direction with strength Bf, we already know that the Dirac field can be 
expanded as [see Eqs. (J3T7j) and (12.1.23|) ] 

ip'(t,r) = ^2 [ c i n j'( r ) ex P( — iwjt) + djVj(r) exp(ivjjt) . (4.2.29) 

3 
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Now, we have seen in the previous Chapter that pairs are much more 
likely created in the presence of a strong and slowly-rotating magnetic field 
with both the electron and the positron in a TGS. Analogously, I assume 
here that all the electrons and the positrons entering the game are in TGSs 
[see Eqs. (j2.1.35j) ]. In the same approximation I sum only on the electron 
and positron TGSs to build the propagator G'(t, r, t', r') that is [96] 

iG'(t, r, t', r') = Y,Mt ~ 0<fc(r)<fc(O ™.p[-ie k (t - t')\ 

n,k (4.2.30) 
- - t)v' n k (r)v n k (r') exp[ie k (t - t')]} 

where the coordinates r' have obviously nothing to do with those introduced 
in Eq. (I4.2.5P that, on the other hand, depended on time. 

I pass now to the second quantization of the radiation field. The presence 
in the interaction Lagrangean density (|4.2.24|) of terms containing the time 
derivative of the radiation field would make the quantization procedure very 
complicated. Nevertheless, I observe that 

1. these additional terms are proportional to the rotational frequency fi; 

2. the matrix elements that I will calculate are already proportional to Q 
through the factor corresponding in Fig. 14.21 to the interaction vertex 
with the time derivative of the external magnetic field. 

For these reasons, since I am not interested in higher order corrections in 
all the other factors in the matrix elements can be evaluated neglecting the 
interaction with the external field. In this way, I can quantize the radiation 
field as it were free and then I can indicate the vector radiation field simply 
as A(t, r) and I have only to expand it into the usual plane- wave basis as 

A(t, r) = -j^= {a qiA exp[-»(wt - q • r)] + a qA exp[i(wt - q • r)]} 

q,A V 10 

(4.2.31) 

where V is the quantization volume, to = |q| is the photon energy and e qi A 
with A = 1, 2 are the polarization versors [97]. 

At this point I have all the quantities I need to calculate the matrix 
elements corresponding to the Feynman diagrams in Fig. 14.21 and this is the 
subject of the next Paragraph. 

4.2.2 Calculation of the photon energy spectrum 

By looking at the interaction Hamiltonian density Jfj it is clear that if 
the final state is the state \e~e +/ y) = \ko, no; k' , n^; q, A) (the initial state is 
obviously the vacuum |0)) then the matrix element at time t corresponding 
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to the Feynman diagram in Fig. 14.2( a) can be written as. 

, fco,no,A;g,n ( ' ) ,q,A( 



SLnnJ*. n>. a.x(t) = I & \ dt' f dr" T dt" 







x n no,fc ( r/ ) aq>(*gfcoO 6 ^2^ A ^P^^*' ~ q ' r ')] 

x r', f, r'')^ 7 ° J. (1/2 %^(0 exp(ie fe ,i") 

(4.2.32) 



where I have pointed out that, while the electromagnetic interaction between 
the Dirac field and the radiation field is always present, the external field 
starts rotating at an arbitrary finite time t" I called zeroQ 

Now, the term in Eq. (|4.2.32|) corresponding to the lower vertex in Fig. 
14.2( a) will be calculated by means of the first-order adiabatic perturbation 
theory. In order to do this, I use the expression (|4,2.30p of the electron 
propagator to write the previous matrix element in the more useful form 

Sk ,n ,k' ,n' ,ci,\(t) 

= - f ^ L £ /I dt ' exp ^ fc o +L °- £k ^ 

x J dr'ul ) ko (r')a z u n k (r')exp(-iq- r') 

X Jo dt " e ^^ £ k'o + £ k) t "] J ^''^k^' y) d t" V n' ,k' {t" 

J dt 1 exp [i (e ko + u> + e k ) t'] 

J — oo 

x J dr / n / n t ofco (r / )a z ^ )fe (r / )exp(-iq- r) 

x j\t"eMKe K -e k )t"\ J dr'^t" ,r")d t ,v n>o K {t" ,v") 

(4.2.33) 

where I used the fact that the TGSs are eigenstates of a z in such a way a x 
and a y can not couple two of them and where I introduced the "rotating" 
states [see Eqs. (|2. 1.461) ] 

Jt, r) = exp(iOt J^)u' fe (r), (4.2.34a) 



6 The Feynman diagram in Fig. I4.2f b) represents the emission of the photon by a 
positron but this process can be taken into account by simply multiplying by two the final 
spectrum of the photons emitted only by an electron. 

7 I remind (see Sects. 11.2.31 and I3.1.ip that the emission of a GRB is supposed to 
happen during the collapse of a magnetar into a black hole or during the formation of 
the magnetar itself. In this framework the instant t" — in Eq. (|4.2.32|l indicates the 
beginning of the magnetar collapse or of its formation. 
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«„,*(*, r) = exp(iQtJ^)v' n>h (r) (4.2.34b) 

that are the instantaneous eigenstates at time t of the one-particle Hamil- 
tonian 

H(t) = a-[V + eA r (t,r))+/3m (4.2.35) 

with Af(t,r) given in Eq. (|4.2.ip . In this way, by applying Eq. (]2.2.8bj) of 
the adiabatic perturbation theory I obtain 



J dr " n i,fc(^> r ")^''^,^(t",r") = -J dr" 



£k' + £fc 



(4.2.36a) 



<fr " u i,fc(*"» r/ ')^t"WnJ ) *s(*"» r ") = - / dr " . __ 

(4.2.36b) 

These kinds of matrix elements can be easily calculated and the result is 

/j // t fj.ll H\a (4.H H\ &£lBf [ -r II 'j I ll\ X a z l f //n 

dr <fc(* > r )dt"V n > oi k> {t ,r ) = ~ £k+£ki J dr <fc( r )^r v n' ,k> ( r ), 

(4.2.37a) 

/j // t ('-1.'/ "\ Q /j." "\ e^lBj- f r II If I ll\ X a z I , ft\ 

dr v^ k {t ,r )d t iiv n i ojK (t ,r ) = - — y v^r )^-v <M) {r ). 

(4.2.37b) 

I observe that in both these matrix elements the integral on the z variable 
gives a conservation of the longitudinal momentum and then of the energy. 
This does not cause any problem in the first matrix element, while the second 
one diverges when k = k' Q . For this reason, this particular matrix element 
will be calculated by writing the left-hand side of Eq. (|4.2.37bp as [see Eq. 
fl4.2.34hjt ] 

dr''< fc ,(t" ! r")^^ !& ,(t" ) r") = -^ / d^u^)^^(f). 

(4.2.38) 

By substituting the explicit expression of the one-particle electron total an- 
gular momentum (I4.2.14ap . I observe that, on the one hand, the term z"Vy 
does not contribute because, by performing the integral on z" from —L z /2 to 
L z /2, it vanishes and, on the other hand, the term a x /2 does not contribute 
either because the TGSs are eigenstates of a z . As a result, we have 



dr' 



(4.2.39) 
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At this point, if I substitute the expressions (I2.1.35P of the TGSs I obtain 
that the matrix elements different from zero are the following ones: 

^<-W^>^*^ = ^y^, (4.2.40a) 

A " t (+11 H\X (+» »\ ^ l eB r( n + 1) /„ n , n ,^ 



-a-' 



f/r 



/ ^'^ +1 , fc , (t",r")^^, fc , (t",r") = "W^- (4.2.40d) 

By inserting the previous matrix elements in Eq. (|4.2.33p and by omit- 
ting the now useless index "0" on ko, k' Q , uq and n' , I obtain the two tran- 
sition amplitudes 

<v(i) (f) 

x / dt' exp[i(e fc — ey + oo)t') \ dt" exp(2ie fc /i") 
J-oo Jo 

e£lk (e q x) z In' f , /+ . , . . . 

J dr <fe( r )«^n'-i,fc'( r )exp(-iq-r ) 



+ 



eBtVuj 



[ dt'ex.p[i(e k + e k < + uj)t'] f dt" 

J-oo Jf 

(4.2.41) 



and 

v(2) 



= -^^y^^5/*'< t (r') a ,<, +1 ,^,(r')exp(- iq .r') 

x / dt' exp[i(efc — + uj)t'] / dt" ex.p(2iek't") 
J-oo Jo 

e»fc / (e qiA ) z I n' + 1 f , n , , 
x I dt' ex.p[i(s k +e k > +oj)t'] f dt" . 

J-oo Jt' 

(4.2.42) 
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Now, by using the expressions (|2. 1.35f) of the TGSs it is easy to show that 

J dr ' u n,k( r ') a ^n'±l,-k'( r ') exp(-iq ■ r) 

(e k + m)(ey + m) ( k k' \ , 

- J 1 n,n'±l,q x ,q y k+q z +k> ,0i 



4:Ek£k' V e fe+ m £ k' + m 



(4.2.43) 

' ) I'n,n'±l,q x ,q v 5k+q z +k> fl 



(e k + m)(ey + m) k k' 



4e fc e fc / V e k + me k > + m j 

(4.2.44) 

where 

4n', 5l)ft = / dxdye'*(x,y)9' n ,(x,y)exp[-i(q x x + q y y)] (4.2.45) 

with the functions 6' n (x,y) given in Eq. (|2.1.36|) . 

I proceed now by calculating the time integrals in Eqs. (|4.2.4ip and 
(I4.2.42p . The results are 



/d£ exp[i(efc + uj — ey)t'] / dt" ex.p(2ieyt") 
oo JO 



1 J exp[i(e k + ey + to + is)t] exp[i(e k — ey + uj + is)t] 
2iey \ i(e k + ey + uj + is) i(e k — ey + uj + is) 

(4.2.46) 



and 

/ ,/ r/ n /i f f , » exp[i(ei, + si./ + w — is)il 
dt'expi {e k + u + ey)t'\ / dt" = * -/i 4.2.47 
-oo 7t' ( £ fc +Sy +U) - is) z 

respectively and the is terms with s — ► + have been added in order to 
make the integrals convergent. Now, it is obvious that e k + ey + uj > 0. 
Also, because of the overall conservation k + q z + k' = of the longitudinal 
momentum coming from Eqs. (I4.2.43P and (|4.2.44jl then, unless the trivial 
case q = 0, it can be shown that e k — ey + uj > 0. In this way, all the 
is terms can be safely eliminated in the final results in Eqs. (|4.2.46p and 
(|4.2.47p . Finally, by substituting the just calculated time integrals and Eqs. 
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(14.2.4;-! I) and (14.2.441) in Eqs. (|4.2.41jl and (|4.2.42jl . I obtain 



Q /(i) m _ e ^( e g,A)z I (e k + m){e k > + m)ri . 

eBj-m ( k k' \ exp[i(£ k + ey + 



8e|, V £ fc + m Ek'+mJ Ek + Ek'+uj 
+ k'exp[i(£ k + E k i +u)t] f l+ k k' 



(£ k +E k i+uj) 2 \ £ k +me k '+m 
eBjm ( k k' \ exp[i(£ k — £ k i + uS) t] 



8e|, \E k + m £ k '+mJ £ k -E k i+u 

(4.2.48) 



and 



<v(2) m _ eO(e q , A ) z / (£fc + m)(e fc / +m){ri + 1) , 

d fc,«,fc',n',q,AW - 4 y efcefc/ewS/F V»'+i, to ,9 t , <) *+?.+*',o 

eBf-m ( k k' \ exp[i(£ k + £ k > + uj)t] 



8e|, \E k + m E k ' + mJ E k +E k i+uj 
k' exp[i(£ k + £ k > + uj)t] f ^ k k! 

(£ k + e k > + to) 2 V e k + m e k t + m 
eBj-m ( k k! \ exp[i(£ k — £ k > + u>)t] 



8e|, \£ k + m £ k i+mJ £ k - e k , + uj 

(4.2.49) 



The probability that a photon is emitted at time t with momentum between 
q and q + dq by an electron or by a positron is given by 



2 00 

y > J J !,A=ln,n'=0 

(4.2.50) 

where the limit of large L z and V is understood and all the momenta are 
intended from now on to be continuous variables. As usual, I am interested in 
macroscopic times t such that mt 3> 1 then I can neglect the oscillating terms 
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coming from the square modulus of SC. \ x (k, k', q, t) and S'jfl, Ah, k' , q, t): 



dp( syn )(q,t) 



Z/ <iq e<? 2 ^ 2 °° 

X] [ n/ | J n,n'-lfe>9 2 /)| 2 + (n / + l)|/n,n'+lfe I 



(2vr) 4 4uj 3 B 
x dk 



n,n'=0 

e(k) + m] [e(k + g 2 ) + m] 



e(k)e(k + & 
eB+m 



8e 3 (k + q z 



k k + q z 



1 2 



e{k) + m e(k + q z ) + m 



[e(k) +uj} 2 + e 2 (k + q, 
[[e{k)+oof -e 2 (k + q z 

{k + q z f 



[e(k) + e(k + q z ) + loY- 



k 



k + q z 



1 2 



e(k) + m e(k + q z ) + m 

(4.2.51) 



where I exploited the longitudinal momentum conservation to perform the 
integration on k' and where I substituted [97] 



ElMq)],| 2 = i-4 = % 

A=l 



(4.2.52) 



with = Qx + 9y- Concerning the sums on n and n', I will calculate them 
together with the integrals I' nn i±i(Qx,Qy)- In fact, 



E Wn ,n'-l 

(q x ,q y )\ 2 + (n' + l)\l' n 

,n'+l 

n,n'=0 

00 

= E (2n' + l)|/; in ,( M )| 2 . 

n,n'=0 

Now, from Eqs. (|4.2.45|) and (|2.1.36|) I can write I' nn i(Qx,Qy) as 
1 



(4.2.53) 



x exp 



where the change of variable 



(4.2.54) 



(4.2.55a) 
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V 



eBj 



-y 



(4.2.55b) 



has been performed. With this expression I can calculate explicitly the sums 
on n and n' , in fact 



n,n'=0 



£ (2n' + l)\l' n)n ,{q x ,q y )\ 2 

j dtdridgdr/ exp[-(£ 2 + if + £' 2 + r]' 2 )] 
x ex P \ HT^tex^' -0+ Qy(v' ~ V)\ \ 



Z — «■ n \ L — ✓ n '\ 

n=0 n'=0 

i j dtdrjdgdn' exp[-(£' - if - (rf - V f} 

x exp J i JL[ qx ^> qy{r] ' _ v)] I [ 2(e _ ir) )g + ir f) + i]. 



If now I put 



v± 



7/ ± ry 



(4.2.56) 

(4.2.57a) 
(4.2.57b) 



I obtain 



oo „ 

£ (2n' + l)K x fe,rf=^ / ^+^-^+^-exp[-2(e 2 +t? 2 )] 



n,n'=0 



x exp 



x [(e+ + ^-) 2 -(^ + ^+) 2 + i]- 

(4.2.58) 

As I have pointed out before, the presence of the external nonuniform elec- 
tric field E^(i,r) = — d t Af(t,r) [see Eq. (|4.2.1|) ] breaks the translational 
symmetry in the plane perpendicular to the magnetic field. In particular, 
as we have seen in the previous Chapter, this fact also made the presence 
probabilities per unit volume diverging in the regions far from the origin 
[see Eqs. (|3.2.2p . (|3.2.3ip and f|3.2.45|) ] . In the present case the divergence 



76 



Section 4.2 



comes from the integrals on the variables £+ and 77+ in Eq. (I4.2,58p . For 
this reason I will keep in (|4.2.58|) only the dominant terms in these integrals 
that is 

00 

£ {2n' + l)\l' n>n ,{q Xy q y )\ 2 ^- d V+ dt; + (e + +V 2 + ) 
ri'—n •* 



n,n'=0 



x / d£-dr]- exp 



2i 



(4.2.59) 



Now, by passing to polar coordinates in the plane, I easily obtain [see 

Eqs. ([4.2.550 and ([4.2.57D ] 



dv+d£+(£+ + vl 



7r / eB 



-R±m 



(4.2.60) 



where R±m is the transverse radius already introduced in Eq. (|3.2.42 j) . 
Instead, the integrals on the variables £_ and r/_ are well-known exponential 
integrals and I only quote the final result: 

E {2n' + l)\l' n ^%)\ 2 ~~ {^-) 2 R\m^V {^^t)j ■ (4.2.61) 

By substituting Eqs. (14.2. 52p and (|4.2.61j> in Eq. (|4.2.51l) I obtain the 
following expression of the probability (dP( syn )(q, t)): 

<J(iP (syn) (q,t) 



2eB r n 2 a em ql y dq T 



00" 



L Z R ±M exp 



7 2 

dxy 



2eBi 



x I dk 



[e(k) + m] [e(k + q z ) + m] 



e(k)e(k + 1 



eBj-m 



8e 3 {k + q z ) 



k 



+ 



k + q z 



e(k) + m e(k + q z ) + m 



[e(k)+u] 2 +e 2 (k + q z ) 



[[e(k) + 



wr — £' 



+ 



'■(k + qz)]' 



2[e{k) + e{k + q z ) + uj\ 



k 



k + q z 



e(k) + m e(k + q z ) + m 



(4.2.62) 



Finally, the photon spectrum per unit volume V = L z irR 2 lM is obtained 
by passing to photon momentum spherical coordinates {to, 6, (ft} and by in- 
tegrating on the angular variables. Only the integral on the variable (ft is 
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trivial then by putting u = cos 9 I obtain 



' 'dN^{u,t)\ a em wm 2 (nR ±M ) 2 B ; 



dudV 
-l 



x 



J du(l — u 2 ) exp 



(47T) 
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cr 2 



^(1-n 2 ) 



x / cfc- 
'o 



2B r 

1 + e(v))[l + e(v + uw)] 



'Bj_ 
8B C 



e(v)e(v + uw) 

2 

r 

+ 



V + W 



l + e(u) 1 + e(v + uw) 

1 [e(u) + ro] 2 + e 2 (u + w) 

e 6 (v + uw) [[ e ( v ) + u,]2 _ e 2( w + : 

(v + uw) 2 
-I ^ 

2[e(v) + e(v + uw) + w}' 



V + uw 



1 2 



1 + e(v) 1 + e(v + uw) 



(4.2.63) 

where I introduced the adimensional quantities w = uj/m, v = k/m and 
e(v) = VI + v 2 . 



Now, the integrals in Eq. (j4.2.63H can not be performed analytically 
and, for this reason, I resort to a numerical integration. Fig. 14.31 shows 
the photon spectrum (j4.2.63H in arbitrary units and with a magnetic field 
strength Bj- = 2.2 x 10 14 gauss. The qualitative form of this spectrum is 
very similar to those shown in Fig. 11.41 In fact, as it is evident from the 
figure, the spectrum shows two different behaviours below and above the 



break energy u 



(syn) 



1-3 MeV. I have checked numerically that, unlike the 



annihilation photon spectrum (see Fig. 14. ID where lu^ ~ m independently 
of the external magnetic field strength, the value of here depends on 

Bf and, in particular, the lower is Bf the lower is uj^ yn \ Nevertheless, it is 
important that by using a typical magnetar magnetic field strength such as 
Bf = 2.2 x 10 14 gauss, then the value of o>^ syn ^ is close to the break energies 
characterizing the experimental GRBs spectra. In fact, as I have said in 
Par. 11.2.11 the experimental break energies are typically below 1 MeV, but 
there are also cases of GRBs with ujb > 1 MeV. On the other hand, I have 
also said in Par. 11.2.11 [see Eq. (jl.2.ip and below] that the experimental 
spectra of GRBs are very well fitted by a function proportional to cu^ 1 in 
the low-energy region and by a function proportional to ui~@ with (3 ~ 2-3 in 
the high-energy region. In the present case, we see from the figure that the 
high-energy part of the spectrum decreases more rapidly than w~ 3 . This is 
be due to the fact that, as I have already pointed out in the previous Section, 
the production of high-energy pairs due to a slowly-rotating magnetic field 
is disfavored [see discussion below Eq. (|4.1.2p ]. Instead, concerning the low- 
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Figure 4.3: Photon spectrum ^eL/V( syn ) (lj, t) / {doodV)} in arbitrary units. The 
magnetic field strength Bf is equal to 2.2 x 10 14 gauss. The dotted curve 
represents a function proportional to u;~ 3 . 



energy region of the spectrum, I want to show analytically that the spectrum 
(I4.2.63P goes just as uj^ 1 in the limit w = u/m <C 1. In fact, all the terms 
in the integrals on u and v that are finite if calculated at w = give a linear 
dependence of the spectrum on the photon energy because of the presence 
of the overall factor proportional to w in Eq. H4.2.63j) . By keeping the only 
diverging term in the low-energy limit, Eq. (|4.2.63|) becomes 

/ dNto*)( u ,t) \ M< i 2 , op ,2 (Br 

\ dudV J ~ ^emWm(nRxM) 



"1 roc „,2 

x / du(l - u 2 ) I dv 



1 



o e 8 (v) [e(v) — e(v + uw) + w] 2 

(4.2.64) 



Now, I can write the diverging factor [e(v) — e(v + uw) +w] 2 in the form 

! i^tJ^ ^ 



[e(v) — e(v + uw) + w] 2 w 2 [e(y) 



vu 
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and, by substituting this expression in Eq. (|4.2.64l) . I finally have 



/ dN^ n \io,t) \ 
\ duodV / 




e 6 (v) [e(v) — vu] 2 
(4.2.66) 



v 2 1 



which is the desired result. In fact, since e(v) — vu > in the integration 
domain, the two integrals are finite and then 



4.3 Direct photon production from vacuum 

In the two previous Sections I have discussed two processes in which photons 
are produced by real electrons and positrons previously created in the pres- 
ence of the magnetic field (I4.1.ip . Instead, I want to study here a process in 
which the photons are still produced in the presence of the slowly-rotating 
magnetic field (|4.1.ip but directly from vacuum that is without an inter- 
mediate electromagnetic interaction with real electrons and positrons. As 
I have done previously, I will treat the dynamics of the system by means 
of the adiabatic perturbation theory. Since I want to analyze the direct 
production of photons, I have to build an Hamiltonian which describes the 
interaction between a quantized and a classical electromagnetic field. The 
use of the effective Lagrangian technique is particularly useful to this scope 
|13| [H] but, as a consequence, the results will be reliable only for photon 
energies much less than the electron mass mJl 

I start with a completely general expression of the effective Lagrangian 
density if of a free electromagnetic field [Et(£, r), By(i, r)]. It is known 
that, in order to be a true Lorentz scalar, the effective Lagrangian density 
must be a function only of the relativistic invariants Ft and G T (because 
Gt is, actually, a pseudoinvariant) with 




(4.2.67) 



F T = -^{B T - E T ) 




(4.3.1b) 



then, in general, 




(4.3.2) 



8 On the other hand, as I have said at the beginning of Sect. 14.11 in order to consider 
the external magnetic field as uniform and slowly-varying then the photon energies are 
always assumed to be much larger than Q. 
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Now, suppose that the total electromagnetic field is the sum of a radia- 
tion field [E(t, r), B(t, r)] which will be quantized and of a classical field 
[E(i, r), B(t, r)] which, at the moment, I assume to be the constant and 
uniform magnetic field B given in Eq. (|2. 1 .38[) : 

Er(t,r)=E(t,r), (4.3.3a) 
B T (t,r) = B(t,r) + B. (4.3.3b) 

Since I am not interested in the interactions of the radiation field with itself, 
I expand the effective Lagrangian density (|4.3.2p up to quadratic terms in 
E(i, r) and B(t, r) and I obtain the expression 

jgfp) = | (E 2 - B 2 ) + ~ [c 2 (n • E) 2 + c 3 (n • B) 2 ] (4.3.4) 

where n = B/5 and where 



ci 



dF T 



E=B=0 



r - B 2 Q2j£ 



E=B=0 



£=B=0 



(4.3.5a) 
(4.3.5b) 
(4.3.5c) 



In the Lagrangian density (|4.3.4p the interaction between the radiation field 
and the classical field is described by the last two terms. The strength of the 
interaction depends on the coefficients c 2 and C3 and I will give them later 
in the case in which the effective Lagrangian density is the Euler-Heisenberg 
Lagrangian density [981 EU E] . 

In order to build up the Hamiltonian density corresponding to the La- 
grangian density (|4.3.4p I have to introduce the fourpotential field A a (t, r) = 
\V(t, r), A(t, r)] relative to the radiation field. If I choose a gauge in which 
V(i, r) = 0|£| the radiation field is given by 

E(i,r) = -<9 t A(i,r), (4.3.6a) 
B(t,r) = 0xA(t,r). (4.3.6b) 

By expressing the Lagrangian density (j4.3.4H in terms of the vector potential 
A(i, r) and of its derivatives, I can calculate the momenta ir/\_{t, r) conju- 
gated to A(t, r) as 

a£?( 2 ) 

ir A (*,r) = = -cit(t,r) - c 2 [n • E(t, r)]n (4.3.7) 



9 It can be shown that another condition is needed to fix unambiguously the radiation 
field gauge [see [55] f° r a more detailed treatment]. 
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and, finally, the Hamiltonian density 
^ 2 W A -cVV-^ {2) 



— + — B 2 — — * -(n • vr A ) 2 - * (n • B) s 

2ci 2 2ci(ci+c 2 ) V A; 2 K ' 



(4.3.* 



At this point in order to calculate the energy spectrum of the pho- 
tons produced directly from vacuum in the presence of the slowly-rotating 
magnetic field (|4.1.1j) . I should substitute the static magnetic field B with 
Br(t) and to apply to the total time-dependent Hamiltonian H^ 2 \t) = 
J drJ4?( 2 ^ (t) the adiabatic perturbation theory. Since, as I have mentioned, 
the production through this mechanism is low I only quote the final result in 
the particular case of the Euler-Heisenberg Lagrangean density [98|. \99\ [5] : 

J^(E-H) = _J_ Ft 
Air 



[°°da 2 
/ — r exp( — im s) 
Jo ^ 



1 f°° ds 
8^2 



(es) 2 |G T |cot I es\j \IF 2 +G T + F T 



x coth ^F 2 + G\ - F T J - 1 + ^(es) 2 F T 

(4.3J9) 

where Ft and Gt are the relativistic invariants of a generic constant and 
uniform total electromagnetic field (Et,Bt). In our case, it can be shown 
that the following asymptotic estimates of the coefficients ci, c 2 and C3 hold 
for external magnetic field, strengths Bj B cr [see Eqs. (|t33H pl: 



(E-H) 


1 

47T 


"em 1 

- Q log 

37T 


(E-H) 

c 2 


"em 

3vr 




(E-H) 


"em 




c 3 


3vr 





(4.3.10b) 

(4.3.10c) 

(E— HI 

where I put approximatively c\ ~ l/(47r) since in every realistic physical 
situation the "log" term is always negligible. By reminding that = 137, I 
can assume that the magnetic field B^(i) is such that 1 <C B^/B cr <C l/a em 
(if Bf = 10 15 gauss then BfjB cr = 22.7). In this approximation it can be 
shown that the final energy spectrum of the photons produced directly from 
vacuum is given by 

/ dN( dh \uj,t)\ 28 n 2 (a em B r \ 2 B 



dudV / 135 (2vr) 3 V B. 



-) ifl«|^«a- m . (4.3.11) 



10 Note that the coefficients ci, C2 and C3 were defined starting from the Lagrangian 
density (I4.3.4|l where the magnetic field does not yet depend on time. 
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As it is evident, the photon energy spectrum does not depend on the photon 
energy in contrast with the experimental GRBs spectra that, in the low- 
energy region lu <^ m are proportional to u . Also, it can easily be checked 
that the direct production of photons from vacuum is low and completely 
negligible with respect to the production through synchrotron radiation (the 
comparison with the production through pair annihilation is harder because 
I have calculated not a photon spectrum but a photon spectrum per unit 
time). In fact, just to give an idea, by assuming Bf = 10 15 gauss and 
R±M = 10 5 cm [see the discussion below Eq. (|3.2.26|) ] then 



{dN( dir \u,t)/(dudV)) 



'dN^){uj,t)/{dujdV)) 



~a em -^(— M -2.oxl()- T2 . 

D cr \ K I 



, , i > ■ ■ ■ x -L M 

u=l MeV 



(4.3.12) 



4.4 Summary and conclusions 

In this Chapter I have presented an analysis of the production of photons 
in the presence of a strong, uniform and slowly-rotating magnetic field. As 
I have said, this investigation has been motivated by suggestions coming 
from the study of GRBs. Nevertheless, the theoretical attitude has been 
to consider some very simplified versions of the processes so that the phe- 
nomenological parameters and the dynamical details could be reduced to 
the minimum. In fact, the external parameters entering the game are only 
two: the magnetic field strength Bf and its rotational frequency 0. The 
analysis has yielded some definite results, making evident, in quantitative 
form, the presence of a production of photons 

1. through the annihilation of pairs previously created in the presence of 
the rotating magnetic field; 

2. as synchrotron radiation by electrons and positrons previously created 
in the presence of the rotating magnetic field; 

3. directly from vacuum. 

In the last two cases I have calculated the energy spectrum of the photons 
produced, while in the first one it resulted more natural to calculate the 
energy spectrum per unit time. In particular, only in the third case an an- 
alytical spectrum has been obtained. To do this, the effective Lagrangean 
technique has been used in such a way the resulting spectrum is reliable 
only for photon energies much less than the electron mass m. In this case 
I found that the spectrum becomes asymptotically independent of the pho- 
ton energy u> [see Eq. (|4.3.11|) ] while the experimental GRBs spectra go in 
the same limit as a; -1 . Nevertheless, I have shown explicitly that the direct 
photon production mechanism is not quantitatively very efficient and it gives 
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an energy spectrum completely negligible with respect to the corresponding 
synchrotron spectrum [see Eq. (|4.3.12j) ] . The comparison with the photon 
production process through pair annihilation is less evident. In fact, in this 
case I have calculated the photon energy spectrum per unit time result- 
ing from the annihilation of two identical "distributions" of electrons and 
positrons created from vacuum in the presence of the slowly-rotating mag- 
netic field (|4.1,1|) [see Eqs. (|4,1.3|) and (|4.1,2|) ], The "annihilation" spectrum 
so obtained (see Fig. I4.ip shows a sharp peak (break energy) around the 
electron mass and a rapid increasing (decreasing) below (above) m and it is 
very different from the experimental GRB energy spectra (see Fig. I1.4D . 

Instead, a much better qualitative agreement has been obtained between 
the synchrotron spectrum and the experimental ones (see Fig. I4.3H . In this 
case, I started from the Lagrangean density of QED in the presence of the 
external magnetic field B f (t) and I transformed it into a more useful "effec- 
tive" Lagrangean density in order to exploit the fact that B^(i) had been 
assumed, in fact, slowly-rotating. By calculating the S'-matrix elements cor- 
responding to the process of pair creation and subsequent photon emission 
by the electron or by the positron (see Fig. I4.2p I finally obtained the photon 
spectrum (I4.2.63P that is shown in Fig. 14.31 when B/ = 2.2 x 10 14 gauss. 
Similarly to the GRBs experimental spectra, the theoretical spectrum shows 
a double decreasing behaviour with two different slopes around the break 
energy. In this case the value of the break energy depends on the magnetic 
field strength and if Bf ~ 10 14 gauss then the break energy lies between 
1 MeV and 3 MeV as some experimental GRBs spectra. Most important, 
the low- energy region of the synchrotron spectrum shows a linear dependence 
on the inverse of the photon energy exactly as the energy spectra of GRBs. 
Finally, it is worth noting that, instead, in the high-energy region above the 
break energy both the annihilation and the synchrotron spectra decrease 
too rapidly with respect to the experimental GRB spectra. As I have said, 
this common feature is due to the fact that in my model the production of 
high-energy pairs (and then of high-energy photons) is not very efficient [see 
Eq. gJLH)]. 
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Pair production in a strong, 
uniform and slowly- varying 
magnetic field: the effect of a 
background gravitational 
field 

The aim of this Chapter is to continue the study of the production of 
electron-positron pairs in the presence of a strong, uniform and slowly- 
varying magnetic field started in Chap. 3. As I have pointed out in Sect. 
11.21 the- physical situation I have in mind is the creation of pairs around 
astrophysical compact objects like magnetars or black holes. But, till now I 
have performed all the calculations in a flat spacetime, that is neglecting the 
effects of the gravitational field created by the compact object. Even if there 
are situations in which this can be safely done [1001 1101] . it is interesting 
to study what happens if the effects of the gravitational field are taken into 
account. This is done in the present Chapter. In particular, after stating 
the general assumptions of the theoretical model (Sect. 15. ip . in Sect. I5.2l the 
structure of the gravitational field is assumed to be such that its effects can 
be calculated perturbatively [102J while, in Sect. 15.31 the effects of a strong 
gravitational field are taken into account nonperturbatively [103J. 

5.1 General assumptions 

In the previous Chapters I have performed all the calculations in Minkowski 
spacetime. By reminding the astrophysical scenario I described in Sect. 11.21 
this means that the effects of the gravitational field of the central engine 
producing a GRB have been neglected. As I have said at the beginning of 
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Par. 13.1.11 if one considers the production of pairs around magnetars, this 
assumptions is very realistic. Nevertheless, we have seen in Par. 11.2.31 that 
many models identify the central engine of GRBs with a massive rotating 
black hole surrounded by an accretion disk. In this case, obviously, the 
strength of the gravitational field can be such large that its effects could be 
relevant. As I will explain in the following Sections, these effects are taken 
into account by working in the framework of quantum field theory in curved 
spacetime (see Sect. I2.3p . Actually, in my model the main responsible of 
the pair creation is still the magnetic field and the fact that it varies with 
time. For this reason, the final presence probabilities will be calculated 
again by using the first-order adiabatic perturbation theory. Nevertheless, 
the presence of the black hole gravitational field will be taken into account 
in the determination of the one-particle electron and positron modes and 
energies and this fact will make the presence probabilities calculated here 
different from those obtained in Minkowski spacetime [102|, 1103] . 

According to what I have just said, I assume that the spatial structure 
and the time evolution of both the gravitational field and the magnetic field 
produced by the black hole are given. Nevertheless, in order to determine 
a realistic form of these fields I should fix the physical properties of the 
source (its mass, its eventual electric charge, its angular momentum and 
so on) and solve the system built up by the Einstein equations and the 
general covariant Maxwell equations. Clearly, this is a hopeless problem and 
a number of approximations have to be done. In particular, I first assume 
that the Einstein equations and the general covariant Maxwell equations are 
disentangled. This corresponds to neglect the gravitational field produced 
by the magnetic field and to assume that the spacetime metric is determined 
only by the black hole. In order to further simplify the problem I assume that 
the spacetime metric is actually that produced by a spherical, uncharged, 
nonrotating black hole and that the corrections to this metric due to its 
eventual charge, to its rotation and to its magnetic field can be neglected. 
Nevertheless, the black hole itself should be capable to produce a time- 
varying magnetic field. This request can be satisfied without changing what 
I have said about the spacetime metric, if the spherical body is collapsing 
but keeping its spherical symmetry and without rotating. In this case, in 
fact, the metric tensor outside the body does not change because of the 
Birkhoff theorem [82j while the magnetic field is found to grow with time to 
compensate for the decreasing of the gravitational energy of the collapsing 
body. 

In these approximations, as I have said, my starting point is the metric 
tensor corresponding to the field created by a spherical body with mass M 
outside the body itself. If I call t, and R = (X, Y, Z) the four coordinates, 
this metric tensor is a function only of the radius R = V X 2 + Y 2 + Z 2 and 
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it can be written as [8 

^ w=dhl! [(lz|)-,_ (l+ ^,_ (l + ^,_ (l+S y 

(5.1-1) 

with re = 2GM the gravitational radius of the body (G is the gravitational 
constant). I have chosen the so-called isotropic metric instead of the usual 
(and equivalent) Schwarzshild metric, because from Eq. (I5.1.ip one sees 
that the spatial distance is proportional to its Euclidean expression and this 
will simplify future calculations. I point out that in this metric the event 
horizon of the black hole is the spherical surface R = ro/4. 

Concerning the magnetic field, we have seen before that the rotating 
magnetic field configuration is very efficient from the point of view of pair 
production then it would be natural to consider here again the same con- 
figuration. Nevertheless, I have evaluated the corrections to the presence 
probabilities due to a weak gravitational field and, actually, there are not 
new interesting qualitative effects to be discussed [in any case they can be 
found in [102J]. Instead, my first task in this Chapter is to show that also the 
presence of a weak gravitational field superimposed to a time-varying mag- 
netic field with fixed direction makes different from zero the probability that 
starting from vacuum a pair is present with the electron and the positron in 
a TGS. For this reason I will consider in this Chapter the following magnetic 
field time evolution 

Bf p (f) = ( ) = ( ) (5.1.2) 

\B^(t)J \Bj + (Bi - B f ) exp(-i/r)/ 

with Bi < Bf and with r a typical macroscopic time characterizing the time 
evolution of the black hole (for example, its collapse duration). 

Since the spacetime is curved and the metric tensor is not simply r) a p, 
one must pay attention in defining the vector potential A^f(t, R) that gives 
rise to B| Xp (i). To this end, I assume that the threedimensional components 
of the magnetic field B| Xp (£) define the spatial-spatial components of the full 
covariant electromagnetic tensor F^J}(t, R) that then are uniform in space: 



F^(t) = -F^(t) = 0, (5.1.3a) 

F™$(t) = -F^(t) ee 0, (5.1.3b) 

F™?(t) = -F^(t) = B^(t). (5.1.3c) 

All the other full covariant components of F^(t, R) obviously vanish while 

the mixed or the full contravariant components are built by means of the 
metric tensor ()5.1.1|) . Now, also in curved spacetime 

F^l{t) = R ) " 9uA^>(t, R), (5.1.4) 
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then by means of Eqs. (I5.1,3p . I can choose the covariant vector A^(t, R) 
as 

^o P (t,R)=0, (5.1.5a) 
A^(t,R) = 1 -YB? p (t), (5.1.5b) 

A^(t,R) = ~XB^(t), (5.1.5c) 
^3 P (t,R) = 0. (5.1.5d) 

These equations define a gauge corresponding to the symmetric gauge in 
Minkowski spacetime [see Eq. (|2.1.4p ]. Finally, it is convenient to define 
also the threedimensional vector potential A| Xp (i,R) as 

A° t xp (t, R) = ~[Rx Bf p (t)] (5.1.6) 

where the minus sign has been inserted to have d x A° xp (t,R) = B° xp (t). 

Now, as I have done previously, I will calculate the probability that a pair 
is present at time t by applying the first-order adiabatic perturbation theory. 
To do this, I have to build up the second-quantized Hamiltonian of a Dirac 
field ^'(t, R) in the presence of the slowly- varying magnetic field (|5.1.2p 
and in the curved spacetime with the static metric tensor (|5.1.1j) and to 
determine its instantaneous eigenstates and eigenenergiesQ The Lagrangian 
density of this system is given by [see Eq. (|2.3.i 



(5.1.7) 



where all the quantities representing the gravitational field are defined in 
Par. 12.3.21 Now, in general, a pair is created in a spatial volume with typical 
length given by the Compton length A = 3.9 x 10~ n cm while for a typical 
10 solar masses black hole re = 3.0 x 10 6 cm. In this way, since A <C tq, I 
am allowed to make some simplifications on the metric tensor (|5.1.ip . In the 
following two Sections I will approximate the metric tensor (]5.1.ip and all 
the connected quantities appearing in Eq. (|5.1.7p in the two cases in which 
the pair is present far from (Sect. I5.2p or near (Sect. I5.3P the event horizon 
of the black hole. While in the first case the gravitational effects will be 
accounted for perturbatively, in the second one a nonperturbative approach 
will be followed. 



In this Chapter, I use the primed notation to indicate the Dirac field and the related 
quantities because I deal with a magnetic field varying only in strength and always lying 
in the Z direction (see the notation used in Sect. 12. ip . 
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5.2 Weak-gravitational field case 



In this Section I assume that the pair is produced in a small neighborhood of 
the space point labeled by the coordinates (X c , 0, 0) with X c > rc/4+ A and 
A > 0. If {X c + x, y, z) is a generic point in this neighborhood then |x| < A, 
\y\ < A and \z\ < A and I can calculate the metric tensor g^ u (X c + x,y,z) 
by keeping only the terms up to first order in x/X c , y/X c and z/X c . It can 
easily be seen that the resulting metric tensor depends only on x and that 
it can be written as 



9 { uXx)=9^+K u {x) 



(0) 



where 



JO) 
h^u(x) 



diag(g t ,-g s ,-g s ,-g s ), 
diag(2/ 4 x, 2f s x, 2f s x, 2f s x) 



with [see Eq. ([5X11 ] 

l-r G /4X c 



9t 



ft 



1 + r G /4X c 
l-rg/AXc r G 
{l + r G /AX c f2Xr 



fs 



1 + 



1 + 



4X C 

r_G 
4X r 



r G 

2XT 



(5.2.1) 

(5.2.2a) 
(5.2.2b) 

(5.2.3a) 
(5.2.3b) 



It is evident that, in order that g^J(x) is a good approximation of 
g^uiX-c + x, y, z), X c can not be chosen to be too close to the critical value 
rc/4. Just to give an idea, it easy to see that, if N S> 1 is a large pure 
number, then 



9w(X c + x, y, z) - g$ (x) 



with fj, = 0, . . . , 3, only if 



g^(X c + x,y,z) 



x c > r -j + Vnx. 



< 



N 



This condition automatically implies that 



h^{x) 



(0) 

g^i 



< 



l 



2VN 



(5.2.4) 



(5.2.5) 



(5.2.6) 



with \i = 0, . . . , 3 and then that h^ u (x) can be considered as a small correc- 
tion of g^v • I n what follows, I assume that the previous inequalities hold 
with sufficiently large N. I want to observe here that even for very large 
values of N, Eq. (|5.2.5|) does not constraint very much the values of X c : in 
this respect the expression "far from the black hole event horizon" is to be 
interpreted as "microscopically far from the black hole event horizon" . 
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Since g$(x) has been split as in Eq. (15.2.11) with the matrix h liu {x) 
much smaller than the matrix gffj , I am allowed to keep in the Lagrangean 
density (|5.1.7f) only the first-order terms in h fJiU (x). To do this, I observe 
that 

g {1) (x) ~g {0) [l + h{x)] (5.2.7) 

where 

g (0) =det(g^) = -g t gl (5.2.8a) 

h(x) = hHx) = 2 ( ^ - ^] x. (5.2.8b) 
M \9t g s J 

Also, being the metric tensor g^J (x) diagonal, I can choose a diagonal tetrad 
[see Eq. (j2.3.2()|) ] with 

4' )0 (x) = -=L= . -L f 1 - ^) , (5.2.9a) 

ef )i (x)= ; 1 ^-U(l + — ) (no sum). (5.2.9b) 

By means of this tetrad it can be shown that the connections hidden in the 
covariant derivatives ^ and S$ ^ are already first-order quantities given by 

r «(s) = ^eJO^ (5.2.IO) 
where is the diagonal zero-order tetrad with 



,(0)0 



(o)i _ - _ - ( no sum )_ (5.2.11b) 



1 1 



(5.2.11a) 



1 



9 



(0) y/9. 



Now, the procedure to calculate the Lagrangian density (|5.1.7p up to 
first order in h fJiU (x) is identical to that used in the weak gravitational field 
approximation |104[ 1105] and I give only its final expression: 



JSf'«(t) 



9 * 'I - / E x)[? 7 °(i^f') - (ido&ft *] 



2 

9s^fg~t 

2 



(1 - f P x) V'Ylidi + eA^(t, r)]V - tt'^ - eA e T f (t, r)] 7 ! *' 



V#sf (1 - f M x)m^ 

(5.2.12) 
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where r = (x, y, z) and where I defined the three coupling^! 



fE 



/m 



3/s 


3r G 


9s 


2X1 


Vs_ 


_ft = 


9s 


9t 




_ft = 


9s 


9t 



1 + 



4X 



IX}. 



1 



4X 



4X 



-i 



1 + 



1 + 



4X 



4X r 



(5.2.13a) 
(5.2.13b) 

(5.2.13c) 



Note that the modifications induced by the metric tensor (|5.2.ip in the 
Lagrangian density (|5.2.12f) are linear in ft and f s (obviously) but nonlinear 
in g t and g s . 

According to what I have said in Sect. 12.31 I define the Hamiltonian 
density of the Dirac field ^?'(t, r) as 



jfW(t) = fY«(do#') + (Sb^')n /(1) -^ ,{1) (t) 



where 



n'W(t,r) = 



"3^(1 - /EX) (J , 

?l n 7 * (* 5 r) 



(5.2.14) 

(5.2.15a) 
(5.2.15b) 



are the first-order conjugated momenta to the fields ^'(t, r) and ^'(t, r) re- 
spectively. By using the previous equations it can easily be shown that, apart 
from derivatives terms, the Hamiltonian density (|5.2.14p can be written as 

jif'W (t) = - f E x)^H' {1) (i)*' (5.2.16) 

where I introduced the one-particle first-order Hamiltonian 

W /(1) (i) = ^g t {-±-\(l- f P x)a ■ [-id + eAf p (t,r)] 

+ a ■ [-id + eA! x p (t,r)](l - f P xj\ + (1 - /mx)/?™ 



-a ■ [-id + eA* xp (t,r)]+ f3m 



(5.2.17) 



Despite its appearance, the previous one particle Hamiltonian (|5.2.17f) is 
an Hermitian operator. In fact, the Hermiticity depends on the definition 



2 I introduced three couplings for later convenience because, actually, only two of them 
are independent. 
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of the scalar product and I remind that the scalar product of two spinors 
^i(t, r) and r ) is defined in curved spacetime as in Eq. (|2.3.37p . By 
choosing the Cauchy surface S as the t = const, hypersurface, then dT, = dr, 
n M (i, r) = (1, 0, 0, 0) [87] and the scalar product (|2.3.37j) becomes up to first- 
order terms 

tf 2 ) (1) = / dr^l - f E x)v[* 2 (5-2.18) 

and the one-particle Hamiltonian (15.2. 17j) results, in fact, to be Hermitian. 
The fact that this one-particle Hamiltonian is a "good" Hamiltonian is also 
corroborated by the fact that the equation of motion of the field *$>'(t,r) 
given, in general, by 

can be written in the present case as 

idoV =H'W (i)W'. (5.2.20) 
Coherently, the total Hamiltonian of the Dirac field is given by 

H'M(t) = J drJT'W(t) = J dry/fi{l - f E x)V*H' m (t)V . (5.2.21) 

Now, this Hamiltonian depends explicitly on time only through the depen- 
dence of Ti'^'it) on the magnetic field B° xp (t) hidden in the vector potential 
A^ xp (t,r) [see Eq. (I5.1.6P ]. then it is a slowly- varying quantity. In this way, 
also in this case, in order to calculate the pair presence probability I can 
use the first-order adiabatic perturbation theory. Since the conceptual steps 
are the same as those followed in Chap. 3, I will not repeat all the details. 
In particular, I have to determine the one-particle electron and positron 
modes and energies of the time-independent Hamiltonian obtained from Eq. 
(|5.2.17|) by substituting the vector potential A^ xp (t,r) with A'(r) given in 
Eq. (|2.1.4I) and corresponding to the static magnetic field B'(0, 0,B) [see 
Eq. (|2.1.ip ]. This one-particle Hamiltonian can be written as the sum 

=n >(o) + z> (5.2.22) 

of the zero-order Hamiltonian 

w /(0) = 191 [a ■ [-id + eA'(r)] + ^g~ s [3m} (5.2.23) 
V 9s 

and of the first-order interaction 

T' = vWp " fM)(3mx - ^{x,H'^} + f E xH'®. (5.2.24) 
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In this way, in the present physical situation the one-particle electron and 
positron modes and energies can be determined perturbatively in the cou- 
plings /e, fp and /m (or, equivalently, ft and f s ) by using the time- 
independent perturbation theory [106J. In particular, in the following Para- 
graph I will determine the zero-order electron and positron modes and the 
corresponding first-order energies. Instead, in Par. 15. 2. 21 1 will calculate the 
first-order corrections only to the zero-order TGSs previously obtained. 



5.2.1 Computation of the one-particle modes up to zero or- 
der and of the one-particle energies up to first order 

If Uj(r) and Vj(r) with j embodying all the needed quantum numbers are 
the electron and positron modes then they satisfy the equations 

H' {1) U' 3 =w 3 U' r (5.2.25a) 
H'MVj = -w 3 Vj (5.2.25b) 

where w 3 and w 3 are the electron and positron one-particle energies. The 
states kj(r) and Vj(r) are assumed to be an orthonormal basis with respect 
to the scalar product (|5.2.18p . i.e. 

(U>,U> f )U = (y;,y/,) (1) = S Jtf , (5.2.26a) 
(L^,y;,) (1) =0. (5.2.26b) 

I first want to determine the zero-order solution of Eqs. (|5.2.25p . Up to 
this order those equations can be written as 

- {a ■ [-id + eA'(r)] + PJgSm} = w^U'^, (5.2.27a) 
9s 

{cx ■ [-id + eA'(r)] + pjgSm] V'^ = -w^V' {0) . (5.2.27b) 
9s ' 

These equations are the eigenvalue equations in Minkowski spacetime of a 
particle with mass ^/g^m and charge — e < in the presence of the magnetic 
field B', then their solutions can be written immediately. In particular, 
j = j = {rid, k, cr, n g } and wf 1 ^ and wf 1 ^ are given by the modified Landau 
levels 



wf = Jgtm 2 + — [k 2 + eB{2n d + 1 + a)], (5.2.28a) 
1 9s 



wf ] = Jg t m 2 + — [k 2 + eB{2n g + 1 - cr)]. (5.2.28b) 



9 



Instead, the corresponding modes will be indicated as u'j (r) and v'j (r) 
respectively and they are given by Eqs. (|2.1.23p with Jg~ s m instead of m, 
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instead of Wj and with instead of Wj. Now, analogously to what 
I have said in Par. 12.1.21 the one-particle electron and positron energies 
(I5.2.28P have two kinds of degenerations. The first one is due to the fact 
that they do not depend on the quantum number n g {rid). The second one 
is due to the fact that the electron (positron) modes with quantum numbers 
3+ = { n d, k, +1, n g } and j_ = {n d + 1, k, -1, n g } (j + = {n d , k, +1, n g + 1} 
and j- = {rid, k, — 1, n g }) have the same energy whatever n g (rid). This 
means, following the time-independent perturbation theory, that the modes 
u'j°\r) and v'j°'(r) will not represent, in general, the correct zero-order 
modes U'^°\r) and Vj (r) and, for this reason, they have been indicated 
by means of the symbols u' (r) and v'j°\r). 

Now, in the following, I will compute explicitly only the zero-order elec- 
tron modes and the first-order electron energies, while the analogous results 
for the positron modes and energies will be only quoted. Following the time- 
independent perturbation theory for degenerate states, I write the zero-order 
solutions of Eq. (|5.2.25ap with a given energy as linear combinations of all 
the degenerate modes u' (r) and v'j°\r) with that energy. Now, it can 
easily be shown that the perturbation I' does not remove the energy degen- 
eracy of the modes characterized by the quantum numbers j- and j+ [see 
|lU2j for further details]. In other words, following the time-independent 
perturbation theory for degenerate states |106j . I have to diagonalize the 
perturbation X 1 inside every subspace labeled by the quantum numbers j 
whatever n g . In this way, the "true" zero-order modes are characterized by 
the quantum numbers {rid, k, cr} and by a new index that, for later conve- 
nience, will be indicated as xn. If I call Jj' \ \ „_ (r) the resulting zero-order 
modes then 



^n d ,k,a,xo^^ ^ ~] P ' n d ,k,cr,xo;n d ,k,a,n g '^'n c i,k,a-,ng^^ (5.2.29) 
n g =0 

where the coefficients , „ .„ u „ solve the secular equation I1U6I 

ri dl K,u,XQ,n d ,rL,(j,ri g B s 

oo 

XT \^n d ,k,a,n g ;n d ,k,<r,n' g ~ € n d ,k,a,x $n g ,n g J P i d \k,a,x ;n d ,k,a,n' g = (5-2.30) 
n' g =0 

with, in general, 

4, = J dv^ s ^{v)l'u' 3l {v) (5.2.31) 

and with e nd kaxo the first-order corrections (to be determined) to the en- 
ergies of the modes ^^^(r). 

Now, one can show that the matrix elements X' h „„ .„ u„„, can be 
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written as 



1' 



n d ,k,a,n g ;n d ,k,a,n' g 



it f \ 9tm* . (o) 

Up ~ fM)—($ + Ue ~ fp)w) 



w 



j dxdy6'* d ng (x,y)x 8' nd n , g (x,y) 

(5.2.32) 



where the operator xq and the functions 0' n (x, y) have been defined in 
Eqs. (|2.1.12ap and (|2. 1.281) respectively. Prom Eqs. (12.1.230 and (I2.1.26|) 
we see that only the transverse functions 9' ndng (x,y) in u'j (r) depend on 
n a , then I have to determine the coefficients Pi°^ , „ „ ._ . „ „ in such a way 

y" n d ,K : cr : XQ\n d ^K,a i ng J 

that the linear combination _ n Pi°\ „ .„ . „9'„ (x, y) diagonalizes 

^— m ff =U n d ^k^G^x^\n d ^k^u^rig n d ,n g \ o 

the operator rcoH This linear combination is given in [73] [it is Eq. (104) in 
Complement Eyi]. The coefficients ^n d k a x n d k a n g resm t independent of 
the quantum numbers {n^, k, a} and they are given by 



,(o) 

n d ,k,a,xo;n d ,k,cT,ng 



PffiO*)) 



2 

(5.2.33) 

where I pointed out that xq is a continuous quantum number and where 



2U g roo 

R nq (z) = ^ ds(z + is) n * exp (-s 2 ) (5.2.34) 

V 7 !" J-oc 



is the ra 9 th-order Hermite polynomial [107] . As a conclusion, the spinors 

U'nk o-( x 0' r ) have the same form of the spinors (r), but with the func- 
tion &' ndjn (%,y) substituted by the function 



e'(x ;x,y) 



oo 

E 

n o =0 



P^oK d ,„>,y) 



(5.2.35) 



3 Now, it is clear why I have called the additional index xo- 
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By using the second expression of 9' nd rig 

(4)^ 4 [e~B fe~B 
T 



, y) in Eq. (12.1.281) I obtain 



nd 

oo 



TT 




exp 



ds exp 



eB_ 
2 

1 



2 x 2 + y 2 
+ o 



1 

7T 



— - (V eBxo + is)VeB(x — iy) 

n. _ 





exp 



=o 9 
eB 



2 X 2 + ?/ 2 
^0 + n 



1 

TT 



x exp [eBxo(x — iy)] 



"a 





exp 



ds exp 



-s + iv eB(x — iy)s 
a iy{x - 2x )] 



(5.2.36) 



From now on, the quantum number n g completely disappears in the 
calculations because it is essentially substituted by the quantum number 
xo- This circumstance allows me to simplify the notation. In fact, I can 
eliminate the index "d" from the remaining quantum number and from 
the related operators ad and a\ [see Eqs. (I2.1.20ap and (|2.1.20bj) ]. Also, 
later calculations will be simplified if I discretize the eigenvalues xq. This 
can be done by requiring the functions @' nd (xo;x, y) to satisfy the periodicity 
condition at x = 



exp 



eBxoL l 



exp 



.eBx$L t 



(5.2.37) 



where L y is the length of the quantization volume in the y direction. In this 
way, the allowed eigenvalues are given by the discrete values 



xo 



2£tt 
eBLi, 



0,±1,. 



(5.2.38) 



I point out that if I imposed the periodicity condition at x' ^ 0, the allowed 
eigenvalues would change. Nevertheless, since at the end of the calculations 
I will perform the continuum limit L y — > oo, I am not interested in the 
exact values of the allowed eigenvalues but only in the eigenstate density 
q(xq) which is 

independently of xq. In this way, the function (|5.2.36|) becomes 



®n,x (x,y) 




exp 



eB 



[(x - x ) 2 - iy{x - 2x )] 



(5.2.40) 
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The numerical factors have been chosen in such a way that if I define [anal- 
ogously to the twodimensional spinors tp'Ar) given in Eq. (|2.1.26ap ] the 
twodimensional spinors 



eMlkZ he' n ^,y) 



'L, 



with 



J = {n, k, a, xq}, 
they result normalized as [see also Eq. (|2.1.31|) ] 



(5.2.41) 
(5.2.42) 



J dr<!>'}(r)& JI (r)=6j,j, 



(5.2.43) 



with Sj t j> = S Tljn i6i ll} ktS <T)(T iS XOta / . Finally, the zero-order spinors U'j°\r) are 
given by [see Eq. (|2.1.23a|) ] 



U'j°\r) 



(°) , i— 



2w 



(0) 




v 



a (0) , ~J 

V s Wj + ^fgtra 



(5.2.44) 



and they are normalized as [see Eq. (|5.2.18jl ] 



dr^ s uf )] {r)uf\r) 



ij,J' 



(5.2.45) 



In Eq. (|5,2.44p I used the fact that the energies (|5,2.28ap do not depend 



neither on n„ nor on xq then w) 



(0) = (0) _ (0) (0) 

J j — W n d ,k,a,n g ~ W n,k,a,x ~ W J • 

Instead, from Eq. (j5.2.32|) the energy corresponding to the mode U'j ^ (r) 



is, up to first order, 



(i) 



(0) , 

W J + 



tf t \9tm f n (o) 

Up ~ fM)—pf + Ue - fp)w) 



w 



Xq 



g t m 2 + — [k 2 + eB(2n + 1 + a) 
9s 



(5.2.46) 



x 1 + 



ft { fs k 2 + eB(2n + 1 + a) 
g t g s g s m 2 + k 2 + eB(2n + 1 + a) 



X Q 



where I used the definitions (15.2. 13)) of the coefficients fp and fu- In 
an analogous way one can write the positron energies up to first order as 



w 



(i) 



g t m 2 + — [k 2 + eB(2n + 1 - a)] 
9s 



x 1 + 



ft { fs k 2 + eB{2n + l-a) 
g t g s g s m 2 + k 2 + eB(2n + 1 - a) 



(5.2.47) 
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I observe that, as in the case of the electron modes and energies, the quantum 
number is here substituted by Xq, and it is pointless to keep the index u g" 
in the remaining quantum number n g because there is no more possibility 
of ambiguity. 

Finally, the zero-order positron modes are given by [see Eq. (12. 1.23b j) ] 



vr(r) 



~c°) , i— 



V 



2w 



(0) 



9s yjj 1 ^ + ^fgim 



(5.2.48) 



with [see Eq. d2.1.26b|> ] 



X j( r ) = 7=f= f-o®n,x ( x ,y) 



and they satisfy the orthonormalization relations 

j dr^ s V'} 0) \v)vf\v) = 5j,j, 
j dv^ s Vj m {v)uf\r) = 0. 



(5.2.49) 

(5.2.50a) 
(5.2.50b) 



5.2.2 Computation of the one-particle TGSs up to first order 

I should pass now to the determination of the one-particle states U'Ar) and 
V1(y) up to first order in the perturbation X'. But, from Eqs. (]5.2.46|) and 
(|5.2.47p we see that the first-order energies of the TGSs, that are given by 



'k,xo 



^0,fe,-l,a;o 



W, 



r,(l) 

0,fc,+l,a;o 



gtm 2 + — k 2 



9s 



1+ lft + fs 



k 



9t 9s g s m 



2 + k 2 



X 



(5.2.51) 



are also independent of the magnetic field strength as in Minkowski space- 
time. As we have seen before, this fact gives the TGSs a particular rele- 
vance in the presence of strong magnetic fields because their energies are 
much smaller than the excited Landau levels. For this reason, in the follow- 
ing I will only calculate the presence probability of a pair with the electron 
and the positron both in a TGS and then I need to compute only the one- 
particle TGSs corrected up to first order. As in the previous Paragraph, 
I will present the calculations of the first-order corrections to the electron 
TGSs and I will quote the analogous results for the positron TGSs. 

The first-order corrections to a given zero-order electron TGS come from 
the coupling of this state with the following classes of states: 



1. the zero-order electron TGSs with the same energy (all but the state 
to be corrected); 
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2. the zero-order electron TGSs with different energy; 

3. the zero-order electron states that are not TGSs; 

4. the zero-order positron states. 

Now, I suppose I want to calculate the first-order corrections to the TGS 
labeled by Jo = {0, k, — 1, xq}. The states in the first class are labeled by the 
quantum numbers {0, k, —1, x' } with x' ^ xq because they have the same 
energy of the state labeled by Jo- But, all the contributions vanish because, 
in general, the perturbation (I5.2,24p can not couple two modes of TC'^ with 
the same and two different xo and xL. Similarly, the states in the second 
class are characterized by the quantum numbers {0, k' , — l,x' } with k ^ k 
and, since [I',V Z ] = they do not give any contribution. Instead, the 
contributions from the states of the remaining two classes are, in general, 
different from zero in such a way the state U'}^ (r) can be written as 

= uT(r) + ?%,UT (r) + £ Q^v}°\r) (5.2.52) 
J' J' 
where the primed sum does not include the TGSs and where [106J 



>(1) 
J ,J> 



1 



.(0) 



w 



(0) 



uT\r)T'uTJ(r 



(o). 



Jo 



Q 



(i) 

Jo,J' 



-(0) 



+ w 



(0) 



J dv^ s V'P\v)X'u'iyv). 



(5.2.53a) 
(5.2.53b) 



In these equations I introduced the zero-order energies of the TGSs e k 
defined as [see Eqs. (15.2. 28|) ] 



(0) 



-(0) 



w, 



(0) 

0,fe, — l.xo 



w, 



(0) 

0,fc,+l,a;o 



gtm 2 + —k 2 . 



9s 



(5.2.54) 



Now, I start by calculating the coefficients ' Jt . From the expression 
(15.2.240 of the interaction Hamiltonian and from Eq. ()5.2.44p we have 



Jo,J' 



1 



-(0) 



(0) \ 

W J> \ 



w 



(0) 



9tm 



2 W f 



(0) , i— 



(0) 



fi 



2e 



(o) , _(o) 
j' i_e fc 



9t(fp - fhl)m 
j dr$'J,(r)x$j (r 

-VqIUp ~ fu)m - l£(w® + ef) + f E el 
1 1 



fp, 



(0) 



9s 



JO) 



(0) 



j dr&},(r)V'xV& Jo (r) 

(5.2.55) 
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By using the ortho-normal properties (|5.2.43|) of the functions $j(r), I can 
write the quantities Pj ji as 



P Jo,J> = { B k,'x 5 n', +1^,-1 ~ iC k,'x S n',oSa',+l ) h',A' , Xo (5.2.56) 
where I defined the coefficients 



■(i) 



B 



« - 1 

k,x — (0) (0) 



4 0) + V9t m 



(0) , /— 



-(0) 



x < V^(/m - fp) 



rn 



V2eB 



1 



ft 



9s (4 0) + v^™)(4 UJ + Vft™)J 



(0) 



\/2ZB 



1 + 



fc 2 



ft 



ft (4°)+^ m )( £ w+^ m ) 



(0) 



(5.2.57a) 



C 



1 



(i) 

Mo £ m_ £ (0)\\ 2g (0) \\ 2£ 



£ k + \/ft m 



(0) , /— 

4 + Vft"i 



^(/M-/p)m-^(4 o) +4 0) )+^--;; 



(0) 



(5.2.57b) 



ft 



fc 



5s (4 0) + Vft^)(4 0) + Vft^) 

with [see Eqs. (I5.2.28|) ] 



,(°) 



c(o) _ ,(o) 

£ fc — U, Q,Jfc,+l,xo ~~ U 'l,Jfc,-l,xo ~~ ^O^j-ljio ~~ W l,k,+l,x 



<7 t m 2 + ^ (fc 2 + 2eB) 
ft 



(5.2.58) 



the zero-order energy of the first-excited Landau level. In the same way, I 
can write the coefficients Qj^j' as 



(1) 



with 



D 



(5.2.59) 
(5.2.60a) 



« 1 

ft, 2:0 / ^q^\ 2 



2 e 



ft(/M - fp)m. —kx , 
V ft- 
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r(l) 

k,xo 



(0) , /— 

£ k + \/5t m 



-(0) 



9s 4 0) + V&m 



(5.2.60b) 



c(°) i r- 



(o) . _ , — , 

£ k + \/^ m /s^ fc 



x < Vat(fM - fp)m 



-(0) 



+ 



5s v^e5 
1 \ 



(0) , i— c(0) , /— 



+ 



f(4 0, -4VM > 



(o) , ^- c(0) i ^~ 



(5.2.60c) 



If I also define the coefficients fi^ as 



k,xo 



fE 



-Xq 



(5.2.61) 



then, the first-order TGS U'J^\r) = U' [1 k '_ l x (r) can be written simply as 



'(i) 



C/' (1) (r) 
u 0,k,-l,x \ r ) 



1 + ^i,x ) ^o!fc,-l,a!o( r ) + ^k,x ^l,k-l,x ( r ) ^k,x ^0,k',+l,x ( r ) 



>(i) rr'(o) 



(i) rr'(°) 



U k,x V 0,-k,+l,x \ l ) ^ tI1 k,x v O,-k,-l,xo\ 1 ) 



(1) t/(0) 



p(l) t/(0) 
r k,x v l 



k,+l,x \ 



(5.2.62) 



(for the sake of clarity the explicit expression of the zero-order states appear- 
ing in the right hand side of the previous equation are given in Appendix 
D). The term A^t/^^r) = /^ <°i_ Mo (r)/2 ha s been added to 
compensate for the factor (1 — Jexq) in the scalar product (15.2. 18)) and then 
to have the states correctly normalized up to first order, as 



(u'j;\u'}P)V = 8, 



Jq,Jo 



(5.2.63) 



where Jo = {0, k, —1, xq} and J' = {0, k', —1, x' }. 

Finally, with analogous calculations it can be shown that the first-order 
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positron TGS Vq^ +1 (t) can be written as 



V' {1) (r) 



1 + A (1) "| V' {0) M + B (1) V' {0) (r) - iC W V' {0) (r) 
K x ^ ^k,xoJ v 0,k,+l,x Q \ 1 J t D k,x v l,k,+l,x Q \ L ) LK -k,xo v O,k,-l,xo\ L J 

- D {1) U' {0) (r) + iE {1) U m (r) - F (1) u' {0) (r) 

u k,x U 0,-k-l,x \ 1 y ^ l ' z k,xo U 0-k,+l,xo\ L > r k,x u l,-k ,-l,x V 1 /• 

(5.2.64) 

and that the satisfy the orthonormalization relations 

(V£\V^)M=6 J0>4 (5.2.65a) 
(Vj^\u'}P)W = 0. (5.2.65b) 



5.2.3 Calculation of the presence probability 

As I have said in the previous Section, I can calculate the pair presence 
probability in the presence of the slowly- varying magnetic field (|5.1.2p and 
of the static gravitational field described by the metric tensor ([5.2, lft by 
means of the adiabatic perturbation theory up to first order in the time 
derivative of the magnetic field. In fact, from now on, the gravitational field 
will not play any further role: I took into account its presence by correcting 
the one-particle electron and positron modes and energies. Nevertheless, 
in order to avoid the possible confusion between the "first-order" relative 
to the adiabatic perturbation theory and the "first-order" relative to the 
gravitational couplings fp and /m, I stress that in what follows I will 
always refer to the second one. In particular, the superscript (1) indicates 
quantities that are first-order in the gravitational couplings. 

Now, I pointed out in Par. 13.1.21 that in Minkowski spacetime if the 
magnetic field changes with time only in strength, the probability that a 
pair is present with both the electron and positron in a TGS is zero. In the 
following we will see that the first-order corrections to the TGSs I calculated 
before will make this probability different from zero in the present physical 
situation. In fact, the gravitational field lies in the x direction [see Eq. 
([5.2. ip ] in such a way the first-order TGSs are not eigenstates of a z as in 
Minkowski spacetime and then the selection rule ([3.1.23P does not hold. 

If the electron and the positron are assumed to be present at time t 
in the TGS state uf\r,t) with J = {0, k, -1, x } and Vj { ^(r,t) with 

Jq = {0, k' , +1, x'q} respectively, then the presence matrix element is given 
by 

^V* ) = Vf^T^ / ^-/^(r.tjtrxaj^^r.t). 

(5.2.66) 
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The presence amplitude at time t can be calculated from this matrix element 
as 



7 (1) - (t) = - 

k,x 1 k',x' n 



1 



dt'H% (0 exp 



l £ 



.(1) 



(5.2.67) 



where I used the fact that the first-order energies £^\ of the TGSs do not 
depend on B^ p (t) and then on time [see Eq. (j5.2.51j) ]. 

Now, the selection rule (|3.1.23p concerning the spin of two TGSs allows 
me to conclude that for the zero-order TGSs the following equalities hold: 



J drU'j^{r,t) (r x a) z U'}?\t,r) = 0, (5.2.68a) 

J drU'j^(t,r)(rx a ) z V^ ) (t,r)=0, (5.2.68b) 

f drV'- m (t, r) (r x a) _ vf ) (t, r) = 0. (5.2.68c) 

By exploiting these equations and by keeping only the terms up to first 
order, I can write Eq. (|5.2.66p as 



+ 



^ei? t exp (t) 



dry/ 9% 



+ ^l^Sli,^ r) - i^ (t)V<%_ hxo {t, r) 

(5.2.69) 

This expression can be further simplified if I use the definitions (I2.1,21ap 
and (|2.1.21cp of the operators x and y and the fact that the operator xq is 
diagonal with respect to the basis {U'j (t, r),Vj (t, r)}. In fact, in this 
way the operator (r x a) becomes 



r x a) 



IOL. 



xo + iyo + \j —gd ] 



ia + | x - iy + \ j —a 



(5.2.70) 
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with a± = (a x ± ia y )/2. Now, from the expressions (jD.laft and (ID. lb}) 
of the zero-order TGSs I conclude that only the operator ia_(xo + iyo) 
gives a nonvanishing contribution in the first integral in Eq. (|5.2.69p and 
that, analogously, only the operator — ia + (xo — iyo) gives a nonvanishing 
contribution in the second one. In conclusion, the matrix element (I5.2.69P 
can be written as [see Eq. (I2.1.15P ] 



l9t -eB?»(t) 



xo + 



1 



eB^{t) 
1 







<l^(t)|a_|V^(f)> 



Xo 



eB° xp (t) 



8 



(U'j^t)\a + \Vp\t)) 



(5.2.71) 



where 



(U'}?(t)\a-\vp>(t)) 



'(i), 



/ 



aT V 9s u o,k-x,x 



(t, r)a_ 

x >^(0^, + iy («»r)-<iL & (*)C.-i^( t ' r ) 

& (t,r)-F^(0E^ ) i | _i^( < » r ) 



0, /c' j-l-l^d 



(5.2.72a) 



<C/^(t)|a + |^(t)) 



'(i), 



dry 7 #f 



-iE {1) (t)V' m (t r) - F (1) (t)V' m (t r) 

tIz k,xo \ l ) V 0,-k,-l,x V 1 ' 1 ) r k,x W v l,-k,+l,x v 1 ' r J 
,'(0) 



xa +V,+i^(*> r )- 



(5.2.72b) 



The calculation of the matrix elements (|5.2.72p is quite tedious but straight- 
forward. In particular, it can be shown that 



(v2Ht)\a-\V$>(t)) = -(uX>(t)\a + \V^>(t)) 



'(i), 



1 



1 



9t(fM - fp)mk5k-k'5 x „ x '. 



(5.2.73) 

For this reason, the terms in Eq. (|5,2.7ip with the derivative with respect 

f(i) ( 

^Bf P W 1 



to xq cancel each other and the final expression of H^}\, (t) is 



^ t cxp w 2 4 0) 



St (/at - fp)mkxo5 k -k>5 Xo y o . (5.2.74) 
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The creation amplitudes at time t can be calculated by means of Eq. 
()5.2,67p and the only one different from zero is equal to 

%,k,-l,xo;0,-k,+l,xo&) 

SO) 



>9t k . , f t/T , , ex P 



gtifu - fp)mx / ds 



a 'A(e®Y Jo exp(-s') - Bf/(B f - Bi) 

(5.2.75) 



where only the first-order terms have been kept and where the time deriva- 
tive of .B| Xp (i) has been substituted. Now, as I have said below Eq. (|5.1.2|) . 
r is a macroscopic time connected with the typical evolution times of a black 
hole, then I can safely assume that ej?"' T This allows me to give an 

asymptotic estimate of the remaining integral for large times t. The result 
is 

(l) fgt k mx B f - Bi 



4 

(5.2.76) 

Finally, by squaring the modulus of this expression and by multiplying the 
result by the number of states for large times [see Eq. (|5.2.39|) ] 

eB ° XP{t)Lv dx ^dk = -^f-dV^dk (5.2.77) 
2-7T 2tt 4^2 g3 

with dV(°' = ^/g^LyLzdxo the "physical" quantization volume up to zero 
order [see Eq. (I5.2.2ap ]. I obtain the differential probability that a pair 
is present with the electron (positron) between xq and xq + dxo and with 
longitudinal momentum between k and k + dk {—k and —k — dk) as 

(5.2.78) 

In this equation the continuum limits L y — > oo and L z — > oo are understood 
and, since the probability is already proportional to f s and our calculations 
are exact up to first order in f s and ft, it is enough to use the zero-order 
"physical" volume dV^°\ Finally, the corresponding probability per unit 
volume and unit longitudinal momentum is given by 

dPM(x ,k;t) ( f s \ 2 y/gS eB f k 2 f Bf - BA 2 2 



dV(°)dk \l§irg s ) 9t (g s m 2 + k 2 ) 3 \ tB, 

(5.2.79) 

It is worth pointing out that the presence of the square of the length xq 
in Eq. (|5.2.79|) is the counterpart of the presence of R±m ° r R± m in the 
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presence probabilities calculated in Minkowski spacetime [see Eqs. (I3.2.2p . 
(|3.2.3ip and (|3.2.45|) ]. Also, note that, as it must be, due to the presence 
of gt in the denominator, the probability (|5.2.79p grows as the expansion 
center X c is moved towards the event horizon rc/4 [see Eqs. ()5.2.3bj) ]. 

Finally, as I reminded in Sect. 13. 2\ in [SS] it was calculated the total 
presence probability of a pair in the presence of a slowly-varying magnetic 



field with fixed direction but in Minkowski spacetime [see Eq. (|3.2.2|) ]n In 

order to have a quantity to be compared to the total probability per unit 
volume (|3.2.2|) . I have to integrate Eq. (|5.2.78|) with respect to k. After 
this integration and indicating the resulting probability per unit volume as 



where I assumed, for simplicity, that B{ ~ Bf, that R±m = and that 
the magnetic field (|3.2.ip is such that Bq = Bi and b = (Bf — Bi)/t. Now, 
as I have mentioned below Eq. (|5.2.6f) . even very large values of N in Eqs. 
(|5.2.4|) - (j5.2.6|) . Eq. (|5.2.5|) suggests that X c can be chosen almost equal to 
rc/4. In particular, I can evaluate the previous inequality approximatively 



by putting g t = 4A(A/r G ) 2 g s = 16 and f s = 64/r G [see Eqs. IJ5T31I ]: 



This quantity is in any case much less than one and then the gravitational 
effect is small in the present weak gravitational field approximation. Nev- 
ertheless, the effect is there and it is reasonable to imagine that it can be 
amplified in the presence of a real gravitational field which is not restricted 
by the present assumptions. In particular, we will see that this is true in the 
strong-field case treated in the following Section. 

5.3 Strong-gravitational field case 

I have said that, as a general aim of my investigation I consider the pro- 
duction of particles by a nonstationary magnetic field, so I am interested in 
situations where the gravitational effects are not the dominant dynamical 
feature. In the previous Section, we have seen that if the electron-positron 
pair is present not too close to the event horizon of a Schwarzshild black 
hole, the gravitational effects may be treated perturbatively. In order to 

4 I remind that in that case either the electron or the positron must be created in a 
state which is not a TGS. 

5 We have to be satisfied of an order-of-magnitude comparison because the time evolu- 
tion of the magnetic field used to obtain Eq. (|3.2.2p is different from that used here. 






(5.2.80) 




(5.2.81) 
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complete the investigation I suppose here that the pair production happens 
near the event horizon where, because of the singularity of the spacetime 
metric, a perturbative approach is inapplicable. The investigation is still 
possible because the isotropic metric (|5.1.ip can be approximated in a form 
[the Rindler metric [108]], where the general covariant Dirac equation is also 
solvable in the presence of a uniform magnetic field, provided the magnetic 
and the gravitational fields are parallel [1091 fllOj. In this sense, the physi- 
cal situation is different from that treated in the previous Section where the 
gravitational and the magnetic field were perpendicular [see Eqs. (|5.1.2|) 
and EH) ]. 

Following what I have just said, I want to consider here the case in which 
the pair is present microscopically speaking near the black hole event horizon 
lying at R = rc/4 [see Eq. (jS.l.ip ]. Also, I choose the reference system in 
such a way the pair is created in a volume centered on the z axis. In this 
way, the same considerations done at the beginning of the previous Section 
allow me to expand the metric tensor (|5.1.ip around the point (0, 0,rc/4). 
In particular, 



g 00 (x, y, r G /4 + z) = + O 



z 



ga(x,y,r G /'i + z) 



-16 + ( — 

JG 



(5.3.1a) 
(5.3.1b) 



with z assumed to be positive. It is clear that I am only interested in the 
pairs created in the (z > 0)-region because those created in the (z < 0)- 
region will fall into the black hole. 

Now, if I keep only the lowest-order nonzero term in g^{x } y,rc/4: + z) 
then the initial metric tensor (|5.1.1|) can be written approximatively in the 
form 



g^ u (x,y,r G /4: + z) ~gffi(z) = dia. 



2zY 
re J 



-16,-16,-16 



(5.3.2) 



This metric tensor has the same form of a Rindler metric tensor describing 
an observer uniformly accelerated in the z direction [108] Actually, the 
physical meaning of the previous coordinates is very different from that of 
the coordinates in Rindler spacetime. For example, while here the coor- 
dinate t is precisely the time coordinate in the region far from the black 
hole, the time coordinate in the Rindler spacetime is a combination of the 
Minkowski time coordinate and of the Minkowski spatial coordinate along 
the acceleration. Nevertheless, the fact that the two metric tensors have the 
same form allows me to conclude that the metric tensor (I5.3.2P describes a 



6 I could have scaled the spatial coordinates in order to have exactly the Rindler metric 
tensor, but I prefer to work with x, y and z that are the Cartesian coordinates at infinity. 



108 



Section 5.3 



constant and uniform gravitational field in the z direction. Observe that no 
assumption is needed about the strength of the gravitational field itself. 

Now, I should pass to the mathematical description of the magnetic field 
that, actually, is identical to that I have done in Sect. 15.11 and it will not be 
repeated here. 

As previously, in order to calculate the pair presence probability, I have 
to build the second quantized Hamiltonian of a Dirac field VP'(t, r) in the 
presence of the just introduced gravitational field and of the magnetic field 
(I5.1.2p . By choosing the diagonal tetrad field e^ >fl (z) with 

4^) = (5.3.3a) 
e^ l (z) = — (no sum), (5.3.3b) 

the spatial connections r (z) vanish while Fq (z) is independent of z and 
it is given by 

rf> = J^yy. (5.3.4) 

In this way, the Lagrangean density (|5.1.7f) becomes 



128z - 
mww. 

r G 

(5.3.5) 

By using the same definition used in the previous Section, the Hamil- 
tonian density of the Dirac field r) can be written, apart from total 
derivative terms, as 

J?'( R ) (t) = QW^H' iR) (i)V (5.3.6) 

with 

n ' {R) ^ = % {T [ ~ idx + eA " P (*,y)] + T [ ~ idy + eA Tv^ + M 

a z {z,d z } 



+ A 

4:1 TQ 



(5.3.7) 



the one-particle Hamiltonian of an electron in the presence of the magnetic 
field (|5.1.2p in the spacetime with the metric tensor (|5.3.2p . 
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Now, the scalar product (12.3.370 between two generic spinors \I/i(i,r) 
and *$>2(t,r) becomes here 

dr —iS\(t,r)^ ^ 2 (t,r) = 64 / drtf{(i,r)tf 2 (i,r) 
r G 2z 7 

(5.3.8) 

and the one-particle Hamiltonian 7i^ R \t) is Hermitian [this definition of the 
scalar product clarifies the presence of the numerical coefficient 64 in Eq. 
(15,3,61) ]. Finally, the total Hamiltonian of the system under study is 

H'W(t) = M j dr^(t,r)H' {R) {t)y'(t,r) (5.3.9) 

and it depends explicitly on time through the time-dependence of the mag- 
netic field [see Eqs. (|5.3.7p and (|5.1.6p ]. Since also in this case I will apply the 
first-order adiabatic perturbation theory to calculate the pair presence prob- 
ability, then in the next Paragraph I will determine the electron and positron 
modes of the time-independent counterpart of the one-particle Hamiltonian 
(lOTp . 

5.3.1 Computation of the one-particle modes and energies 

In this Section, in order to compute the electron and positron one-particle 
modes and energies, I assume the magnetic field to be static, lying in the z 
direction and, in particular, to be given by Eq. (|2.1.1|) . All the quantities 
that depend on time through the magnetic field B < f xp (t), except the vector 
potential A° xp (t, r) and the magnetic field B^ xp (i) itself that will be sub- 
stituted by A'(r) and B' respectively [see Eqs. (I2.1.4|) and (|2.1.1I) ]. will be 
indicated here with the same symbol used in the previous Paragraph but, 
of course, omitting the time-dependence. 

In the following, I limit myself to the determination of the electron one- 
particle modes U' 3 {r) and energies Wj where j embodies all the needed quan- 
tum numbers. Since the magnetic field B' is parallel to the z axis, the 
eigenvalue equation 

H'^U^WjU'j (5.3.10) 

that is [see Eq. (jO?7|) ] 
| + eA'M] + ^Hd y + eA' y ( X )] + /3m} + V, 

= w j U j 
(5.3.11) 

can be solved exactly [109J. In order to determine unambiguously the spinor 
basis, I require that the functions U' 3 {r) are also eigenstates of the conserved 
spin operator [7H I109j 

Sz = Y~8^ { ax ^~ idy + eA 'y {x ^ ~ ay ^~ idx + eA ' x{y) ^ } ' (5 - 3 ' 12) 
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It is useful to write the corresponding eigenvalue equation in the form 



SJJ' 



k , 



a 



u' 



2m 3 



(5.3.13) 



where a = ±1 and k 3 > is a real parameter that, in general, depends on 
the various quantum numbers. 

Now, Eqs. (15.3. lip and (|5.3.13p have been solved together in [109] and in 
|110j (where the effect of the anomalous magnetic moment of the electron is 
also taken into account). Actually, in those papers another electromagnetic 
gauge for the vector potential A'(r) is used but the calculations can be 
adapted straightforwardly to the present case. I only sketch the procedure 
to determine the electron modes U'j(r) by quoting the relevant steps. The 
first goal is to decouple Eq. (|5.3.1ip into a longitudinal part depending only 
on z and a transverse part depending only on x and y. This is achieved 
by multiplying Eq. (|5,3.1ip by a z (3 and by exploiting Eqs. (|5.3.12j) and 
(I5.3.13p . The resulting "longitudinal" equation is 



Aakj — io z ^i d z 



. PVT. 
2z 



z ) 



0. 



(5.3.14) 



Since a z commutes with 7 and 7 s and since {7°, 7^} = 2rj a ^ , if I square 
the previous equation I obtain 



d z + ^d z - 16k 2 - Q - 2iw J r G a z - 4^ 2 r^ 

In order to satisfy this equation I write the spinor U'^r) as 
Ufa) = N 3 [P-ML^z) + P+M'+^z)] Z' 3 (x,y) 



u; = 0. 



where N 3 is a normalization factor, 



P± 



(5.3.15) 



(5.3.16) 



(5.3.17) 



are two 4x4 projectors, H^(x, y) is a spinor depending only on the transverse 
coordinates and Mj_ Az) are two functions to be determined. By substituting 
Eq. (|5.3.16p in Eq. (|5.3.15|) I obtain the following equations for the functions 
M' ±J (z): 



d 2 ld_ 

dz 2 z dz 



whose general solution is given by 

M^ tJ (z) = a I I l/2 ±2i W3 r G {^ J z) +a K K 1/2±2iWjr . G (4k :) z) 



(5.3.18) 



(5.3.19) 
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with and K\(£) the modified Bessel functions [89]. Now, the functions 
diverge exponentially as £ — > oo while the functions K\(£) go to zero 
exponentially in the same limit. Now, my model is reliable only for z <^tq 
then I have to choose the solution that coherently is very small in the region 
of large z. For this reason I put aj = and ax = 1 in Eq. (15.3. 19|) and I 
obtain_| 

M^(z) = K 1/2±2iWjra (Ak 3 z). (5.3.20) 

In order to determine the transverse spinors r! 3 {x,y) I observe that 
[P±,S Z ] = 0, then by substituting the spinor (|5.3.16|) in Eq. (|5.3. 13|) I 
have 

[Ama z - if} [a x [-id y + eA' y (x)] - a y [-id x + eA' x (y)] }] r! d = 4akjE' r 

(5.3.21) 

I observe that the energy eigenvalue Wj does not appear in this equation 
then, looking also at Eqs. (15.3. 18j) and (I5.3.20p . any continuous value = 
E > is acceptable. The fact that in the presence of the gravitational field 
described by the Rindler metric tensor §5. 3.2\) the energy of the electron has 
continuous eigenvalues from zero to infinity that do not depend on the other 
quantum numbers is the most relevant difference with the case in which no 
gravitational field is present. The physical origin of this difference lies on the 
fact that in the present case the linear momentum along the gravitational 
field is not a constant of motion. In other words, the "longitudinal" energy 
of the electron contains not only, as in absence of the gravitational field, the 
rest energy and the kinetic energy but also a negative gravitational potential 
energy. 

In order to solve Eq. (I5.3.2ip I square it, then 




(5.3.22) 

with C z the z component of the electron orbital angular momentum. The 
solutions of this equation are well known [7'6\ [73] . Two nonnegative integer 
quantum numbers nd and n g have to be introduced and 

kj = k nd = \j m2 + ^ L ( 5 - 3 - 23 ) 

can be interpreted as a sort of "transverse" energy of the electron in the 
spacetime with the metric (I5.3.2p . With this definition the spinor S'Jx,y) 



7 I point out that, since the behaviour of the functions /a(£) is regular near 5 = 0, they 
will be used again later [see Eq. (I5.3.36P ], 
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is given by 



"j(. x > V) ~ ^n d ,n g ,tj( x i y) 



2\/k 



f ^/k nd + ma6' nd _ ltng (x,y) ^ 
y/k 

icry/Kj + ma9' nd _ l ng (x, y) 



\ ia^k nd -mo6' ndng (x,y) J 

(5.3.24) 

where the functions 9' ndJlg (x,y) are given in Eq. (I2.1.28p . The functions 
6'_ l ng (x, y) are not defined but this does not cause any problem because, by 
solving step by step Eq. (|5.3.22|) . one finds that if = then a must be 
equal to —1 and the corresponding coefficient \/k Ud + ma vanishes. Finally, 
I quote that the coefficients in Eq. (|5.3.24|) have been chosen in such a way 
the spinors E' nd ^ ng ^(x, y) result normalized as 



/ 



dxdyZ'£ d)ng a 



(5.3.25) 



At this point I have to determine only the normalization factor in Eq. 
(15.3.16p . As I have said, E is a continuous eigenvalue. For this reason if I 
require that the functions U' 3 {r) = U' ndn a (E;r) are normalized as [see Eq. 

K d ^ a (E), U' n>>l (E'))^ = 5(E - E')6 ndtn , d 6 ngK 6 a y (5.3.26) 

then, it is easy to show that the final form of the electron modes of the 
one-particle Hamiltonian is [see [109J for a more detailed derivation of the 
normalization factor! 



*4 (£;r) 



k n j G cosh(27r£r G ) 



4vr 2 

x [P-K 1 /2 + 2iEr G ( 4k n d z) + P+K 1/2 _ 2iErG (4k nd z)]E' nd ^ g)(T (x,y). 

(5.3.27) 

The physical meaning of the quantum numbers E and a is clear from 
Eqs. (|5.3.10p and (15.3. 131) . In order to understand the physical meaning of 
the remaining quantum numbers and n g I introduce the z component 



(1/2) 



of the total angular momentum operator and the operator 

2 / -n \ 2" 



R 



xy 



16 



eB 



1 2 eB 



(5.3.28) 



corresponding to the operator (|2.1.13l) in Minkowski spacetime. Now, it can 
easily be shown that the two previous operators commute between them, 
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with the time-independent form of the one-particle Hamiltonian (|5.3.7f) and 
with S z . Also, it can be shown that 

rt l,2) K d ,n g AE) = L d -n g -^j K d>ng>a (E), (5.3.29a) 
R'U'JE) = n2n C l) U'A E )- (5.3.29b) 



Moreover, in the following I will use the operator p^ y that is defined as 



2 
Pxy 



16 



X Vy 

2 eB 



+ 



2 eB 



(5.3.30) 



and that corresponds to the square of the radius of the helix along which a 
classical electron moves in the presence of the magnetic field B' given in Eq. 
(I2.1.ip and in the spacetime with metric tensor (|5.3.2I) [see |73j]. It can be 
shown that the states (|5.3.27|) are not eigenstates of , but that if S> 1 
then 



Pxy^n d ,n g ,a( E ) — g ^ ^n d ,n g ,a( E ) 



n d > 1. 



(5.3.31) 



In conclusion, the previous eigenvalue equations allow me to conclude that, 
as in Minkowski spacetime, the quantum numbers and n g are connected 
with the motion of the electron and of the positron in the plane perpendic- 
ular to the magnetic field. 

The positron modes can be built in an analogous way and the final result 

is 

r) 



kn g rc cosh(27r-Er G ) 



4vr 2 

n g z y\" 'n d ,n g ,o(% •> V) 

(5.3.32) 



with 



( iay/K g - ma6' ng _ hnd (x,y)\ 

x,y) 



2 -\ / kit-) 



1(7 



\JKi g + maO. 



n g ,n d \ 

y/k^-maO ' 1>nd (x, y) 



(5.3.33) 



V \/ k n g + rna9' ng , nd (x,y) J 

These modes satisfy the eigenvalue equations 

n^V^ ngta (E) = -EV^ ngia (E), 
k 

S zVn dt n 3 A E ) = ~ a ^ V n d ,n g ,a( E ), 



(5.3.34a) 
(5.3.34b) 
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J z {1/2) V^ ng ,AE) = -[n d -n g + V^ ng>a {E), (5.3.34c) 
Rl y V^ a (E) = n2n " + 1) V^ ding>a (E), (5.3.34d) 



eB 
g 

eB 



32n 

PlyVn d ,n 9 A E ) ~ -JTT V n d ,n 3 ,A E ) n 9 » 1 (5-3.34e) 



with the constraint that if n g = then a = +1 and they are such that 

{K d ,n a A E )X' n> AE')) {R) = 0. (5.3.35b) 



Finally, in Appendix E I show that the set of spinors U' nd ^ ng ^{E;Y) and 
V n d ,n a ,A E '> r ) is complete. 

As usual, it is preferable to deal with normalizable modes then I have to 
find a convenient boundary condition at a given surface z = b that discretizes 
the energies E. Since the procedure is identical for the electron and the 
positron modes, I will consider only the electron modes. Now, as I have said, 
the functions Ki/ 2 ±2iEr G (^ n d z ) S° exponentially to zero for large values of 
k nd z and go to infinity as (k nd z) ~ 1 ^ 2 for small values of k nd z |89| . For this 
reason, it is clear that 



1. the modes U' ndrLg (J {E]r) can not satisfy a "zero" condition at a given 
kn d b <C 1 or a canonical periodicity condition between two points 
k nd bi < 1 and k nd b 2 » 1; 



2. if we want to build modes with a finite normalization integral we have 
to modify the functions K 1 /2±2iEr G {^k rLd z) in the region with k rid z <C 
1. 



On the other hand, we already know that the modified Bessel functions 
I 1 /2±2iEr G {^kn d z) have a regular behaviour in the region k nd z <C 1 and that 
they also satisfy the longitudinal equation (|5.3.18[) [see Eq. ([5.3. 19j) ] . For 
this reason, I consider an arbitrary fixed surface z = b such that k nd b 1 
and assume that the electron modes of the one-particle Hamiltonian and of 
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the spin operator (|5.3. 12f) are the spinors Uj(r) defined as 



U'j(r) 



N 



«) 



K d r G cosh(27r£ , n|rad r G ) 
4tt 2 

P-h/2+2iE n , nd r G (^k nd z) 



l/2-2iE n>n r G 



(4k 



"n d ,n a 



Jx,y) \iz<b 



N 



(» 



k nd r G cosh(2irE nind r G ) 



4vr 2 

P - K l/2+2iE n 
+P+K\/2-2iE n 



r G ( 4 *V, 



, ro 



(4k 



'id' 



^n d ,n g 



if z > b 

(5.3.36) 

where J = {n, rid, n g, °} with n a new integer quantum number characteriz- 
ing the discrete energies (as I will see the discrete energies will also depend 
on the quantum number rid) an d where Nj and Nj are two real normal- 
ization factors to be determined. It is evident that the spinor Uj(r) satisfies 
Eqs. (|5.3.1ip and (|5.3.13j) in the regions z < b and z > b [see also Eq. 
(I5.3.19|) ]. Also, since Eq. (15.3. lip is a first-order equation in the variable z 
I only require that the spinor Uj(r) is continuous at z = b. By means of 
this condition and by requiring that the norm of Uj(r) is unit, I make the 

energies discrete and determine the normalization factors Nj and Nj- 
The details of the calculations are given in Appendix F and here I only quote 
the final expression of the discrete energies [see Eq. (jF.16j) ] 



En. 



I'd 



7VK 

2r^ 



log 1 (k nd b) 



k nd b 







(5.3.37) 



and of the coefficients iYj <} and iV) >j [see Eqs. (TOaj) and (fETHl) ]: 



7T y/1 + (4E n , nd r G ) 2 1 



N ( T <] = iV«) - 
3 n ' nd 8k nd b cosh(27r E n>nd r G ) 



Ar(» = N (» 

J n d 



1 



k nd b -> 0, 



K d b -» 



(5.3.38a) 
(5.3.38b) 



with [see Eq. (lFl7j) ] 



Qn d = -—log(k nd b) 

7T 



k nd b 



(5.3.39) 



the density of the energy levels. Obviously, all these quantities will be used 
in the intermediate calculations but at the end I have to perform the limit 
b — > and the physically relevant results must be independent of b. 
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5.3.2 Calculation of the presence probability 

We already know that, in order to obtain in the framework of the adia- 
batic perturbation theory the presence probability of a pair in the presence 
of the slowly- varying magnetic field (|5.1.2|) . I have to calculate the corre- 
sponding presence matrix element. If the pair is present at time t with the 
electron in the state Uj(t,r) and the positron in the state Vj,(t,r) with 
J = {n, rid, n g , a} and J' = {n' ,n' d ,n' g ,a'}, the presence matrix element is 
given by 

H'}f)(t) = (JJ>(t)\H>W(t)\0(t)) 

16e£? T cxp (t) f ,__ rr/t ^ ^ . (5 ' 3 - 40) 



re 



J drU'J(t,r)z(xa y - ya x ) Vj,(t,r) 



[see Eqs. (I5.3.9p . (|5.3.6f) and (I5.3.7P ], A more useful form of the previous 
matrix element can be given by using the matrices a± and the expressions 
(|2.1.21a|) and (|2.1.21c|) of the operators x and y: 



H'Sit) 



32ieB^ p (t) 
r G J2eB^(t) 



J drU'](t,r)z[a-(a d + a\) - a + (a g + a\)]Vj,(t, r). 

(5.3.41) 



Now, as I have said at the end of the previous Section, in order to cal- 
culate these matrix elements I should use the expression (|5,3.36p of Uj(t,r) 
with Nn^n d (t) and N^\t) given by Eqs. (|5.3.38j) and an analogous expres- 
sion of V},(t,r)E Actually, an easy power counting shows that the contribu- 
tion of the integral on the variable z from to b goes to in the limit b — > 0. 
In fact, each spinor contains a factor [k nd (t)b-\/\og(k nd (t)b)]~ 1 coming from 

N n%(t) [see Eqs. (|5.3.38a|) and (|5.3.39l) ]. Also, from Eq. (|R2al) we see that 
the modified Bessel functions I\/2+2iE n nd (t)r G (^ri d (t)z) behave as \/k nd (t)z 
in the integration domain < z < b. Finally, because of the presence of 
the z factor in the matrix element (|5.3.4ip the result of the integral on z 
depends on b as k nd (t)blog~ 1 (k nd (t)b) and then it goes to zero in the limit 
6^0. In this way, since at the end of the calculations the limit b — > has to 
be performed, the matrix element (|5.3.4ip can be calculated by using in the 
whole region z > the expressions of the spinors U'j(t, r) and Vj,(t, r) valid 
in the region z > b. Actually, I can use directly the spinors U' nd ^ ng ^{E;r) 

and V , , ,(E';r) multiplied by N^\t) and (t) respectively because 
the presence of the factor z in the matrix element (|5.3.41|) makes finite the 



8 I remind that these quantities are now time-dependent because the magnetic field 
depends on time. 
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resulting integral from to oo [see also the general formula (IF.10|) ]: 



H [H) , , ,(E,E';t)- 



x J drU^ ng ^(E;t,r)z[a.(a d + oj) - a + (a g + a j d )]V^ n ^ a ,(E' ;t,r). 

(5.3.42) 

At this point, I have to substitute Eqs. (I5.3.27P and (15.3.320 [with the time- 
varying magnetic field B^ xp (t) instead of B] in place of Un d ,n g ,cr{E; t, r) and 
V i i ,{E'\ t, r) and apply the various operators. By using the intermediate 

n d' n g' <T 

matrix elements 



/ [kn d (t) — ma] [k n > (t) — ma'] 

= (1 ± ia)(l T ia ) \ — . . , g . . Ly_ii„ „/ , 

V A ; V Mk nd (t)k n , g (t) nd ^ 'W 

(5.3.43a) 

M . W1 . . I[ k n d (t) + ma] [k n ,{t) + ma'\ 

= lTKJ 1±W W r— ; „'i „' 

; V uk nd {t)k nlg {t) 

(5.3.43b) 



and 



y dxdy~!^ d Jl g)U {t, x, y)P±a-Z' n ^ n , ga ,(t, x, y) 



(1 ± «t)(1 t «r ) t / ^ ^ 5 ndjn , g 6 ngK , 



Q4k nd (t)k K (t) 



(5.3.44a) 



64fc nd (*)£;„, (t) 

(5.3.44b) 



(1 T »<7)(1 ± ia')\l — ^ <Wn^n g ,n d 
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that can be easily checked, the matrix element (I5.3.42p can be written as 

ieB^ p (t) /cosh(27r£V G ) cosh (2irE'r G ) 



H ,{R) , , ,(E,E';t) 

n d ,n g ,a;n' d ,n'a'\ ' ' / 



7T^ 



2[eB° xp (t)] 



+ 



(ra g + 1) [k nd (t) - ma] [k nd+1 (t) - ma'] 
Qn d {t)en d +i{t) 

Re [(1 - ia)(l + ia')I nd , nd+ i(E, E';t)] 5n d ,n g -i5 ng +i,n> d 
/ n d [k nd (t) - ma] [k nd (t) - ma 1 ] 



QUt) 



Re [(1 - ia)(l + ia')I nd , nd (E, E';t)] S nd y g S ng 
jUg [k nd (t) + ma] [kn d -i(t) + ma'] 



Qn d {t)Qn d -l{t) 

Re [(1 +ia){l - ia')I nd , nd ^(E,E';t)] 6. 



n d -l,n' g On g -l,n' d 



l n d [k nd (t) + ma] \k nd (t) + ma' 



Re [(1 + M7)(1 - ia')I ndtnd (E,E';t)] 5 nd y5 ng y 



(5.3.45) 



where the adimensional function 
Il,l>(E,E';t) 



dssK l / 2 -2iE 



'jTQ 



4fcj(t) 



2eB^ p {t) 



K 



l/2+2iE'r G 



4ki>(t) 



2eB^ p {t) 



(5.3.46) 



has been introduced. 

Before continuing, I want to point out that from Eq. (|5.3.45l) it can be 
seen that the total angular momentum of the field ^'(^r) is conserved in 
the transition. In fact, in any case [see Eqs. (|5.3.29ap and (15. 3. 34c |) ] 



1 / / 1 / 

nd ~ n g - - + n d - n g + - = n d - n g + n d 



n' g = 0. 



(5.3.47) 



Of course, this selection rule is a consequence of the fact that the time 
evolution of the magnetic field does not break the rotational symmetry of 
the system around the z axis or, in other words, of the fact that Jz and 
H'( R '(t) commute. 

Now, since I am interested only in the strong magnetic field regime in 
which B^ xp (t) 3> B cr , I can simplify the expression of the transition ma- 
trix element (|5.3.45p by taking into account only those transitions whose 
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probabilities are proportional to the lowest power of B cr /B^ xp (t). In the 
framework of the adiabatic perturbation theory the first-order transition 
amplitude in B° xp (t) of the presence of a pair at time t in the state with 
quantum numbers {E, n^, n g , cr; E', n' d , n' g , a'} is given by 



7 W , , AE,E';t) 

i r t (5.3.48) 
= — - — / dt'H' {R) , , ,(E,E':lf)exp[i(E + E')t'} 

E + E' Jq n d ,n g ,a;n d ,n' g ,a'\ ' ' J 

and the corresponding probability is the square modulus of this quantity. 
It is evident that, since the energies E do not depend on B^ xp (t), we 
can perform the (B cr /B^ xp (t))-powei counting directly on the matrix el- 
ement (]5.3.45p . To this end I need the general behaviour of two particular 
classes of the integral (|5.3.46p that is of Iq^i (E, E';t) with n' d > and of 
I ndjTld (E , E' ; t) with rid > 0. By reminding the expression (|5.3.23p with the 
time-dependent magnetic field for k nd {t) and by using the general formula 
(IF.lOh . it can easily be seen that 



Io, n > 4 {E, E'; t) ~ ^ — if n ' d > o and B™ p (t) » B cr , (5.3.49a) 

In d ,n d (E, E'; t) if n d > and B° xp {t) » B cr (5.3.49b) 

where, for later convenience, I also pointed out the dependence on the quan- 
tum numbers and n' d . Obviously, the integrals I nd ,n d ±i(E, E'; t) behave 
as the integral (|5.3.49b|) and then this criterion allows me to neglect the tran- 
sitions in which the electron is in a (rid = 0,a = — l)-state or the positron 
in a (rig = 0, a = +l)-state [see Eq. (I5.3.45I) ]. 

Another criterion I will use to select only the most probable transitions is 
the dependence of the corresponding probabilities on the quantum numbers 
rid an d rig. As previously, I can work directly on the matrix element (|5.3.45p 
by keeping in mind that at the end I will sum the probabilities with different 
values of n d and n g . Now, according to what I have said at the end of Chap. 
3, the internal consistency of the model requires here that the sum on n g 
(and on rid) cannot be extended up to infinity but that they must be stopped 
up to a certain N^ xp (t) corresponding through the relation 

Nr(t) - ~%^R1m (5-3.50) 

to a fixed R ±M [see Eqs. (|5.3.29bl) . (I5.3.31P . ([5.3.34dl) and d5.3.34eP ], By 
reminding the physical meaning of the operators R\ y and p%. y , it is clear 
that the transverse motion of a classical electron (positron) is confined in a 
circle with radius 2R±m, in such a way this quantity can be assumed as the 
radius of the quantization cylinder whose axis is parallel to the z axis. 
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Coming back to the matrix element (|5.3.45p . we see that it contains two 
kinds of terms, the first one being proportional essentially to ^Jn^ and the 
second one to y/nd. By taking into account Eq. (15.3.49bj) it is easy to see 
that the first kind of terms gives rise to final probabilities proportional to 
(iV| Xp (i)) 2 log iVf xp (t) , while the second one to final probabilities propor- 
tional to (iV T exp 0)) 2 . For all these reasons I can consider only the transitions 
to states with > and n' g > in such a way only the following four kinds 
of transitions amplitudes result different from zero: 

7n d ,n 9 ,cr;n d +l,n 9 +l,cr(^'i & i 




(n g + l)\Jn d {n d + 1) cosh(27r£'r G ) cosh(27r.E / r G ) , 
Qn d (t)g nd+1 (t) 
Re(I nd: n d+ i(E,E')) , 

x eTe> F (^>*)> 



ln d ,n g ,(r;n d ^-l,ng-\-l,-a^-> E ! *) 



id / (n fl + 1) \/ n d( n d + 1) cosh(27r£V G ) cosh(27ri? / rG) 

= ^V Qn d (t) end+l (t) ( 5 - 3 - 51b ) 

Im{I nd>nd+1 (E,E')) 

x eTe' F ^ E ^h 

7n (i +l,n 9 +l,o-;n d ,n 9 ,o-(-^'' ^ ' 




ra 5 + l)\/rid(nd + 1) cosh(27r £V G ) cosh(27r£ 7 r G ) 
6n d (t)Qn d +l(t) 

^ iIn f^ E>)) F(E,E'-,t), 

(5.3.51c) 

„(«) 




n 5 + l)\/nd(nd + 1) cosh(27rii>G!) cosh(27i\E'r G ) 
6n d (t)Qn d +l(t) 



where 



(5.3.51d) 

f* B cxp (t') 

F(E,E';t) = ^ dt'^^exp[i(E + E')t'} (5.3.52) 

and where I pointed out that in the strong magnetic field regime, if > 
the integrals I nd ,n d ±i{E, E') do not depend on time. By squaring these 
amplitudes, by summing on the polarization a and by multiplying by the 
number of electronic states g nd (t)dE with energies between E and E + dE 
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and by the number of positronic states g nd +i(t)dE' with energies between 
E' and E' + dE', I obtain the differential probabilities 



dR 



n d ,n g ;n d +l,n g 



n g+ i(E,E';t) 



Y~i( n a + 1) \Jn d {n d + 1) cosh(27r£V G ) cosh(27r£V G ) 



(5.3.53a) 



\l 



n d ,n d +l 



(E,E>) 



A2 



{E + E') 



\F(E,E';t)\ dEdE', 



dP i d +l,n g +l;n d ,n g ( E i E ' 'i *) 



2vr 4 



(n 9 + 1) \Jn d {n d + 1) cosh(27i\Er G ) cosh(27r£;V G ) (5.3.53b) 



I/, 



n d +l,n d 



(E,E>) 



(E + E>) 2 



\F(E,E';t)\ 2 dEdE' 



that, as expected, do not depend on the unphysical parameter b. Now, 
I want to calculate the probability dP^ R \E , E' ; t) that a pair is present 
at time t with the electron with energy between E and E + dE and the 
positron with energy between E' and E' + To do this I have to sum 
on the remaining quantum numbers n d and n g . As we already know, both 
the series on n d and n g are diverging, then I can perform the summations 
by assuming n g ~ n g + 1 and n d ~ n d + 1 because the most relevant terms 
are those with n g 3> 1 and ^ 1. Starting from Eqs. f|5.3.53|) we have 



dP {R) {E,E'-t) 



7T J 



'7V T exp (t) 



cosh(27r Ere) cosh(2ir E'r G ) 



E 



iV T exp (t) 
n d =\ 



n d \I nd n .(E,E') 



{E + E') 2 



\F(E,E';t)\ dEdE' 
(5.3.54) 



where N° xp (t) has been defined in Eq. (|5.3.5Uft . 

The next step is the explicit calculation of the functions I ndt n d (E, E') 
and F(E,E';t) defined in Eqs. (|5. 3.46ft and (|5. 3.52ft . By using the general 
formula (jF.lOft and the properties of the T function (|F.5ft and 



r(i + iOr(i-i^) = |r(i + ^ 



Tit; 



with £ € 



sinh(-7r£) 
it can easily be shown that 

7T 1 - i(E - E')r G tt(E + E')r G 



(5.3.55) 



In d ,n d (E,E') 



2n d cosh [tt(E - E')r G ] sinh [ir(E + E^g] 



(5.3.56) 



122 



Section 5.3 



Instead, after some calculations the following expression of F(E,E';t) 
can be obtained: 

F(E,E';t)= / ds'- — ^ r~~ -. — - exp \-s' + i{E + E')ts'] 

Jo + exp (-s') J 

(5.3.57) 

with s' = t' )t. Now, as usual I am interested only in energetic electrons 
and positrons such that Et ^> 1 and E't 3> 1 and in large times. For this 
reason, I can give the following asymptotic estimate of the integral (|5.3.57|) : 

By substituting this expression and Eq. (|5.3.56j) in Eq. (|5.3.54j) , I can write 
the asymptotic value of the probability dP( R > (E, E' ;t) as 



r, 



G" 



r G [(£ + £'> G ] 2 v ' 

(5.3.59) 

where I made the substitutions valid for large limes [see Eq. (|5.3.50p ] 

1 , ^fiirrw.^O 1 / eBfR 



^ n 9 = - 2 {^{t)?= l -(^l^L\ , (5.3.60a) 

JVf p (t) 

£ i =1 °g( Ar t XP W)=log( £ ^ M ) (5.3.60b) 



and where the function 

^(rg)f^ E >) = gZg cosh(27r£r G )cosh(27r£'r G ) 

1 ' ; 2 cosh 2 [7r(£-£> G ]sinh 2 [7r(£ + .E> G ] ' ' 

has been introduced. I pointed out the dependence of 6^ rG \E,E') on the 
Schwarzshild radius r G because from a physical point of view I am interested 
in energies E and E' such that Ere 3> 1 and E'r G 3> 1. In this energy 
region the function 5^- rG \E, E') strongly depends even on small changes of 
E and E' through the hyperbolic functions. This can be seen more clearly 
by writing Eq. (|5.3.61l) as 

S (ra) (E E >) = nr G cosh[2^(E - E')r G ] + cosh[2^(£ + E')r G ] 
1 ' ' 4 cos\i 2 [k(E - E')r G ]s\nh 2 [iT{E + E')r G ] 



irr G 



+ 



2 [ sinh 2 [ir(E + E')r G ] cosh 2 [ir(E - E')r G ] 



(5.3.62) 
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From this expression and by reminding that E, E' > 0, I obtain 
lim 5 {rG) (E,E') 

lim <! — „ ~ — — — - + 



r-c-oo 2 \smh 2 [Tr(E + E')r G ] cosh 2 [it {E - E')r G ) ) (5.3.63) 

vrr G [O ftE^E' 

= lim ^ = < 

r-c^oo 2cosh 2 [vr(£ - E')r G ] loo if E = E' . 

Finally, by observing that 

vr f°° r G dE' f°° dq 



r = 1, (5.3.64) 

JO COSh 77 



2 7-oo cosh [tt(E — E')r G ] Jo cosh r\ 
I can conclude that 

lim <5 (rc) (£, = S(E - E') (5.3.65) 

and then that 

<5 (rc) (£?, ~ 5(£ - £') if > and Er G , E'r G > 1. (5.3.66) 

With this result and by integrating Eq. (15.3.591) with respect to the positron 
energy E' , I finally obtain the probability that an electron is present at large 
times with an energy between E and E+dE such that Er 3> 1 and Er G ^> 1 
in the form 

(5.3.67) 

In order to obtain a probability per unit volume I have to give an estimate of 
the height of the quantization cylinder. Now, I have said that the modified 
Bessel functions Ki/ 2 ±2iEr G (^n d (t)z) (I refer to the electron wave functions 
but an identical conclusion can be drawn for the positron ones) are exponen- 
tially decreasing as k nd (t)z ^> 1. In particular, it can easily be shown that 
if Er G 3> 1 the exponential behaviour of the function Ki/ 2 ±2iEr G (.^n d (t)z) 
for large times starts at 

2Er G 2Er G 
z = — — < G . (5.3.68) 



For this reason I can assume L z = 2Er G / y / 2ei?/ and 



V = <2R±M) 72H7 = ^Bj (5 - 3 - 69) 
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as the volume of the quantization cylinder in such a way 



d^dE vr V 128 J to V 32 / V rB i J E 3 

This probability does not depend on the electron mass nor on the gravita- 
tional radius of the black hole but this is due only to the fact that I am 
working in the strong magnetic field regime and in the high-energy region. 
Instead, it is not so obvious how to explain the dependence of the final pres- 
ence probability (|5.3.70p on the logarithm of R±m- Now, before comparing 
Eq. (|5.3.70p to the analogous result Eq. (13.2.2)) where no gravitational field 
effects were taken into account I have to integrate Eq. (|5.3.70p with respect 
to the electron energy. Nevertheless, I stress that also in this case the com- 
parison can be only qualitative because a linear dependence on time of the 
magnetic field was used to obtain Eq. (I3.2.2p . Before doing that, I note that 
the presence probability scales here as E~ 3 in the high-energy region while 
in the preceding case the corresponding probability per unit longitudinal 
momentum k behaved as A; -4 [see the details in [88]]. This means that the 
production of high-energy electrons (positrons) is favoured in the presence 
of the gravitational field. Now, to be coherent with the approximations I 
have made, I perform the integration of Eq. (|5.3.70p from E m = 100 r^ 1 to 
infinity in fact, by assuming r = 1 s and re = 3.0 x 10 6 cm as for a 10 solar 
masses black hole, it also results E m r 3> 1. By indicating the resulting total 
probability per unit volume as dP^ R \t)/dV, I obtain 

dpC'Ht) 1 f eBf \ 5 ^ 2 . / eBfR 2 LM \ /B/-B.V R\ M 



In general, the presence probability per unit volume depends here on the 
(5/2)-power of the magnetic field strength while in Eq. (|3.2.2j) it depended 
on its (3/2)-power [remind that in Eq. (l3~2~2]> b = Bf n (i)]. To give a 
more quantitative estimate I use the typical values R±m = 10 5 cm, Bi ~ 
Bf = 10 15 gauss and I assume the magnetic field (13.2. 1 f) to be such that 
B^ n (t) ~ Bq = B{ and b = {Bf — B{)/t. In this way I obtain [see also Eq. 
(l3X2j) ] 

dP^(t) El _g 7xl0 -29 (5 3 72) 



dP( R )(t) eB f 

that clearly allows me to conclude that the effects of the gravitational field in 
the pair production process are really relevant and they can not be neglected 
at all. 



5.4 Summary and conclusions 



This Chapter has been devoted to the study of the effects that the presence 
of a gravitational field can have on the production of electron-positron pairs 
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in the presence of a strong, uniform and slowly- varying magnetic field. The 
motivation of this analysis comes from the fact that in the model I have in 
mind the pair production is assumed to happen near astrophysical compact 
objects such as magnetars or black holes. Then, it is reasonable that, espe- 
cially in the last case, the gravitational effects can be relevant. Two different 
physical situations have been treated here: in the first one the pair is imag- 
ined to be produced far from the black hole event horizon (Sect. I5.2|) and 
in the second one it is assumed to be produced near the black hole event 
horizon (Sect. I5.3p . In both cases the presence of the gravitational field 
has been taken into account only in the determination of the electron and 
positron one-particle modes and energies while the calculation of the pair 
presence probabilities has been performed by using the first-order adiabatic 
perturbation theory. 

In particular, if the pair is produced far from the black hole event hori- 
zon, the effects of the gravitational field have been treated perturbatively. 
I have shown how the modified one-particle modes and energies of the elec- 
tron and of the positron reflect on the pair presence probabilities [see Eq. 
(|5.2.79|) ]. In particular, I have examined the case of the production of a pair 
in the presence of a magnetic field varying only in strength and always per- 
pendicular to the gravitational field. Firstly, I have found that the presence 
probability contains a factor g^ 1 that makes it growing and growing as one 
gets closer and closer to the event horizon of the black hole. More impor- 
tant, I have also obtained a new new qualitative result: even in the presence 
of a weak gravitational field, even if only the strength of the magnetic field 
changes with time it is possible to create a pair with both the electron and 
the positron in a TGS. Actually, this probability is a small quantity with re- 
spect to the total probability that a pair is created in Minkowski spacetime 
in the presence of a time- varying magnetic field with fixed direction [see Eq. 
(|5.2.81j) ]. but this result is a consequence of the fact that the gravitational 
field has been treated perturbatively. From this point of view, the informa- 
tion we have gained is that in the presence of a gravitational field this new 
effect is there. 

Instead, if the pair is produced near the black hole event horizon, the 
situation is very different mostly because the one-particle energy of the elec- 
tron (positron) is an independent continuous nonnegative quantum number. 
Also in this case, I have considered a situation in which the magnetic field 
does not change its direction with time but, here, always remaining parallel 
to the gravitational field. In this case [see Eq. (|5.3.70l) ] we have seen that 
the presence probability depends on the (5/2)-power of the magnetic field 
strength while the analogous quantity in Minkowski spacetime depended 
only on its (3/2)-power. Also, the presence probability scales here as 
in the high-energy region while in the flat-spacetime case the probability 
behaved as k~ A with k (—k) the longitudinal momentum of the electron 
(positron). In this way, the production of high-energy electrons (positrons) 
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is strongly favoured in the presence of the gravitational field. Moreover, the 
ratio (|5.3,72p allowed me to safely conclude that the effects of a strong grav- 
itational field in the pair production process are dramatically important and 
they can not be neglected at all. A final observation concerns the fact that 
the presence of a strong gravitational field makes possible the creation of 
pairs that can not fly to infinity because they do not have enough energy 
(as I have said, the one-particle energy spectrum of the electrons and of 
positrons extends down to zero). These electrons (positrons) created with 
such energies annihilate inside the gravitational field producing low-energy 
photons which may fly away and, eventually, contribute to the low-energy 
part of the GRBs spectra. 
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In this Appendix I study the general features of the presence amplitudes 
of pairs in which the electron and/or the positron is not in a TGS in the 
presence of the rotating magnetic field (|3,2.3|) . Firstly, I write the corre- 
sponding presence matrix elements by using the compact notation Hjji{t) = 
(jj' (t)\H(t)\0(t)} [see Eq. (|3.1.14|) ]. In this way since in the presence of a 
purely rotating magnetic field these matrix elements and the energies Wj(t) 
and Wj'(t) are actually time- independent, the corresponding amplitudes are 
given by [see Eq. (|3.1.13|) ] 



Tii'W = 2- —2 

( Wj + Wj,) 



exp 



- (wj + Wj,) t 



sm 



2 ( 



- i + WjA t 



(A.l) 



The presence matrix elements different from zero are 



Hn d ,k,a,n g ;n g +l, — k,—a,n d (t) ^n^fcjO^ng+lfrig,— k,— <7,n,j CO 



crQN. 



33 



eB 7 



(n g + 1) 



l + a 



W 2 - k 2 



33 



+ 2n d + 1 



Hn d +l,k,a,n g ;n g ,~k,—a,n d {t) — H rld ] tang - ng .— i t — and J r \(t) 



oVLN jr \l^{n d + l) 



3 + 0- 



W 2 - k 2 



33 



eBr 



+ 2n d + 2 



-QNjj'\/ —^-{n d + 



Hn d ,k,a,n g ;n g ,k' ,(T,n d (0 



nN jf 



(A.2a) 



(A.2b) 



Hrifa k, (T,ng;n g , k', — CTjU^+l (0 Hn d +l,k,a,n g ;n g ,k' ,—a,n d (t) 

FTbT, ~2TTAk, {Lz) ( A - 2c ) 



<t (fctf - W 2 f ) ^4 Lz) - ik (2n d + l)5 K . k , 



H nd ,k,-l,n g ;n g +l,-k,-l,n d -l(t) — Hrid,k — l,ng+l;n g> — k,— l,n<j— 1 (0 



-2ikQNjj>\/n d (ng + 1), 
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(A.2d) 



(A.2e) 
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H nd ,k,+l,n g ;n g +l,-k,+l,n d +l{t) — -ffn d ,fe,+l,n 9 +l;n g ,— fc,+l,n d +l (t) 

(A.2f) 



= 2iknN jf J(n d + l)(n g + 1), 

H nd ,k,+l,n g ;n a ,k' ,+l,n d +2(t) = H nd+ 2,k,-l,n a ;n g ,k> ,-l,n d (t) 

= -2ClN jf (ik5 K _ k , + eB;^5f^ y/(n d + 2)(n d + 1). 



(A.2g) 



In order to simplify the previous formulas I defined the quantities 



Ak = k' - k, (A.3a) 

K = k' + k, (A.3b) 

W?-/ = (toj + m)(uy + m), (A. 3c) 

N 3jl = ^,l^± m ^'t m . (A.3d) 

•" 21l j ; . \ 2,r ; 2w f V 



Also, by reminding the quantization conditions ([2.1. 3Uj) . I note that the 
'k 



function S'jf"^ in Eqs. (IA.2jl defined as 



AL z) If*. , , 10 if A = 



= — / dz zexp(-ikz) = { e (A.4) 

becomes the derivative of the 5 function in the limit L z — > oo. 

All the previous matrix elements can be divided into two groups: the ones 
related to transitions in which the longitudinal linear momentum conserves 
and the others characterized by the presence of the function S'jf^ in which 
it does not. Of course, this fact is due to the dependence of the transition 
operators T'-, y (t) on z [see Eq. ([3.2. 5p ]. 

Finally, I observe that if one sums the probabilities corresponding to the 
previous matrix elements by means of Eq. (jA.lj) with respect to the quantum 
number n d , all the series converge. Only the series corresponding to the 
matrix elements (|A.2bj) diverges logarithmically and it is not so obvious 
how to give a physical interpretation of such a kind of divergence. However, 
we can understand qualitatively why the probability of creating a pair with 
larger and larger n d and then with larger and larger energy decreases so 
slowly. In fact, the quantum number rid is also connected with the radius 
p± of the helix along which a classical electron performs its motion (see Fig. 
I2.ip . In particular, it can be shown that p\ ~ n d [73J. From this point 
of view, while creating, for example, an electron with larger and larger n d 
needs an amount of energy proportional to [see Eq. (|2.1.19a|) ]. the 

electron wave function extends over a volume that also increases with n d , in 
such a way that the magnetic energy available for the electron creation also 
increases with n d - 
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I want to show here with some detail why the function f(k,t)dk given in 
Eq. (|4.1.2|) can be interpreted as the mean number of electrons (positrons) 
per unit volume present at time t with a longitudinal momentum between 
k and k + dk. 

The argument I propose is closely related to the treatment given in 
in dealing with multiple soft photon production. Firstly, I remind that since 
I am considering the production from vacuum of electron-positron pairs up 
to first order in the adiabatic perturbation theory, the vacuum persistence 
probability is one [see the discussion below Eq. (|2.2.1U|) ]. Also, following 
the adiabatic perturbation theory, f(k,t)dk represents the probability per 
unit volume that a pair is present at time t with the electron with a longitu- 
dinal momentum between k and k + dk and the positron with a longitudinal 
momentum between —k and —k — dk. Since there is a sharp correlation be- 
tween the electron and the positron quantum numbers, f(k,t)dk is also the 
probability per unit volume that an electron is present at time t with a longi- 
tudinal momentum between k and k + dk or, symmetrically, that a positron 
is present at time t with a longitudinal momentum between k and k + dk 
[note from Eq. (I4.1.2P that f(k,t) = f(—k,t)]. Then, the total probability 
that an electron is present at time t is given by the integral J dkdV f(k,t) 
where dV = L z dA± = L z 7rdR± M . Since no interaction is introduced among 
the particles produced, I would conclude that the probability of finding at 
time t two electrons with a longitudinal momenta between k and k + dk 
and between k' and k' + dk' respectively is f(k,t)f(k',t)dkdk'dVdV' and 
so on. But, really, I must take into account the overall conservation of the 
probability, then, by summing all the terms and by taking into account that 
the particles are indistinguishable, I get the normalization factor 



N(t) 
that is 



1 + J dkdVf(k,t) + ^ J dkdk'dVdV f(k,t)f(k' ,t) + . . . 



N(t) = exp 



dkdVf(k,t) 



: 1, 
(B.l) 

(B.2) 



In conclusion, the total probability that r electrons are present at time t is 
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given by 



rl 



dkdVf(k,t) 



1 r 

exp 



dkdVf(k,t) 



(B.3) 



which is a Poissonian distribution in the number of electrons present. This 
allows me to interpret the quantity J dkdV f(k,t) as the mean number of 
electrons present at time t and then f(k,t)dk as the mean number of elec- 
trons per unit volume present at time t with a longitudinal momentum 
between k and k + dk. 

In the previous description the pair is treated as an effective boson. This 
is allowed when the production rate is low, so there is a little chance of having 
two electrons in the same cell of the phase space, where the Pauli princi- 
ple would play the dominant role (the two electrons would be necessarily 
correlated). In order to check the consistency of this treatment I compare 
the mean number of electrons produced with the number of available quan- 
tum states in the same conditions. The number of the electrons produced 
per unit volume and unit longitudinal momentum is given, according to the 
previous discussion, by the function f(k,t), while the number of available 
quantum states per unit volume and unit longitudinal momentum is [106J 

n(k) = -l^t = ^ (B.4) 
y J 2tt 2vr 4tt 2 v ; 

where I have not considered the spin factor because in the present problem 
the electrons created have fixed spin direction [see Eq. (|3.2.7p ]. In this way, 
the ratio f(k,t)/n(k) is always less than 

1 /„ _ B r x 2 



T2{ nR ^'tj ■ (B ' 5) 

But, by following the same technique used to obtain the strong inequality 
(I3.2.43p . one sees that the quantity UR±mB^ / B cr must be much less than 
one in order that the first-order adiabatic perturbation theory can be safely 
applied, then all the previous description is coherent. 
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In this Appendix I want to calculate the matrix element uQv [see Eq. (|4.1.5p ] 
of the pair annihilation into two photons process by using not the electron 
propagator in vacuum as in the main text but the so-called Schwinger prop- 
agator [5] that is the electron propagator in the presence of a constant and 
uniform magnetic field that I will indicate as B^ = (0, 0, Bf). The structure 
of the Schwinger propagator is, in general, very complicated but it simplifies 
in the strong field approximation (B+ 3> B cr ) in which I am working. In 
this case the propagator is obtained as a sum only over all the TGSs of the 
electron and of the positron and, working in the symmetric gauge (I2.1,42p . 
it is given by [96J: 



G (t, r, t , r ) = exp 



eB 



i——{xy' - x'y) 



G'{t 



t'r 



r') (C.l) 



where (t, r) and (£', r') are two fourpoints and where 

eBi 



G'(T,R) =^exp 



-(X 2 + Y 2 



dWdK Wj° - Kj 3 + 



(2tt) 2 W 2 - K 2 - m 2 



n_ exp[-i(WT - KZ)} 
(C.2) 



with R = (X, Y, Z) and with ITL = (1 — <7 2 )/2 is the spin-down projector. 

If I use the usual notation for the photon field [see for example [97j], 
the transition amplitude of the pair annihilation into two photons can be 
written as 

^n'^.A'^^'i^q') = -4vra em J dtdrdt'dr' 

x V,fc'( r ) exp(-ie't)7 Q iG'(t, r, t' , r )^u' n>k (r' ) exp(-iet') 
■ exp[t{ujt — q • rjj — j^=- exp[i(uj t — q • r JJ 



. (e q ,A)a r., ,i ,^ (2q',A')/3 r./ /, / Nil 

H — __ exp[z(u;t — q • r )J — exp[a(u; t — q • rjj > 

y2Vuj' J 



v 7 ^ 



(C.3) 
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where the initial electron and positron are assumed to be in the TGSs 
u' nk (r') with energy e = e(k) = \/m 2 + k 2 and v' nk (r') with energy e' = 

e(k') = V m 2 + k' 2 respectively and where the final photon states are those 
described in the main text before Eq. (|4.1.5p (V is the quantization volume 
and the limit of large V is understood). In the following, I will calculate 
only the first amplitude in Eq. (|C.3j) corresponding to the creation of the 
photon (q, A) in (i, r) and of the photon (q',A') in (t',r') and I will call it 

S'n n' A A'(^' I; l')- The second amplitude can be calculated in an analo- 
gous way. 

Now, from Eqs. (|2.1.35|) one notes that the difference of the TGSs with 
those of the electrons (positrons) freely propagating along the z axis with 
spin down (up) lies only in the dependence on the transverse coordinates 
[97] . By inserting Eqs. (|2.1.35f) and the Schwinger propagator in Eq. (|C.3j) . 
I observe that the integrals on the time variables and on the longitudinal 
variables can be performed exactly. In fact, they give two 5 functions that 
allow to calculate the integrals on W and on K in the Schwinger propagator 
and that guarantee the energy and the longitudinal momentum conservation. 
In this way, the amplitude S'^, A x ,(k, k' , q, q') can be written as 



s: 



'(1) , 

n,n',X,X' [ 



c,/c',q,q) 



2ira e 
L z VVlolo 
x (2vr) 2 (5(A; + k! — q z — q' z )5(e + e' 



(C.4) 



to 



where L z is the length of the quantization volume in the z direction. In this 
expression I defined the transition matrix element 



(e + m)(e' + m) , 

^n/fe, q y ,q x ,q y ) 



x 



2e2e' 
k' 



0-1 



x n_7 /3 (e q / ;A /) /3 



e' + m 

7°(o; — e') — 7 3 (g z — k') +m 



{uj — e') 2 — (q z — k') 2 — m 2 

( \ 
1 


k 

\ e + m/ 



(C.5) 
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where the function 



1 



1 



d£dqdgdr{(g - ir/) n (£ + irj) n ' 



7T 2 yjnln'l 

x exp - irf) - ir](£' + it/)] 



x exp 



eB. 



■(Qx€ + QyV + qW + QyV'^ 



(C.6) 



depends only on the transverse momentum variables q x , q y , q' x and q' y . 

The transition amplitude (lC.4j) is to be compared with the analogous 
calculated in the vacuum [see for example Eq. (25.9) in [95J] 



S$ iV ,(k,k',q ) q') 



27ra P 



7^,A,v(k,k',q,q') 



F 2 vW 

x (2^) 4 <5(k + k' - q - q')5(e + e'-u-u') 



(C.7) 



with, for an electron moving along the z axis with spin down and a positron 
moving along the same axis with spin up, 

A A ,(k J k / ,q,q') = Qi^(fc,A/,q,q') 



(e + m)(e' + m) 
2e2e' 



A:' 



e' + m 



- 1 



X T 7^ ~ 2 ~2 T2 ~ ~2 ~2 TK^,\>)p 



(u - a') 



q 2 



(q z - k 



( 



\ 



k 



\ e + m/ 



(Ci 



Now, suppose that the transverse momenta of the outgoing photons are 
much smaller than the electron mass m. In this approximation the matrix 
element in the vacuum has a weak dependence on them. In the same way, 
since m 2 /e = B cr <C Bj- then q x <C *JeBf, q y <C >JeBf, q' x -C v/eB/ and 
q'y <C \J eBj and from Eq. (|C.6p we see that in this case N' n n ,(q x , q y , q' x , q' y ) 

[and then ^, y(k, k', q, q')] does not depend on the magnetic field. In 
particular, 

N n,n'(qx,qy,q'x,qy) - <W (C9) 

that is, the transverse structure of the electron and positron states does not 
influence the matrix element Q'zli \ y(k, k', q, q'). Also, by comparing Eqs. 
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(IC.5j) and (1C.8|) . one sees that in this approximation the matrix element in 
the presence of the magnetic field B^ and the matrix element in the vacuum 
have a very similar structure. These observations give me the possibility to 
conclude that one is allowed to use the approximated treatment leading to 
Eq. lji4-l-5\ ) when the photons have small transverse momenta compared with 
the electron mass m. This situation is verified trivially when the electron and 
positron momenta are small compared with m but also when the incoming 
particles are in the ultrarelativistic regime: in fact, in this case, the photon 
production shows a pronounced peak in forward and backward directions 
|112j . Instead, when the transverse momenta of the photons are not small, 
Eqs. (|C.6p and (|0.8|) suggest that the corrections to the matrix element 
in the presence of the magnetic field B^ are proportional to the quantities 
q x j \fe~Bf, q y j yJeBf, q' x / yj eBf and q' y / yj eBf while those to the matrix 
element in the vacuum are proportional to q x /m and q y /m. In particular, by 
performing the Fourier transform of the Schwinger propagator one sees that 
it contains exponential terms in the squared transverse momenta and this 
implies that, by using the "vacuum" quantities to calculate the differential 
cross section da(k, k' , uj)/dw in Eq. (|4.1.3[) . one overestimates the number 
of photons emitted with large transverse momenta. 
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I want to give here the explicit expression of the zero-order electron and 
positron TGSs and the zero-order electron and positron states corresponding 
to the first-excited Landau levels. I remind that I used them to calculate 
the transition matrix elements (|5.2,72|) . These states can be easily obtained 
by substituting Eq. (|2.1.25l) in Eqs. (15.2.441) and (|5.2.48|) 
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U' m (r) 

U l,k,-hx n \. T ) 



V' {0) (r) 



V' {0) (r) 



4 0) + Vdim 



2£ 



(0) 







V Vi.^ +VSim \ kQ' hxo (x,y) J J 



' 9t 



exp(ikz) 



(D.ld) 



4 0) + V9i m 



(. 



(0) 



i f -J&m&v) V 







exp(— i/cz) 



(D.le) 



4 0) + 



2.5 



(0) 







exp(— z/cz) 



1, 



(D.lf) 



In these expressions I have used the definitions (I5.2,4ip and (|5.2.49p of the 
twodimensional spinors ^j(r) and Xj(r) and the definitions (|5.2.54j) and 

115. 2.58ft of the energies ei°^ and 
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In this Appendix I want to show that the set of spinors U' nd ^ ng a (E;r) and 
V^ d>n )0 .(E;r) with E > is complete. Since it is equivalent, but math- 
ematically easier, I will show that the set of spinors U' nd ^ n jCT (E; r) with 
— oo < E < oo is complete. In practice I have to show that 

/°° 00 _j_ 1 

dE E E K,n 9)CT (^r)]J^ J n giCT (^rO] 6 =g i V(r-rO 
"°° n d ,n g =0cr=-l 

(E.l) 

where a,b = 1, . . . , 4 are two spinor indices and where the definition (|5.3.8[) 
of the scalar product between two spinors has been taken into account. By 
using the general expression (|5.3.27j) of the spinors U' ni rig {E\ r) and the fact 
that the projectors P± are two real and symmetric matrices [see Eq. (15.3. 17j) 
and remind that I work in the Dirac representation of the 7 matrices] , I can 
write the previous equation as 



J —oo „ . „ „ 7r 



x n d ,n g ,cr 



n d Z )\ac 



X [ ,=; n d ,n 9 ,o-( X )y)] c [' = 'n d ,n 9 ,o-( :r > D )] d 

x [P-K l/2 _ iu (Ak ni z') + P + K 1/2+iy (Ak nd z')\ db = 5 a , b 5(r - r') 

(E.2) 

where v = IEtq and where the summation on the spinor indices is under- 
stood. Now, by using the integral representation 



1 f 00 

K x (0 = , / \ I cfecos(£sinhs)cosh(As) if |Re(A)| < 1 and f > 

cos(A7r/2) J Q 

(E.3) 

of the modified Bessel functions, it can easily be shown that 

r ^ 8 ^ co f(^) ^ 1/2Tw(4 fc radZ )K 1/2±w (4fc red zO = 25(z - z>). (E.4) 

J -co 
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In the same way it can be seen that the integrals 

_o ' Kl/2±iv (4fcn d z)Kl /2±iu {<±K d Z') 



8fc nd C OSh(7T^) 



7T" 



(E.5) 



vanish by using the integral representation (IE.3j) for one of the Bessel func- 
tions and the following 

1 r°° 

K x (0 = — 7T , . / dssin(£sinhs)sinh(As) if Re(A)| < 1 and £ > 
sm(A7r/2) J 

(E.6) 



for the other. With these results Eq. (|E.2j) is true if 



n d ,n g ,a 



+ [P + ] ac [S^fo y)] c 2/')] d [P+U ( E - 7 ) 



-5 a>b 5(x - x')5(y - y'). 



Now, by using the expression (|5.3.24p of the spinors E' nd>n j(7 (x,y) it can be 

seen that the matrix M(x,y,x',y') = Y,n d ,n g ,a E n d ,n g ,a( x >y"> E ™d,n g A x ',y') 
has the general structure 



M(x,y,x',y') 

/E(x,y,x',y') A(x,y,x',y') -B(x, y, x' , y') \ 

D(x,y,x',y') E(x,y,x',y') C(x,y,x',y') 

B(x,y,x',y') E(x,y,x',y') A(x,y,x',y') 

-C(x,y,x',y') D(x,y,x',y') E(x,y,x' ,y')J 

(E.8) 



V 



with 



A(x,y,x',y' 



B(x,y,x',y' 



C(x,y,x',y' 



^ [^.rig.ffl 1 ) U)] i [ l=, nd,n ff ,o-( a; > V )] 2 
i,n g ,cr 

["n ( i,n 9 ,0-( a '> 2/)]s [^Tid.ngjO-^ '2/ )] 45 

i,n g ,cr 

^ Kd,n 9 ,ff( x ) 2/)] 3 [ ,= 'ra d ,n 9 ,o-( x ) 2/ )] j 

[^Hd.ng.ffl 1 ! 2/)] l ["nd^gjO-C 35 > 2/ )] 3' 
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■ 9 ,cr 
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n d ,n g ,a 



(E.9a) 



(E.9b) 



(E.9c) 



n d ,n g ,cr 
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D(x,y,x',y' ) = ^ [Z' n ^ ng ^(x, y)] 2 [E'* djng ^(x' , y')] 1 

nd,n„,<T 

, (E.9d) 

n d ,n g ,a 

E(x,y,x',y') = ^ [H / ridinfl>CT (x,y)] a [H^ )n9i(T (x / ,y / )] a a = l,...,4. 

(E.9e) 

If one performs the products among the matrices P± and M(x,y,x' ,y') in 
Eq. (|E.7h . it can be shown that 

[P-\ac[M(^x\y')] cd [PA db 

+ [P + ] QC [M(x,y,x',y')] cd [P + ] db = 6 a , b E a (x, y, x' , y'). 



In this way, it is evident that the exact expressions of the four functions 
A(x, y, x' , y'), . . . , D(x, y, x' , y') are not needed. Finally, by using the com- 
pleteness of the set of spinors z>' ndtn a (x, y) one sees that 



E a (x, y, x', y') = -S(x - x')5(y - y') a = l,...,4 (E.ll) 

and then that Eq. (IE.7j) is, actually, an identity. 
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In this Appendix I will impose that the spinor (|5.3.36j) is continuous at z = b 
and that its norm is unit. As a result, I will discretize the energies E and 
determine the two factors Nj and A r j > ^ with J = {n, rid, n g , a} appearing 
in Eq. (|5.3.36j) . The continuity condition is satisfied if 

^k nd b) = A^^l/2+2^^^ {±K d b) (F.l) 

where I pointed out that N^n d and N^„ d can not depend on n g and a 
and that the energies depend on a new integer quantum number n. Since 
k nd b « 1, I can use the approximated expressions of the modified Bessel 
functions near the origin [89j 

K^E^adknA ~ F(1/2 + 2 2 i£ "~ ro) (2fc„»-^-^-.,- C (F.2b) 
to write Eq. (|F. 1[) in the form 

-m ( - + 2iE nind r G ) T 2 ( - + 2zE n , nd r G ) = 1 (F.3) 

^n,n d ^ ' ^ ' 

where the property T(z+1) = zT(z) has been used. By equating the modulus 
and the phase of the left and right hand sides of Eq. (|F.3j) , I obtain the two 
real conditions 



{<) jk_ VI + {^ nd r G f (>) 

8k nd b co S h(^E nind r G ) w [ ■ aJ 



arctan(4^^, G ) + arg (r fl + 2rE n , nd r G )) = 

(F.4b) 



where n = 0, ±1, . . . and where the following property of the V function has 
been used |89|: 
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2 

with (GK. (F.5) 



COSh(-7T^) 
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The condition (IF,4bj) determines the allowed discrete energies while, in or- 
der to determine N^n d , I have to require that the following normalization 
condition holds [see the expression (15.3.8)) of the scalar product]: 

64 J dvU'l ni ^ a {v)U' nM ^(v) = 1. (F.6) 

By exploiting the orthonormalization condition (I5.3.25P on the transverse 
spinors S' n ^ n t(T (x, y), it can be seen that the previous condition is equivalent 
to require that 



64 



k nd r G cosh(2£' n 

,n d r G ) 

4^2 



dz 



l/2+2i£„, n ,r G \^n d 



dz 



K 



l/2+2iE n , n .r G 



= 1. 

(F.7) 



By using the approximated expressions (|F.2ap calculated in 4:k nd z, the first 
integral gives 



dz 



hi 



2+2i£ n , n ,r G (^k nd z) 



2 1 cosh(2nE n , nd r G ) {k nd bf 



k nd 1 + (4£ , „ i „ d r G ) 2 7T 
The second integral can be evaluated by using the following identity 

2 



(F.8) 



dz 



Kl/2+2tE n , nd r G ^ k n d z) 



lim 



d' 

oo 



dz(4k nd zY K 1/2+2iE G (4k nd z) 



(F.9) 



dz{Ak nd z) e K x /%\ 



-2iE n 



(4kn d z) 



The first integral on the right hand side of this equation is a particular case 
of the general formula |107| 



dss~ p K x (as)K fl (bs) 
(\ - 

x r 

x F 



r(i - p) 

1 — p — X + p\ ^ f 1 — p + A — p 



1 — p — A — p, 



1- p + A + p 1- p- A + p b 2 
2 ' 2 



(F.10) 



where a, b, p, A and p are complex numbers such that Re(a + b) > and 
Re(p) < 1 — |Re(A)| — |Re(//)| and where F(r, s;u;z) is the hypergeometric 
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function. Instead, in the second integral on the right hand side of Eq. 
(|F,9p the approximated expression (|F.2b|) of the modified Bessel function 
calculated in 4k nd z can be used. Obviously, even if these integrals are both 
diverging in the limit e — > 0, their divergences must cancel each other because 
the left hand side of Eq. (1F.9D is finite. In fact, 



lim 



dz(4k nd z) e K 1/2+2iEntnd r G (4k nd z 
b 

dz(4k nd z) e K 1/2+2iEnndrG (Ak nd z) 



1 



lim 

Ak nd e^O 



-2+e 



r(i + e) 



2 + e 



1 + e 



+ 2iE„ n ,r G 
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2e 



r ( - + 2iE n , nd r G 



r ( - + 2iE n>nd r G 



lim 



27 ~ ^ exp ( elog ( Akn d b )) 



8k 



>i d 



r ( - + 2i£ n , nd r G 



log (4/c nd 6) 



(F.ll) 



where I used the property r(e/2) = 2r(l + e/2)/e. Finally, by exploiting 
Eq. (|R5|> . then 



If 



l/2+2j£; n ,„ d r G (4fen d 2) 



7T log(4fc Wd 6) 
8fc n<j cosh(2£; r! , in(1 r G 



(F.12) 



and, by substituting Eqs. (|R4al) . (|R8|) and (lF~T2"1) in Eq. (|R7l) . I obtain 
the following expression of N^} d 



47T 



•n,n d 



n d 



r G (l-81og(4fc nd 6)) 



(F.13) 



Since, at the end of the calculations the limit b — > will be performed, I 
give the expression of in this limit: 



7T 



2r G log {k nd b) 



k nd b -> 0. 



(F.14) 



In the same limit an easy expression of the density of the energy levels 
Q(E n ,n d ) can De obtained. In fact, this quantity is defined as 



Q{E n , 



dn 
dE n 



(F.15) 
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Now, if k nd b -> then Eq. p^Ibl) gives simply 

E n .n d = ^~ log" 1 (k na b) k Ud b -> (F.16) 

and the density of the energy levels does not depend on the energy itself: 
q {E n>nd ) = g nd = ^ log (k nd b) k nd b 0. (F.17) 

Finally, with this definition the normalization factor can be written in 
the limit k nd b — > simply as 

N$ = -±= k nd b - 0. (F.18) 
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